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RADAR  SCENE  GENERATION  FOR  TACTICAL  DECISION  AIDS 


This  work  was  supported  by  the  U.S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion,  Vicksburg,  Miss,  as  part  of  the  Air  Land  Battle  Field  Environment  (ALBE)  thrust 
program  with  point  of  contact  Dr.  Richard  Weiss.  We  ha’.e  published  67  refereed  jour¬ 
nal  and  conference  papers  for  the  research  on  the  Radar  Scene  Generation  For  Tactical 
Decision  Aids. 

EMSARS  (Electromagnetic  Model  of  Scattering  Applied  to  Remote  Sensing)  is  an 
electromagnetic  background  scattering  model  developed  at  the  Massachusetts  Institute  of 
Technology.  The  physical  characterization  of  the  terrain  background  is  in  terms  of  layered 
continuous  random  medium  and  discrete  scatterer  models.  Various  theories  and  mathe¬ 
matical  techniques  such  as  the  analytical  Wave  Theory  with  the  Born  Approximation  and 
the  Radiative  Transfer  Theory  are  used  to  solve  the  integral  equations  characterizing  the 
electromagnetic  scattering  effects.  We  recently  presented  applications  of  the  EMSARS 
model  to  millimeter  wave  scattering  by  vegetation.  Tree  is  modeled  as  discrete  scatterers 
with  multiple  size  for  the  trunk,  the  branches  and  leaves,  respectively.  The  leaves  of  conif¬ 
erous  trees,  the  branches,  and  the  trunks  are  modeled  as  circular  cylinders.  For  deciduous 
trees,  the  leaves  are  modeled  as  thin  circular  disks.  The  physical  optics  approximation 
is  used  to  calculate  the  millimeter  wave  returns  from  circular  disks  while  finite  cylinder 
approximation  is  used  for  circulsu’  cylinders.  Then,  the  branching  structure  and  their  ori¬ 
entation  distributions  are  applied  to  obtain  the  radar  returns  from  vegetation.  The  trunk 
is  model  as  a  circular  cylinde '  with  random  surface  roughness  where  geometrical  optics 
approximation  is  used  to  calculate  the  radar  return  from  the  surface  of  trunk.  Except  for 
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the  case  in  which  the  ground  is  illuminated,  we  only  consider  single  scattering  returns. 
Transmission  loss  factors,  due  to  absorption  and  scattering,  are  included  when  computing 
the  contribution  &om  interior  branches  and  leaves.  Theoretical  predictions  are  compared 
with  the  measurement  data  collected  in  Bend,  Oregon  by  NVEOD  using  a  94GHz  radar. 


Polarimetric  radar  calibration  algorithms  using  a  combination  of  point  targets  and 
reciprocal  distributed  targets  are  developed.  From  distributed  target,  a  matrix  equation 
is  derived  which  can  be  converted  into  an  equivalent  point  target  response.  The  equiva¬ 
lent  point  target  corresponds  physically  to  a  90-degree  polarization  rotator  which  is  non¬ 
reciprocal  and  rotation-invariant.  Due  to  this  equivalent-point-target  interpretation,  the 
cases  of  polarimetric  calibration  using  two  point  targets  and  one  distributed  target  then 
reduce  to  those  using  three  point  targets,  which  have  been  solved  in  previous  research. 
Regarding  the  calibration  using  one  point  target  and  one  reciprocal  distributed  target, 
two  cases  are  analyzed  with  the  point  target  being  a  sphere  (trihedral  reflector)  or  a  po¬ 
larimetric  active  radar  calibrator  (PARC).  For  both  cases,  the  general  solutions  of  system 
distortion  matrices  are  written  as  a  product  of  a  particular  solution  and  a  matrix  with  one 
free  psirameter,  and  then  an  additionrd  assumption  about  the  distributed  target  is  made 
to  determine  the  free  parameter.  For  the  trihedral-reflector  case,,  when  the  particular  so¬ 
lution  is  applied  for  calibration,  the  measured  polarimetric  data  can  be  crJibrated  to  the 
level  that  only  rotation  errors  remain,  and  once  azimuthal  symmetry  is  assumed  for  the 
distributed  target,  an  iterative  scheme  is  devised  to  solve  the  rotation  parameter.  For  the 
PARC  case,  the  residual  error  is  like  the  channel  imbalance  after  the  data  are  calibrated  by 
the  particular  solution,  and  the  free  parameter  can  be  determined  by  knowing  one  ratio  of 
two  covariance  matrix  elements  of  the  distributed  target.  Numerical  results  were  simuJated  ^ 
to  demonstrate  the  validity  of  the  algorithms  developed. 
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To  develop  an  algorithm  for  retrieving  forest  biomass  using  SAR  data,  a  study  is 
performed,  which  includes:  1)  observations  on  the  relation  between  backscatter  and  forest 
biomass,  2)  interpretation  of  the  observations  by  modeUing  and  3)  simulation  of  radar 
response  for  a  wide  range  of  forest  and  ground  conditions.  We  have  obtained  preliminaiy 
results  of  P-band  backscatter  simulations  as  functions  of  radar  incidence  angle,  ground 
slope,  undergrowth  vegetation  state,  vegetation  moisture  content,  and  crown  structure. 

A  branching  model  is  proposed  for  the  remote  sensing  of  vegetation.  The  frequency 
and  angular  responses  of  a  two-scale  cylinder  cluster  are  calculated  to  demonstrate  the  sig¬ 
nificance  of  vegetation  architecture.  The  results  indicate  that  it  is  necessary  for  theoretical 
models  to  take  into  account  the  architecture  of  vegetation  which  plays  an  important  role 
in  determining  the  observed  coherent  effects.  A  two-scale  branching  model  is  implemented 
for  soybean  with  its  internal  structure  and  the  resulting  clustering  effects  considered.  Fur¬ 
thermore,  at  the  scale  of  soybean  fields,  the  relative  location  of  soybean  plants  is  described 
by  a  pair  of  distribution  function.  The  polarimetric  backscattering  coefficients  are  obtained 
in  terms  of  the  scattering  properties  of  soybean  plants  and  the  pair  distribution  function. 
Theoretic^ll  backscattering  coefficients  evaluated  using  the  hole- correction  pair  distribution 
are  in  good  agreement  with  extensive  data  collected  from  soybean  fields.  Compared  with 
the  independent-scatterer  pair  distribution,  it  is  found  that  the  hole- correction  approxima¬ 
tion,  which  prevents  two  soybean  plants  from  overlapping  each  other,  is  more  realistic  and 
improves  the  agreement  between  the  model  calculation  and  experimental  data  near  normal 
incidence.  Extension  to  a  multiscale  branching  model  can  be  achieved  by  recursion. 
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Historically,  the  volume  scattering  effects  stimulated  the  development  of  the  contin¬ 
uous  random  medium  model  and  the  random  discrete  scattering  model  for  the  description 
of  the  media.  Theoretical  treatments  were  developed  along  two  different  paths.  Invoking 
the  principle  of  reciprocity,  the  wave  theory  based  on  Maxwell’s  equations  has  been  used 
to  calculate  the  emissivity.  The  other  approach  was  to  start  with  the  radiative  transfer 
equations  and  solved  for  the  brightness  temperatures  directly.  Attempts  have  been  made 
to  derive  the  radiative  transfer  theory  from  the  wave  theory.  At  the  setme  time,  both  the¬ 
oretical  approaches  have  been  used  to  calculate  the  radiometric  emissions  and  to  interpret 
experimentally  measured  data.  The  successful  interpretation  of  the  Cosmos  243  data  was 
perhaps  the  first  most  important  step  towards  a  serious  development  of  the  continuous 
random  medium  model  to  account  for  the  volume  scattering  effects  of  snow  ice  fields. 
Subsequent  interpretation  of  measurement  results  from  snow  field  with  both  passive  ra¬ 
diometers  and  active  radar  systems  established  a  unique  position  for  its  description  of  earth 
terrain  media.  Recent  efforts  in  classifying  sea  ice  with  correlation  function  characteriza¬ 
tion  are  demonstration  of  acceptance  of  this  model.  Future  inverse  scattering  developments 
will  perhaps  rely  heavily  on  this  model.  The  emissivity  calcxUations  for  periodiccd  rough 
surfaces  are  also  performed  and  checked  against  measured  data.  Recent  development  in 
polarimetric  active  remote  sensing  with  synthetic  aperture  radar  has  created  significant 
theoretic2il  results  and  practical  applications.  In  passive  remote  sensing,  the  third  and  the 
fourth  Stokes  parameters  for  earth  remote  sensing  have  not  received  much  attention  in 
the  past  partly  due  to  the  expected  small  values  any  measurement  can  yield.  We  have 
made  initial  calculations  and  experimental  measurements  to  show  that  at  least  the  third 
Stokes  parameter  can  give  appreciable  number  in  both  theoretical  prediction  and  actu2Ll 
experimentation  in  a  plowed  field.  Such  results  may  have  practical  implications  in  measur¬ 
ing  wind  directions  in  ocean  waves,  for  instance.  We  believe  polarimetric  passive  remote 
sensing  is  a  viable  field  which  should  be  explored  in  light  of  its  potential  applications  that 
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may  derived  from  the  full  Stokes  vector  instead  of  its  first  two  parameters,  i.e.,ters,  i.e., 
the  horizontally  and  vertically  polarized  brightness  temperature  components. 

A  multivariate  K-distribution,  well  supported  by  experimental  data,  is  proposed  to 
model  the  statistics  of  fully  polarimetric  radar  clutter  of  earth  terrain.  In  this  approach, 
correlated  polarizations  of  backscattered  radar  returns  are  chciracterized  by  a  covariance 
matrix  and  homogeneity  of  terrain  scatterers  is  characterized  by  a  parameter  a.  As  com¬ 
pared  with  C-,  L-,  and  P-  band  polarimetric  SAR  image  simultaneously  measured  by  Jet 
Propulsion  Laboratory  (JPL)  on  Mt.  Shasta,  it  is  found  that  a  appears  to  decrease  from 
C-  to  P-band  for  both  the  forest  and  burned  areas. 

The  concept  of  polarimetry  in  active  remote  sensing  is  extended  to  passive  remote 
sensing.  The  potential  use  of  the  third  and  fourth  Stokes  parameters  U  and  V,  which 
play  an  important  role  in  polarimetric  active  remote  sensing,  is  demonstrated  for  passive 
remote  sensing.  It  is  shown  that,  by  the  use  of  the  reciprocity  principle,  the  polarimetric 
parameters  of  p2issive  remote  sensing  can  be  obtained  through  the  solution  of  the  associated 
direct  scattering  problem.  These  ideas  are  applied  to  study  polarimetric  passive  remote 
sensing  of  periodic  surfaces.  The  solution  of  the  direct  scattering  problem  is  obtained  by 
an  integral  equation  formulation  which  involves  evaluation  of  periodic  Green’s  functions 
and  normal  derivative  of  those  on  the  surface.  Rapid  evaluation  of  the  slowly  convergent 
series  associated  with  these  functions  is  observed  to  be  critical  for  the  fe^lsibility  of  the 
method.  New  formxilas,  which  are  rapidly  convergent,  are  derived  for  the  calculation  of 
these  series.  The  study  has  shown  that  the  brightness  temperature  of  the  Stokes  parameter 
U  can  be  significant  in  passive  remote  sensing.  Values  as  high  as  50  K  are  observed  for 
certain  configurations. 
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Classification  of  sea  ice  types  using  polarimetric  radar  is  an  area  of  considerable 
current  interest  and  research.  In  this  paper,  a  new  clzissification  technique  is  applied 
to  determine  sea  ice  types  in  polarimetric  and  multi-frequency  SAR  images,  utilizing  an 
unsupervised  neural  network  to  provide  automatic  classification,  and  employing  an  iterative 
algorithm  to  improve  the  performzmce.  Learning  Vector  Quantization  (LVQ)  is  first  applied 
for  the  unsupervised  classification  of  SAR  images,  and  the  results  are  compared  with 
those  of  the  Migrating  Means  method.  Results  show  that  LVQ  outperforms  the  Migrating 
Means  method,  which  classified  most  of  the  pixels  into  two  classes.  To  further  improve  the 
performance,  an  iterative  algorithm  is  applied  where  the  SAR  image  is  reclassified  using 
the  Maximum  Likelihood  (ML)  classifier.  It  is  shown  that  this  algorithm  converges,  and 
significcintly  improves  the  performance  of  the  unsupervised  LVQ  method  while  preserving 
the  advantages  of  automatic  operation  inherent  in  unsupervised  techniques. 

Monte  Carlo  simulation  of  random  rough  surface  is  carried  out  using  the  Finite- 
Difference  Time-Domain  (FDTD)  method.  The  solution  domain  is  truncated  using  a 
periodic  and  an  absorbing  boundary  condition.  The  advantages  of  using  FDTD  include 
(1)  matrix  inversion  is  eliminated.  (2)  time-domain  response  is  readily  available,  and  (3) 
scattering  intensities  for  a  range  of  frequencies  can  be  obtained  simultaneously. 

A  random  medium  model  hzis  been  applied  to  interpret  the  radar  backscattering  co¬ 
efficients  measured  over  a  soybean  canopy.  The  vegetation  is  modeled  a  random  medium 
containing  non  spherical  scatterers  which  are  soybean  leaves.  The  theoretical  results  are 
in  good  agreement  with  experimental  data.  An  interesting  observation  is  the  high  correla¬ 
tion  between  the  canopy  fractional  volume  and  the  correlation  length  of  the  permittivity 
fluctuations  in  the  horizontal  plane.  The  model  is  then  used  to  simulate  backscattering 
coefficients  of  various  observation  conditions.  The  resiilts  define  the  different  sensitivity 
domains  of  radar  data  to  soil  moisture  and  vegetation  biomass  for  inversion  consideration. 


7 


The  multiple  species  in  the  random  medium  are  considered  as  randomly  oriented  el¬ 
lipsoids  and  described  by  miiltiple  three-dimensional  ellipsoidal  correlation  functions.  The 
variances  and  correlation  lengths  of  the  correlation  functions  characterize  the  fluctuation 
strengths  and  the  physical  geometry  of  the  scatterers,  respectively.  The  strong  permittiv¬ 
ity  fluctuation  theory  is  extended  to  account  for  the  multiple  species  of  ellipsoidal  shape. 
In  the  rJindom  medium,  a  coincidence  of  an  observation  point  with  a  source  point  gives 
rise  to  the  singularity  of  the  dyadic  Green’s  function  which  is  properly  taken  into  account 
with  exclusion  volumes  of  the  scatterers.  Polarimetric  scattering  properties  of  a  remotely 
sensed  medium  are  depicted  with  a  covariance  matrix  whose  elements  are  polarimetric 
scattering  coefficients.  The  medium  has  a  layer  configuration:  the  top  layer  such  as  air 
is  considered  as  homogeneous,  the  middle  layer  such  ais  a  vegetation  canopy  is  random, 
and  the  underlying  layer  such  as  soil  is  a  homogeneous  half  space.  More  random  medium 
layers  cam  ailso  be  introduced  in  the  configuration  to  account  for  weather  effect  such  as  fog 
cover.  The  distorted  Born  approximation  is  then  used  with  the  effective  permittivity  to 
obtain  the  polarimetric  covariance  matrix.  The  resxilt  for  the  cross-polarized  return 
is  non- zero  even  in  the  first  order  approximation.  Due  to  the  non- spherical  shape  and  the 
random  orientation  of  the  scatterers,  the  correlation  coefficient  between  the  HH  and  VV 
returns  has  a  magnitude  different  from  unity  and  a  small  phztse.  The  scattering  coefficients 
are  also  used  to  c^Jculate  the  Mueller  matrix  for  synthesis  of  polarization  signatures.  The 
copolarized  signature  of  the  random  medium  has  a  rather  straight  distortion  track  and  a 
recognizable  pedestal. 
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In  passive  remote  sensing  of  earth  terrain,  the  Radiative  Transfer  (RT)  theory  is 
widely  applied  with  various  types  of  scattering  models  to  interpret  observation.  Part  of  the 
advantages  of  the  RT  theory  can  be  attributed  to  its  simplicity  in  formulating  the  reflection 
and  transmission  in  the  presence  of  boundaries  as  well  as  in  treating  the  scattering  effect 
by  the  incoherent  sum  of  the  scattering  from  each  individual  scatterer.  However,  being 
restricted  by  this  incoherent  nature,  the  RT  theory  cannot  properly  take  into  account 
the  coherent  effect  which  becomes  significant  at  lower  frequency  in  the  case  of  remote 
sensing  of  a  structured  vegetation  canopy.  Also,  in  the  past,  the  volume  scattering  and 
the  surface  scattering  have  been  treated  separately.  A  rigorous  composite  model  is  thus 
required  to  correctly  evaluate  the  volume-surface  interactions  when  the  rough  surface  is 
visible  by  the  radiometer.  We  formulated  the  vector  radiative  transfer  equation  for  passive 
microwave  remote  sensing  of  a  vegetation  canopy  overlying  a  soil  half-space,  and  calculated 
brightness  temperatures  resulting  from  microwave  thermal  emission  based  on  the  first  order 
iterative  solution.  The  rough  surface  boundary  condition  has  been  integrated  into  this 
model  through  the  coupling  matrix  for  the  boundary  between  the  vegetation  layer  and  soil 
such  that  the  interactions  between  volume  scattering  and  surface  scattering  are  rigorously 
included.  A  randomly  distributed  stem-leave  model  is  employed  to  construct  the  phase 
matrix  such  that  the  effects  of  coherence  and  multi-scsde  can  be  properly  accounted  for. 

We  have  applied  the  radiative  transfer  theory  to  vegetation  with  clustered  struc¬ 
tures.  To  take  into  account  vegetation  structure  in  the  radiative  transfer  theory,  the  phase 
matrix  of  a  vegetation  cluster  is  calcxilated  by  incorporating  the  phase  interference  of  scat¬ 
tered  fields  from  every  components.  Subsequently,  the  resulting  phase  matrix  is  used  in  the 
radiative  transfer  equations  to  evaluate  the  polarimetric  backscattering  coefficients  from  a 
layer  of  medium  embedded  with  vegetation  clusters.  Theoretical  results  are  illustrated  for 
various  kinds  of  vegetation  clusters.  It  is  found  that  the  simulated  polarization,  frequency, 
and  angular  responses  carry  significant  information  regarding  the  structure  of  vegetation 
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clusters.  The  theoretical  resxilts  will  also  be  compared  with  the  signatures  observed  in 
multi-irequency  polarimetric  synthetic  aperture  radar  images. 

Classification  of  terrain  cover  using  polarimetric  radar  is  an  area  of  considerable 
current  interest  and  research.  A  number  of  methods  have  been  developed  to  classify  ground 
terrdn  types  from  fully  polarimetric  Synthetic  Aperture  Radar  (SAR)  images,  and  these 
techniques  are  often  grouped  into  supervised  and  unsupervised  approaches.  Supervised 
methods,  including  both  conventional  Maximum  Likelihood  (ML)  and  more  recent  Multi¬ 
layer  Perceptron  classifiers,  have  yielded  higher  accuracy  than  unsupervised  techniques,  but 
suffer  from  the  need  for  human  interaction  to  predetermine  classes  and  training  regions. 
In  contrast,  unsupervised  methods  determine  classes  automatically,  but  generally  show 
limited  ability  to  accurately  divide  terrain  into  natural  classes.  In  this  paper,  a  new 
terrsdn  classification  technique  is  introduced,  utilizing  unsupervised  neural  networks  to 
provide  automatic  classification,  but  employing  an  iterative  algorithm  vhich  overcomes 
the  poor  accuracy  of  other  unsupervised  techniques. 

In  remote  sensing,  the  encountered  geophysical  media  such  as  agricultural  canopy, 
forest,  snow,  or  ice  are  inhomogeneous  and  contain  scatterers  in  a  random  manner.  Further¬ 
more,  weather  conditions  such  as  fog,  mist,  or  snow  cover  can  intervene  the  electromagnetic 
observation  of  the  remotely  sensed  media.  In  the  modelling  of  such  media  accounting  for 
the  weather  effects,  a  multi-layer  random  medium  model  has  been  developed.  The  scatter¬ 
ing  effects  of  the  random  media  are  described  by  three-dimensional  correlation  functions 
with  variances  and  correlation  lengths  corresponding  to  the  fluctuation  strengths  and  the 
physical  geometry  of  the  inhomogeneities,  respectively.  With  proper  consideration  of  the 
dyadic  Green’s  function  and  its  singularities,  the  strong  fluctuation  theory  is  used  to  cal¬ 
culate  the  effective  permittivities  which  account  for  the  modification  of  the  wave  speed 
and  attenuation  in  the  presence  of  the  scatterers.  The  distorted  Born  approximation  is 
then  applied  to  obtain  the  correlations  of  the  scattered  fields.  From  the  correlation  of  the 
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scattered  field,  calculated  are  the  complete  set  of  scattering  coefficients  for  polarimetric 
radar  observation  or  brightness  temperature  in  passive  radiometer  applications. 

Polarimetric  calibration  algorithms  using  combination  of  point  tairgets  and  recipro¬ 
cal  distributed  targets  are  developed.  From  distributed  target  reciprocity,  derived  is  an 
equivalent  point  tairget  response.  Then  the  problem  of  polarimetric  calibration  using  two 
point  targets  and  one  distributed  target  reduces  to  that  using  three  point  targets,  which 
has  been  solved  in  previous  research.  Regarding  the  calibration  using  one  point  target 
and  one  reciprocal  distributed  target,  analyzed  are  two  ceises  with  the  point  target  being  a 
trihedral  reflector  or  a  Polarimetric  Active  Radar  Calibrator  (PARC).  For  both  cases,  the 
general  solutions  of  the  system  distortion  matrices  are  written  as  a  product  of  a  particular 
solution  and  a  matrix  with  one  free  parameter.  For  the  trihedral-reflector  case,  this  free 
parameter  is  determined  by  assuming  azimuthal  symmetry  for  the  distributed  target.  For 
the  PARC  case,  knowing  one  ratio  of  two  covariance  matrix  elements  of  the  distributed 
target  is  required  to  solve  the  free  parameter. 

Finite  difference  techniques  have  been  applied  to  analyze  various  electromagnetic 
problems  in  both  frequency  and  time  domains.  Examples  of  these  applications  include 
scattering  and  radiation,  microwave  and  millimeter  wave  circuits,  and  hyperthermia.  Be¬ 
sides  the  widespread  applications,  the  latest  research  efforts  have  been  directed  at  achieving 
more  accurate  discretization  schemes,  improving  absorbing  boundary  conditions  for  open 
region  problems,  and  efficient  implementations  on  supercomputers. 

We  also  developed  a  Finite-Difference  Time-Domain  (FD-TD)  scheme  on  triangular 
grids.  The  discretization  scheme  is  based  on  the  combination  of  the  finite  difference  and 
control  region  approximations.  The  flexibility  of  the  triangular  grid  is  utilized  to  provide 
more  accurate  target  modeling  capability  than  the  traditional  FD-TD  technique  of  the 
rectangular  grid. 
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Fully  polarimetric  scattering  of  electromagnetic  waves  from  snow  and  sea  ice  is 
studied  with  a  layered  random  medium  model  and  applied  to  interpret  experiment2d  data 
obtained  under  laboratory  controlled  conditions.  The  snow  layer  is  modeled  as  an  isotropic 
random  medium.  The  sea  ice  is  described  as  an  anisotropic  random  medium  due  to  the 
elongated  form  of  brine  inclusions.  The  underl3ring  sea  water  is  considered  as  a  homo¬ 
geneous  half-space.  The  scattering  effects  of  both  random  media  are  described  by  three- 
dimensional  correlation  functions  with  variances  and  correlation  lengths  corresponding  to 
the  fluctuation  strengths  and  the  physical  geometry  of  the  inhomogeneities,  respectively. 
The  strong  fluctuation  theory  is  used  to  calculate  the  effective  permittivities  of  the  random 
media.  The  distorted  Born  approximation  is  then  applied  to  obtain  the  covariance  matrix 
which  describes  the  fully  polarimetric  scattering  properties  of  the  remotely  sensed  media. 

Polarimetric  scattering  measurements  have  been  from  the  simulated  bare  sea  ice 
in  CRREL,  and  the  two-layer  configuration  is  used  to  model  this  experimental  set-up. 
The  distinction  on  the  characteristics  of  the  media  are  investigated  with  the  conventional 
backscattering  coefficients  and  the  complex  correlation  coefficient  p  between  and  cryv- 
For  ice-type  identification,  the  measured  covariance  matrices  are  studied  with  the  model 
to  infer  the  physical  characteristics  pertaining  to  the  different  ice  types.  The  three-layer 
configuration  is  then  used  to  investigate  the  effects  on  fully  polarimetric  radar  returns  from 
snow  covered  sea  ice. 

There  is  also  considerable  interest  in  identifying  and  classifying  ice  types  by  using 
polarimetric  scattering  data.  Due  to  differences  in  structure  and  composition,  ice  of  differ¬ 
ent  types  such  as  frazil,  first-year,  or  multi-year  can  have  different  polarimetric  scattering 
behaviors.  To  study  the  polarimetric  response  of  sea  ice,  the  layered  random  medium 
model  is  used.  In  this  model,  the  sea-ice  layer  is  described  as  an  anisotropic  random 
medium  composed  of  a  host  medium  with  randomly  embedded  inhomogeneities,  such  as 
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elongated  brine  inclusions,  which  can  have  preferred  orientation  direction.  The  under¬ 
lying  sea-water  layer  is  considered  as  a  homogenous  half  space.  The  scattering  effect  of 
the  inhomogeneities  in  the  sea  ice  are  characterized  by  three-dimensional  correlation  func¬ 
tion  with  variance  and  correlation  lengths  respectively  corresponding  to  the  fluctuation 
strength  and  the  physical  geometry  of  the  scatterers.  The  effective  permittivity  of  the  sea 
ice  is  calculated  with  the  strong  fluctuation  theory  and  the  polarimetric  backscattering 
coefficients  are  obtained  under  the  distorted  Born  approximation.  The  distinction  on  the 
characteristics  of  different  ice  types  are  investigated  with  the  conventional  backscattering 
coefficients  and  the  complex  correlation  coefficient  p  between  and  <Tvv-  The  correlation 
coefficient  p  contains  additional  information  on  the  sea-ice  structure  and  can  be  useful  in 
the  identification  of  the  ice  types.  By  relating  to  the  covariance  matrices,  the  model  is  used 
to  explain  the  polarization  signatures  of  different  ice  types.  In  the  case  of  snow-covered  sea 
ice,  the  snow  layer  is  modeled  as  an  isotropic  random  medium  and  the  obtained  solution 
accounts  for  the  effect  of  snow  cover  on  polarimetric  scattering  properties  of  sea  ice. 

The  layered  random  medium  model  is  used  to  investigate  the  ftilly  polarimetric  scat¬ 
tering  of  electromagnetic  waves  from  vegetation.  The  vegetation  canopy  is  modeled  as  an 
anisotropic  random  medium  containing  nonspherical  scatterers  with  preferred  alignment. 
The  underlying  medium  is  considered  as  a  homogeneous  half  space.  The  scattering  effect 
of  the  vegetation  canopy  are  characterized  by  three-dimensional  correlation  functions  with 
variances  and  correlation  lengths  respectively  corresponding  to  the  fluctuation  strengths 
and  the  physical  geometries  of  the  scatterers.  The  strong  fluctuation  theory  is  used  to 
calculate  the  anisotropic  effective  permittivity  tensor  of  the  random  medium  and  the  dis¬ 
torted  Born  approximation  is  then  applied  to  obtain  the  covariance  matrix  which  describes 
the  fully  polarimetric  scattering  properties  of  the  vegetation  field.  This  model  accounts  for 
all  the  interaction  processes  between  the  boundaries  and  the  scatterers  and  includes  all  the 
coherent  effects  due  to  wave  propagation  in  different  directions  such  as  the  constructive 
and  destructive  interferences.  For  a  vegetation  canopy  with  low  attenuation,  the  boundary 
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between  the  vegetation  and  the  underlying  medium  can  give  rise  to  significant  coherent 
effects. 


The  model  is  used  to  interpret  the  measured  data  for  vegetation  field  such  as  rice, 
wheat,  or  soybean  over  water  or  soil.  The  temporal  variation  of  and  aw  of  the  X- 
band  SAR  data  of  rice  fields  shows  a  wide  range  of  responses  at  different  growth  stages. 
From  the  data  of  wheat,  recognizable  chainges  of  the  angular  and  polarization  behaviour  of 
the  backscattering  coefficients  are  observed  at  X-band  before  and  after  the  heading  of  the 
wheat.  For  soybean,  the  data  collected  during  the  growing  season  shows  similar  results 
for  both  h-  and  v-polarizations.  The  observed  effects  on  backscattering  coeffic'ents  of  the 
vegetation  structural  and  moisture  conditions  at  different  growth  stages  can  be  explained 
by  analyzing  the  different  interaction  processes  pointed  out  by  the  model. 

Accurate  calibration  of  polarimetric  radar  systems  is  essential  for  the  polarimetric 
remote  sensing  of  earth  terrain.  A  polarimetric  calibration  algorithm  using  three  arbitrary 
in-scene  reflectors  is  developed.  The  transmitting  and  receiving  ports  of  the  polarimetric 
radar  are  modeled  by  two  unknown  polarization  transfer  matrices.  These  unknown  matri¬ 
ces  are  determined  using  the  the  measured  scattering  matrices  from  the  calibration  targets. 
A  Polarization-Basis  Transformation  technique  is  introduced  to  convert  the  scattering  ma¬ 
trices  of  the  calibration  targets  into  one  of  the  six  sets  of  targets  with  simpler  scattering 
matrices.  Then,  the  solution  to  the  original  problem  can  be  expressed  in  terms  of  the 
solution  obtained  using  the  simpler  scattering  matrices.  The  uniqueness  of  polarimetric 
calibration  using  three  targets  is  addressed  for  all  possible  combinations  of  calibration 
targets.  The  effect  of  misalignment  of  the  calibration  targets  and  the  sensitivity  of  the 
polarimetric  calibration  algorithm  to  the  noise  are  illustrated  by  investigating  several  sets 
of  calibration  targets  in  detail. 
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In  the  interpretation  of  active  and  passive  microwave  remote  sensing  data  from 
earth  terrain,  the  random  medium  model  has  been  shown  to  be  quite  successful.  In  the 
random  medium  model,  a  correlation  function  is  used  to  describe  the  random  permittivity 
fluctuations  with  associated  mean  and  variance.  In  the  past,  the  correlation  functions 
used  were  either  assumed  to  be  of  certain  form  or  calculated  from  cross  sectional  pictures 
of  scattering  media.  We  calculate  the  correlation  function  for  a  random  collection  of 
discrete  scatterers  imbedded  in  a  background  medium  of  constant  permittivity.  Correlation 
functions  axe  first  calculated  for  the  simple  cases  of  the  uniform  distribution  of  scatterers 
and  the  uniform  distribution  with  the  hole  correction.  Then,  the  correlation  function  for 
a  more  realistic  case  is  outained  using  the  Percus-Yevik  pair  distribution  function.  Once 
the  correlation  function  ij  obtaiined,  the  strong  fluctuation  theory  is  used  to  calculate 
the  effective  permittivities.  Then,  the  distorted  Born  approximation  is  used  to  calculate 
the  backscattering  coefficients  from  a  halfspace  configuration.  The  theoretical  results  are 
illustrated  by  comparing  the  effective  permittivities  and  the  backscattering  coefficients 
with  the  results  obtained  with  the  discrete  scatterer  theory. 

Strong  permittivity  fluctuation  theory  has  been  used  to  solve  the  problem  of  scat¬ 
tering  from  a  medium  composed  of  completely  randomly  oriented  scatterers  under  the  low 
frequency  limit.  Based  on  Finkel’berg’s  approach,  Gaussian  statistics  is  not  assumed  for 
the  renormalized  scattering  sources.  The  effective  permittivity  is  obtained  under  the  low 
frequency  limit  and  the  resvdt  is  shown  to  be  isotropic  due  to  no  preferred  direction  in  the 
orientation  of  the  scatterers.  Numerical  results  of  the  effective  permittivity  are  illustrated 
for  oblate  and  prolate  spheroidal  scatterers  and  compared  with  the  results  for  spherical 
scatterers.  The  results  derived  are  shown  to  be  consistent  with  the  discrete  scatterer  the¬ 
ory.  The  effective  permittivity  of  random  medium  embedded  with  nonspherical  scatterers 
shows  a  higher  imaginary  part  than  that  of  spherical  scatterer  case  with  equal  correlation 
volume.  Under  the  distorted  Born  approximation,  the  polsuimetric  covariance  matrix  for 
the  backscattered  electric  field  is  calculated  for  the  half-space  randomly  oriented  scatterers. 
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The  nonspherical  geometry  of  the  scatterers  shows  significant  effects  on  the  cross-polarized 
backscattering  returns  <Tf^^  eind  the  correlation  coefficient  p  between  HH  and  VV  returns. 
The  polairimetric  backscattering  scattering  coefficients  can  provide  useful  information  in 
distinguishing  the  geometry  of  scatterers.  A  polarimetric  radu  calibration  algorithm  us¬ 
ing  three  in-scene  reflectors  is  developed  based  on  the  exact  solution  for  general  target 
choices.  The  transmitting  and  receiving  ports  of  the  polarimetric  radar  are  modeled  by 
two  unknown  polarization  transfer  matrices.  These  transfer  matrices  are  solved  in  terms  of 
measurements  from  three  independent  calibration  targets  with  known  scattering  matrices. 
The  solutions  for  several  sets  of  calibration  targets  with  simple  scattering  matrices  are 
presented  first.  Then,  when  at  least  two  of  the  target  scattering  matrices  can  be  simulta¬ 
neously  diagonalized,  polarimetric  calibration  is  derived  using  the  method  of  simultaneous 
diagonalization  of  two  matrices.  A  transformation  matrix  is  found  to  convert  the  general 
scattering  matrices  into  the  simple  cases,  and  the  problem  is  solved  in  the  transformed 
domain.  The  solution  to  the  original  problem  then  can  be  expressed  in  terms  of  the  solu¬ 
tions  obtained  for  the  simple  scattering  matrices.  All  possible  combinations  of  calibration 
targets  are  discussed  and  the  solutions  are  presented  for  the  cases  that  at  least  two  of  the 
scattering  matrices  can  be  simultaneously  diagonalized. 

Conventional  classification  techniques  for  identification  of  vehicle  types  from  their 
range  profiles,  or  pulse  responses,  have  been  shown  to  be  limited  in  their  practical  ability 
to  distinguish  targets  of  interest.  These  limitations  arise  from  the  need  for  large  signature 
libraries  and  time  consuming  processing  for  profile  matching  algorithms,  and  from  the  sis- 
sumptions  made  toward  the  statistics  of  extracted  features  for  parametric  methods.  To 
overcome  the  practical  constraints  of  existing  techniques,  a  new  method  of  target  recogni¬ 
tion  is  examined  which  utilizes  neural  nets.  The  effectiveness  of  this  neural  net  classifier 
is  demonstrated  with  synthetically  generated  range  profiles  for  two  sets  of  geometries,  ais 
produced  using  RCS  prediction  techniques.  The  first  set  consists  of  three  simple  canonical 
geometries  for  which  RCS  predictions  can  be  done  directly.  For  these  targets,  two  neural 
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net  configurations  are  compared,  and  the  effects  of  varied  aspect  sampling  density  for  the 
training  profiles  and  noise  corruption  in  the  test  profiles  are  demonstrated.  Comparisons 
are  made  between  the  neural  net  classifier  and  several  conventional  techniques  to  deter¬ 
mine  the  relative  performance  and  cost  of  each  algorithm.  A  similar  set  of  comparisons  is 
performed  for  the  second  group  of  targets  consisting  of  more  realistic  air  vehicle  models, 
each  composed  from  a  collection  of  canonical  shapes.  In  both  cases,  the  neural  net  classifier 
is  shown  to  match  or  exceed  the  performance  of  conventional  algorithms  while  offering  a 
more  computationally  efficient  implementation. 

Earth  terrains  are  modeled  by  a  two-layer  configuration  to  investigate  the  polari- 
metric  scattering  properties  of  the  remotely  sensed  media.  The  scattering  layer  is  a  random 
medium  characterized  by  a  three-dimensional  correlation  function  with  correlation  lengths 
and  variances  respectively  related  to  the  scatter  sizes  and  the  permittivity  fluctuation 
strengths.  Based  on  the  wave  theory  with  Born  approximations  carried  to  the  second  or¬ 
der,  this  model  is  applied  to  derive  the  Mueller  and  the  covariance  matrices  which  fully 
describe  the  polarimetric  scattering  characteristics  of  the  media.  Physically,  the  first-  and 
second-order  Born  approximations  account  for  the  single  and  double  scattering  processes. 
For  an  isotropic  scattering  layer,  the  five  depolarization  elements  of  the  covariance  matrix 
are  zero  under  the  first-order  Born  approximation.  For  the  uniaxial  tilted  permittivity  case, 
the  covariance  matrix  does  not  contain  any  zero  elements.  To  account  for  the  randomness 
in  the  azimuthal  growth  direction  of  leaves  in  vegetation,  the  backscattering  coefficients 
are  azimuth2dly  averaged.  In  this  case,  the  covariance  matrix  contains  four  zero  elements 
although  the  tilt  angle  is  not  zero.  Under  the  second-order  Born  approximation,  the  co- 
variance  matrix  is  derived  for  the  isotropic  and  the  uniaxial  untilted  random  permittivity 
configurations.  The  results  show  that  the  covariance  matrix  has  four  zero  elements  and 
a  depolarization  factor  is  obtained  even  for  the  isotropic  case.  To  describe  the  effect  of 
the  random  medium  on  electromagnetic  waves,  the  strong  permittivity  fluctuation  theory, 
which  accounts  for  the  losses  due  to  both  of  the  absorption  and  the  scattering,  is  used  to 
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compute  the  effective  permittivity  of  the  medium.  For  a  mixture  of  two  components,  the 
frequency,  the  correlation  lengths,  the  fractional  volume,  and  the  permittivities  of  the  two 
constituents  are  needed  to  obtain  the  polarimetric  backscattering  coefficients.  Theoretical 
predictions  are  illustrated  by  comparing  the  results  with  experimental  data  for  vegetation 
fields  and  sea  ice. 

The  random  medium  model  with  three-layer  configuration  is  developed  to  study 
fully  polarimetric  scattering  of  electromagnetic  waves  from  geophysical  media.  This  model 
can  account  for  the  effects  on  wave  scattering  due  to  weather,  diurnal  and  seasonal  varia¬ 
tions,  and  atmospheric  conditions  such  as  ice  under  snow,  meadow  under  fog,  and  forest 
under  mist.  The  top  scattering  layer  is  modeled  as  an  isotropic  random  medium  which 
is  chciracterized  by  a  sccdar  permittivity.  The  middle  scattering  layer  is  modeled  as  an 
anisotropic  random  medium  with  a  symmetric  permittivity  tensor  whose  optic  axis  can 
be  tilted  due  to  the  preferred  alignment  of  the  embedded  scatterers.  The  bottom  layer  is 
considered  as  a  homogeneous  half-space.  Volume  scattering  effects  of  both  random  media 
are  described  by  three-dimensional  correlation  functions  with  variances  and  correlation 
lengths  corresponding  to  the  strengths  of  the  permittivity  fluctuations  and  the  physical 
sizes  of  the  inhomogeneities,  respectively.  The  strong  fluctuation  theory  is  used  to  derive 
the  mean  fields  in  the  random  media  under  the  bilocal  approximation  with  singularities 
of  the  dyadic  Green’s  functions  properly  taken  into  account  and  effective  permittivities 
of  the  random  media  are  calciilated  with  two-phase  mixing  formulas.  The  distorted  Born 
approximation  is  then  applied  to  obtain  the  covariance  matrix  which  describes  the  fully 
polarimetric  scattering  properties  of  the  remotely  sensed  media. 
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The  three-layer  configuration  is  first  reduced  to  two-layers  to  observe  fully  polaxi- 
metric  scattering  directly  from  geophysical  media  such  as  snow,  ice,  and  vegetation.  Such 
media  exhibit  reciprocity  eis  experimentally  manifested  in  the  close  proximity  of  the  mea¬ 
sured  backscattering  rad^lr  cross  sections  and  and  theoretically  established  in  the 
random  medium  model  with  symmetric  permittivity  tensors.  The  theory  is  used  to  inves¬ 
tigate  the  signatures  of  isotropic  and  anisotropic  random  media  on  the  complex  correlation 
coefficient  p  between  and  a  function  of  incident  angle.  For  the  isotropic  ran¬ 

dom  medium,  p  has  the  value  of  approximately  1.0.  For  the  untilted  anisotropic  random 
medium,  p  has  complex  values  with  both  the  real  and  imaginary  parts  decreased  as  the 
incident  angle  is  increased.  The  correlation  coefficient  p  is  shown  to  contain  information 
about  the  tilt  of  the  optic  axis  in  the  anisotropic  random  medium.  As  the  tilted  angle  be¬ 
comes  larger,  the  magnitude  of  p  is  maximized  at  a  larger  incident  angle  where  the  phase 
of  p  changes  its  sign.  It  should  be  noted  that  the  tilt  of  the  optic  axis  is  aJso  related  to  the 
nonzero  depolarization  terms  in  the  covariance  matrix. 

The  effects  on  polarimetric  wave  scattering  due  to  the  top  layer  are  identified  by 
comparing  the  three-layer  results  with  those  obtained  from  the  two-layer  configuration. 
The  theory  is  used  to  investigate  the  effects  on  polarimetric  radar  returns  due  to  a  low-loss 
and  a  lossy  dry-snow  layers  covering  a  sheet  of  thick  first-year  sea  ice.  For  the  low-loss 
snow  cover,  both  a/j/j  and  <tvv  are  enhanced  compared  to  those  observed  from  bare  sea  ice. 
Furthermore,  the  boundary  effect  is  manifested  in  the  form  of  the  oscillation  on  and 
Cvv-  The  oscillation  can  also  be  seen  on  the  real  and  imaginary  parts  of  the  correlation 
coefficient  p.  The  magnitude  of  p,  however,  does  not  exhibit  the  oscillation  while  clearly 
retaining  the  same  characteristics  as  observed  directly  from  the  uncovered  sea  ice.  In  con¬ 
trast  to  the  low-loss  case,  the  lossy  top  layer  can  diminish  both  and  <rvv  and  depress 
the  boundary-effect  oscillation.  When  the  thickness  of  the  lossy  top  layer  increases,  the 
behavior  of  the  correlation  coefficient  p  becomes  more  and  more  similar  to  the  isotropic 
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case  signifying  that  the  information  from  the  lower  anisotropic  layer  is  masked.  At  appro¬ 
priate  frequency,  the  fully  polarimetric  volume  scattering  effects  can  reveal  the  information 
attributed  to  the  lower  layer  even  if  it  is  covered  under  another  scattering  layer.  Due  to  the 
physical  base,  the  random  medium  model  renders  the  polarimetric  scattering  information 
useful  in  the  identification,  classification,  and  radar  image  simulation  of  geophysical  media. 

Polarimetric  radar  backscatter  data  observed  with  satellite  and  airborne  synthetic 
aperture  radars  (SAR)  have  demonstrated  potential  applications  in  geologic  mapping  and 
terrain  cover  classification.  Accurate  calibration  of  such  polarimetric  radar  systems  is 
essential  for  polarimetric  remote  sensing  of  earth  terrain.  A  polarimetric  calibration  algo¬ 
rithm  using  three  in-scene  reflectors  has  been  developed  which  will  be  a  useful  tool  in  the 
radar  image  interpretation. 

The  transmitting  and  receiving  ports  of  the  polarimetric  radar  are  modeled  by  two 
unknown  polarization  transfer  matrices.  The  measured  scattering  matrix  is  equal  to  the 
product  of  the  transfer  matrix  of  the  receiving  port,  scattering  matrix  of  the  illuminated 
target,  the  transfer  matrix  of  the  transmitting  port,  and  a  common  phase  factor.  The 
objective  of  polarimetric  radar  calibration  is  to  determine  these  two  unknown  polarization 
transfer  matrices  using  measurements  from  targets  with  known  scattering  matrices. 

The  transfer  matrices  for  the  transmitting  and  receiving  ports  are  solved  in  terms 
of  measurements  from  three  in-scene  reflectors  with  arbitrary  known  scattering  matrices. 
The  solutions  for  several  sets  of  calibration  targets  with  simple  scattering  matrices  are  first 
presented.  Then,  the  polarimetric  calibration  using  three  targets  with  general  arbitrary 
scattering  matrices  is  derived  using  the  method  of  simultaneous  diagonalization  of  two 
matrices.  A  transformation  matrix  is  found  to  convert  the  general  scattering  matrices  into 
the  simple  cases,  and  the  problem  is  solved  in  the  transformed  domain.  The  solutions  to 
the  original  problem  then  can  be  expressed  in  terms  of  the  solutions  obtained  for  the  simple 
scattering  matrices.  All  possible  combinations  of  calibration  targets  axe  discussed  and  the 


20 


solutions  of  each  cases  are  presented.  Thus,  if  three  scatterers  with  known  scattering 
matrices  are  known  to  exist  within  a  radar  image,  then  the  whole  image  can  be  calibrated 
using  the  exact  solution  presented.  The  effects  of  misalignment  of  calibration  targets  and 
of  receiver  noise  are  also  illustrated  for  several  sets  of  calibration  targets. 

A  multivariate  K-distribution  is  proposed  to  model  the  statistics  of  fully  polari- 
metric  radar  data  from  earth  terrain  with  polaudzations  HH,  HV,  VH,  and  VV.  In  this 
approach,  correlated  polarizations  of  radar  signals,  as  characterized  by  a  covariance  ma¬ 
trix,  sire  treated  as  the  sum  of  N  n-dimensional  random  vectors;  N  obeys  the  negative 
binomial  distribution  with  a  parauneter  a  and  mean  N.  Subsequently,  an  n-dimensionai 
K-distribution,  with  either  zero  or  nonzero  mean,  is  developed  in  the  limit  of  infinite  N 
or  illuminated  area.  The  Probability  Density  Function  (PDF)  of  the  K-distributed  vector 
normalized  by  its  Euclidean  norm  is  independent  of  the  parameter  a  and  is  the  same  as  that 
derived  from  a  zero-mean  Gaussian-distributed  random  vector.  The  above  model  is  well 
supported  by  experimental  data  provided  by  MIT  Lincoln  Laboratory  and  the  Jet  Propul¬ 
sion  Laboratory  in  the  form  of  polarimetric  measurements.  The  results  are  illustrated  by 
comparing  the  higher-order  normalized  intensity  moments  and  Cumulative  Density  Func¬ 
tions  (CDF)  of  the  experimental  data  with  theoretical  results  of  the  K-distribution. 

The  three-layer  random  medium  model  is  developed  for  microwave  remote  sensing 
of  snow-covered  sea  ice.  The  electromagnetic  wave  theory  and  strong  fluctuation  theory 
are  employed  to  study  the  propagation  and  volume  scattering  of  electromagnetic  waves  in 
the  medium.  With  the  application  of  the  Feynman  diagrammatic  technique  and  the  renor¬ 
malization  method,  mean  fields  for  the  isotropic  and  anisotropic  random  media  are  derived 
under  the  bilocal  approximation.  Then,  the  effective  permittivities  for  both  random  media 
are  obtained  from  the  dispersion  relations  of  the  mean  fields.  Further,  with  the  discrete- 
scatterer  concept  for  two-phase  mixtures,  the  scattering  parts  of  effective  permittivities 
are  computed,  in  the  low-frequency  limit,  for  both  isotropic  and  anisotropic  random  media 
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with  specified  correlation  functions.  The  distorted  Born  approximation  is  then  used  to 
compute  the  co-polarized  and  cross-polarized  backscattering  coefficients  which  are  com¬ 
pared  with  scatterometer  data  at  9  and  13  GHz  for  bare  and  dry-snow  covered  thick  first 
yesLT  sea  ice  taken  at  Point  Barrow,  NWT. 

For  active  and  passive  microwave  remote  sensing,  the  correlation  function  used  in 
the  random  medium  model  provides  a  direct  link  between  electrical  behaviors  and  physical 
properties  of  geophysical  media.  The  distribution,  shape,  size,  and  orientation  of  embedded 
inhomogeneities,  such  as  ice  grains  in  snow  and  brine  inclusions  in  sea  ice,  can  be  charac¬ 
terized  by  the  functional  form  of  the  correlation  function,  the  variance,  and  the  correlation 
lengths.  Based  on  the  probability  theory,  anadytical  expressions  of  correlation  functions 
for  two-phase  mixtures  with  randomly  distributed  inclusions  of  spherical  and  spheroidal 
shapes  are  derived.  It  is  shown  that  the  functional  form  of  the  correlation  function  is 
determined  by  the  shape  and  orientation  of  inclusions  while  correlation  lengths  are  related 
to  the  fractional  volume  of  the  scatterers  and  the  total  common  surface  area. 

Supervised  and  unsupervised  classification  procedures  are  developed  and  applied  to 
Synthetic  Aperture  Radar  (SAR)  polarimetric  images  in  order  to  identify  their  various 
earth  terrain  components.  For  supervised  classification  processing,  the  Bayes  technique 
is  used  to  classify  fully  polarimetric  and  normalized  polarimetric  SAR  data.  Simpler  po¬ 
larimetric  discriminates,  such  as  the  absolute  and  norm^dized  magnitude  response  of  the 
individual  receiver  channel  returns,  in  addition  to  the  phase  difference  between  the  receiver 
channels,  are  also  considered.  Another  processing  algorithm,  based  on  comparing  general 
properties  of  the  Stokes  parameters  of  the  scattered  wave  to  that  of  simple  scattering 
models,  is  also  discussed.  This  algorithm,  which  is  an  unsupervised  technique,  classifies 
terrain  elements  based  on  the  relationship  between  the  orientation  angle  and  handedness 
of  the  transmitting  and  receiving  polarization  states.  These  classification  procedures  have 
been  applied  to  San  Francisco  Bay  and  Traverse  City  SAR  images,  supplied  by  the  Jet 
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Propulsion  Laboratory.  It  is  shown  that  supervised  classification  yields  the  best  overall 
performsince  when  accurate  classifier  training  data  is  used,  whereas  unsupervised  classifi¬ 
cation  is  applicable  when  training  data  is  not  avzulable. 

There  is  considerable  interest  in  determining  the  optimal  polarizations  that  maxi¬ 
mize  contrast  between  two  scattering  classes  in  polarimetric  radar  images.  A  systematic 
approach  is  developed  to  obtain  the  optimal  polarimetric  matched  filter,  i.e.,  that  filter 
which  produces  maximum  contrast  between  two  scattering  classes.  The  meiximization  pro¬ 
cedure  involves  solving  an  eigenvalue  problem  where  the  eigenvector  corresponding  to  the 
mzLximum  contrast  ratio  is  optimal  polarimetric  matched  filter.  To  exhibit  the  physical 
significsmce  of  this  filter,  it  is  transformed  into  its  associated  transmitting  and  receiving 
polarization  states,  written  in  terms  of  horizontal  and  vertical  vector  components.  For  the 
special  case  where  the  transmitting  polarization  is  fixed,  the  receiving  polarization  which 
maximizes  the  contrast  ratio  is  also  obtained.  Polarimetric  filtering  is  then  applied  to  syn¬ 
thetic  aperture  radar  images  obtained  from  the  Jet  Propulsion  Laboratory.  It  is  shown, 
both  numerically  and  through  the  use  of  radau'  imagery,  that  maximum  image  contrast  can 
be  realized  when  data  is  processed  with  the  optimal  polarimetric  matched  filter. 

Earth  terraiins  are  modeled  by  a  two-layer  configuration  with  a  random  permittivity 
described  by  a  three-dimensional  correlation  function  horizontad  and  vertical  correlation 
lengths  and  variances.  Using  the  wave  theory  with  Born  approximations  carried  to  the 
second  order,  this  model  is  applied  to  derive  the  polarimetric  backscattering  coefficients 
of  the  Mueller  and  covariance  matrices.  From  a  physical  point  of  view,  the  Born  first- 
and  second-order  approximations  account  for  the  single  and  double  scattering  processes, 
respectively. 


For  the  isotropic  random  permittivity  configuration,  five  out  of  the  nine  elements 
of  the  covariance  matrix  are  zero  under  the  Born  first-order  approximation.  For  the  uni¬ 
axial  tilted  random  permittivity  case,  the  covariance  matrix  does  not  contain  any  zero 
elements.  Furthermore,  in  order  to  account  for  different  directions  of  scattering  especially 
for  vegetation,  an  averaging  scheme  is  developed  in  which  the  backscattering  coefficients 
are  averaged  over  the  azimuthal  direction.  In  this  ca^e,  the  covariance  matrix  is  described 
by  four  zero  elements  though  the  tilt  angle  is  different  from  zero.  Applying  the  Born  ap¬ 
proximation  to  the  second-order,  the  covariance  matrix  is  computed  for  the  isotropic  and 
the  uniaxial  untilted  random  permittivity  configurations.  The  covariance  matrix  has  four 
zero  elements,  and  a  depolarization  factor  is  obtained,  even  for  the  isotropic  case. 

To  describe  the  effect  of  the  random  medium  on  electromagnetic  waves,  the  strong 
permittivity  fluctuation  theory,  which  accounts  for  the  losses  due  to  both  of  the  absorption 
and  the  scattering,  is  used  to  compute  the  effective  permittivity  of  the  medium.  For 
a  mixture  of  two  components,  only  the  frequency,  the  correlation  lengths,  the  fractional 
volume,  and  the  permittivities  of  the  two  constituents  are  needed  to  obtain  the  polarimetric 
backscattering  coefficients.  Theoretical  predictions  are  illustrated  by  comparing  the  results 
with  experimental  data  for  vegetation  fields  and  sea  ice. 

A  mathematically  rigorous  and  fully  polarimetric  radar  clutter  model  is  developed 
to  ev£iluate  the  radar  backscatter  from  various  types  of  terrain  clutter  such  as  forested 
areas,  vegetation  canopies,  snow  covered  terrains,  or  ice  fields.  With  this  model,  we  can 
calculate  the  radar  backscattering  coefficients  for  the  multi-channel  polarimetric  radar 

returns,  in  addition  to  the  complex  cross  correlation  coefficients  between  elements  of  the 
polarimetric  measurement  vector.  The  complete  polarization  covariance  matrix  can  be 
computed  and  the  scattering  properties  of  the  clutter  environment  characterized  over  a 
broad  range  of  incident  angles  and  frequencies. 
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Algorithms  have  been  developed  to  identify  and  classify  different  features  in  po- 
larimetric  SAR  images.  Optimal  schemes  are  designed  to  enhance  radar  target  detec¬ 
tion  and  discrimination.  The  fully  polcirimetric,  multi-frequency,  multi-incident  angle, and 
multi-layer  random  medium  model  are  employed  in  the  evaluation  of  synthetic  clutter 
signatures  that  can  be  used  in  a  hardware-in-the-loop  test  system.  With  the  Generzil- 
ized  K-distribution,  better  description  of  the  statistics  of  the  polarimetric  SAR  data  are 
obtciined.  The  identification  and  classification  procedures  has  been  applied  to  JPL  po¬ 
larimetric  SAR  images  containing  different  features.  RCS  codes  are  improved  to  predict 
radar  cross  sections  from  various  targets.  PracticcJly,  the  developments  reported  in  this 
document  are  highly  applicable  to  the  target-clutter  detection  problems  such  as  land  mines 
disguised  under  vegetation  or  soil  covers. 

Earth  terrain  media,  such  as  vegetation,  forest,  snow,  and  ice  exhibit  strong  vol¬ 
ume  scattering  effects.  To  study  polarimetric  radar  backscatter  of  earth  terrains  from  the 
point  of  view  of  electromagnetic  wave  theory,  we  use  a  layered  random  medium  model  to 
characterize  the  terrain  clutter.  The  random  medium  is  described  by  a  background  per¬ 
mittivity  with  a  fluctuating  component;  the  randomness  of  the  fluctuation  is  characterized 
by  a  three-dimensional  correlation  function  with  a  variance  and  horizontal  and  vertical 
correlation  lengths.  The  variance  corresponds  to  the  strength  of  the  fluctuation,  whereas 
the  correlation  lengths  coincide  with  the  geometrical  size  of  the  basic  scattering  elements. 
The  polarimetric  backscattering  coefficients  can  be  obtained  from  the  electromagnetic  wave 
theory  by  calculating  the  covariance  matrix  of  the  polarimetric  meeisurement  vector. 
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The  problem  of  scattering  by  a  layer  of  random  medium  (e.g.  snow  field,  meadow, 
etc.)  can  be  solved  with  a  number  of  different  techniques.  We  developed  a  wave  approach 
by  applying  the  Born  approximation.  An  integrjd  equation  is  formulated  for  the  electric 
field  by  using  the  unperturbed  Green’s  function  for  a  layered  remdom  medium  in  the 
absence  of  permittivity  fluctuations.  The  random  fluctuations  are  treated  as  induced 
scattering  sources  in  the  integral  equation,  which  will  then  be  solved  by  iteration  to  obtain 
a  Neumann  series.  To  include  the  depolarization  effects,  the  Born  approximation  is  carried 
out  to  the  second  order.  Physically,  the  Born  first-order  ^lnd  second-ordei  approximations 
account,  respectively,  for  a  single  and  double  scattering  process.  The  calculated  covariance 
matrix  for  a  layered  isotropic  random  medium  has  four  of  its  dements  equal  to  zero, 
indicating  absence  of  correlation  between  the  HV  and  HF  terms  and  HV  and  VV  terms. 
The  theoretical  results  are  shown  to  be  consis'  nt  with  measurement  data  obtained  from 
the  MIT  Lincoln  Laboratory. 

With  the  same  approach,  the  po-  j  .metric  scattering  effects  from  a  layer  of 
anisotropic  random  medium  (e.g.  sea  ice)  or  from  a  layer  of  isotropic  random  medium 
on  top  of  another  layer  of  anisotropic  random  medium  (e.g.  tree  canopy,  snow-covered 
sea  ice,  etc.)  are  investigated  and  the  complete  covariance  matrices  with  nine  non-zero 
elements  are  computed.  Theoretically  c€dculated  covariance  matrices  wiU  be  compared 
with  experimentzil  polarimetric  HH,  HV,  and  VV  clutter  data  collected  from  the  various 
terrains.  Physical  interpretations  will  be  illustrated  along  with  the  other  properties  of  the 
covariance  matrix  elements. 
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The  evaluation  of  clutter  backseat  ter  coefficients  for  HH,  HV  and  VV  polarimetric 
returns  is  useful  in  the  design  and  analysis  of  optimal  radar  tsu-get  detection,  discrim¬ 
ination  and  classification  schemes.  Another  important  application  is  the  generation  of 
random  clutter  returns  for  Monte-Carlo  simulations.  The  application  of  the  fully  polari¬ 
metric  clutter  model  in  the  terrain  classification  will  be  illustrated.  The  optimal  (Bayes) 
fully  polarimetric  and  normalized  polarimetric  classifiers  as  well  as  simpler  polarimetric 
discriminants  will  also  be  presented.  In  addition,  the  optimal  polarimetric  matched  filter 
will  be  introduced  as  a  means  of  enhancing  the  contrast  between  two  areas  within  radar 
polarimetry.  Fully  polarimetric  synthetic  aperture  radar  images,  obtained  from  the  Jet 
Propulsion  Laboratory,  will  be  used  to  demonstrate  various  terrain  classification  schemes 
and  the  contrast  enhancement  obtained  using  the  polarization  matched  filter. 

A  clutter  model  is  used  to  simulate  fully  polarimetric  returns  for  a  coherent,  stepped- 
frequency  radar.  The  objective  is  to  create  site- dependent  synthetic  clutter  signatures 
that  can  be  utilized  in  a  hardware-in-the-loop  test  system.  The  fully  polarimetric,  multi¬ 
frequency,  and  multi-incident  angle  two-layer  anisotropic  random  medium  model  is  em¬ 
ployed  to  compute  the  normalized  backscatter  coefficients  of  terrain  clutter.  Polarization 
covariance  matrices  are  calculated  for  each  of  N  high  resolution  range  bins,  at  each  of  the 
M  discrete  frequencies  that  comprise  the  stepped-frequency  bandwidth.  The  covariance 
matrices  are  decomposed,  multiplied  by  complex  Gaussian  noise,  and  weighted  according 
to  the  random  product  model  to  account  for  the  spatial  variability  of  clutter.  This  gen¬ 
erates  the  normalized  electric  field  backscattered  from  each  of  the  N  range  bins,  at  each 
of  the  M  discrete  frequencies.  These  fields  are  coherently  added,  taking  into  account  the 
effects  of  terrain  shadowing  and  overlay,  in  order  to  realize  the  single-frequency  polarimet¬ 
ric  return  that  a  radar  would  measure  from  the  specified  terrain.  The  radar  return  for 
each  of  the  other  discrete  frequencies  is  computed  in  a  similar  manner.  In  addition,  the 
real-time  implementation  of  this  algorithm  is  considered  in  the  context  of  synthetic  wide¬ 
band  processing.  The  result  of  this  procedure  yields  the  clutter’s  coherent  phase-history 
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profile.  The  clutter  phase-history  returns  can  then  be  coherently  superimposed  on  the 
target  phase-history  returns.  The  combined  (or  clutter  only)  returns  can  be  reduced  to 
obtain  either  (1)  the  coherent  high  resolution  range  profile  (HRRP)  or  (2)  the  noncoherent 
autocorrelation  range  profile. 

Existing  radar  cross  section  (RCS)  prediction  algorithms  using  the  Physical  Optics 
and  the  Physical  Theory  of  Diffraction  techniques  have  been  shown  to  give  erroneous 
results  for  cases  in  which  multiple  scattering  between  portions  of  the  target  is  significant. 
In  particular,  double  reflections  between  flat  plate  surfaces  are  important  for  geometries 
resembling  the  corner  reflector,  and  this  effect  is  neglected  by  current  algorithms  which 
include  only  single  surface  interactions.  To  overcome  this  limitation  a  method  is  presented 
of  extending  the  conventional  results  obtained  using  Physical  Optics  to  include  a  field  due 
to  double  reflections  between  two  polygonal  plate  surfaces. 

An  integrcil  expression  is  derived,  applying  the  tangent  plane  approximation  at  each 
surface,  and  correctly  accounting  for  the  potentially  near-field  interactions  between  the 
two  perfectly  conducting  plates.  The  initial  complexity  of  this  expression  necessitates 
numerical  evaluation  of  the  integrals  and  results  in  a  large  computational  burden.  The 
requirement  of  numerical  integration  is  eliminated,  however,  by  applying  the  method  of 
Stationary  Phase  and  expanding  quadratically  the  exact  inter-plate  phase  term.  The  two 
surface  integrations  are  replaced  by  line  integrals  using  Stokes’  Theorem,  and  for  polygonal 
plates,  the  contour  is  subdivided  into  linear  segments  on  which  the  integration  is  performed 
analytically.  In  the  case  where  the  conducting  plates  are  loaded  with  one  or  more  layers 
of  dielectric,  the  equivalent  electric  and  magnetic  surface  currents  are  derived  using  the 
reflection  coefRcients  at  each  plate. 
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The  resulting  field  expressions  are  derived  with  sufficient  generality  so  as  to  permit 
incorporation  in  standard  RCS  prediction  codes  allowing  arbitrary  geometry  and  illumi¬ 
nation.  Comparisons  with  measurements  of  several  simple  two-plate  targets  confirm  the 
accuracy  of  the  predicted  field  and  demonstrate  a  significant  improvement  over  existing  al¬ 
gorithms  which  neglect  multiple  reflections  or  which  etssume  a  fair  field  interaction  between 
the  surfaces. 

Polaximetric  radar  backscatter  data  have  been  used  extensively  to  classify  terrain 
cover.  Since  it  is  difficult  to  calibrate  out  the  effects  of  amplitude  and  phase  errors  induced 
by  atmospheric  effects,  path  loss,  etc.,  absolute  amplitude  and  phase  of  radar  returns  are 
not  reliable  features  for  terrain  classification  purposes.  The  use  of  normalized  polarimetric 
data  is  proposed  such  that  only  the  relative  magnitudes  and  phases  will  be  utilized  to 
discriminate  different  terrain  elements.  It  is  shown  that  the  Bayes  clzissification  error 
does  not  depend  on  the  form  of  the  normalization  function  if  the  unknown  radar  system 
calibration  factor  is  modeled  as  a  multiplicative  term  in  the  received  signal.  Assuming  a 
multivariate  Gaussian  distribution  for  the  un-normalized  polarimetric  data,  the  Probability 
Density  Function  (PDF)  of  the  normalized  data  and  the  corresponding  Bayes  classifier 
distance  measure  for  the  normalized  data  are  derived.  Furthermore,  by  assuming  a  specific 
form  of  the  covariance  matrix  for  the  polarimetric  data,  exact  PDFs  are  given  for  HH,  HV, 
VV  and  span  type  normalization  schemes.  Corresponding  cleissification  errors  axe  evaluated 
to  verify  their  independence  from  all  normalization  functions. 

The  scattering  of  electromagnetic  waves  from  a  randomly  perturbed  periodic  sur¬ 
face  is  formulated  by  the  Extended  Boundary  Condition  (EBC)  method  and  solved  by  the 
Small  Perturbation  Method  (SPM).  The  scattering  from  periodic  surface  is  solved  exactly 
and  this  solution  is  used  in  the  SPM  to  solve  for  the  surface  currents  and  scattered  fields 
up  to  the  second  order.  The  random  perturbation  is  modeled  as  a  Gaussian  random  pro¬ 
cess.  The  theoretical  results  are  illustrated  by  calculating  the  bistatic  and  backscattering 
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coefHcients.  It  is  shown  that  as  the  correlation  length  of  the  rauidom  roughness  increases, 
the  bistatic  scattering  pattern  of  the  scattered  fields  show  several  beams  associated  with 
each  Bragg  diffraction  direction  of  the  periodic  surface.  When  the  correlation  length  be¬ 
comes  smaller,  then  the  shape  of  the  beams  become  broader.  The  results  obtained  using 
the  EBC/SPM  method  is  also  compared  with  the  results  obtained  using  the  Kirchhoff 
approximation.  It  is  shown  that  the  Kirchhoff  approximation  results  show  quite  a  good 
agreement  with  EBC/SPM  method  results  for  the  hh  and  vv  polarized  backscattering  co¬ 
efficients  for  small  angles  of  incidence.  However,  the  Kirchhoff  approximation  does  not 
give  depolarized  returns  in  the  backscattering  direction  whereas  the  results  obtained  using 
the  EBC/SPM  method  give  significant  depolarized  returns  when  the  incident  direction  is 
not  perpendicular  to  the  row  direction  of  the  periodic  surface. 

A  general  mixing  formula  is  derived  for  discrete  scatterers  immersed  in  a  host 
medium.  The  inclusion  particles  are  assumed  to  be  ellipsoidal.  The  electric  field  inside 
the  scatterers  is  determined  by  quasistatic  analysis,  assuming  the  diameter  of  the  inclusion 
particles  to  be  much  smaller  than  the  wavelength.  The  results  are  applicable  to  general 
multiphase  mixtures,  and  the  scattering  ellipsoids  of  the  different  phases  can  have  different 
sizes  and  arbitrary  ellipticity  distribution  and  axis  orientation,  i.e.,  the  mixture  may  be 
isotropic  or  anisotropic.  The  resulting  mixing  formula  is  nonlinear  and  implicit  for  the  ef¬ 
fective  complex  dielectric  constant,  because  the  approach  in  calculating  the  internal  field  of 
scatterers  is  self-consistent.  Still,  the  form  is  especially  smtable  for  iterative  solution.  The 
formula  contains  a  quantity  called  the  apparent  permittivity,  and  with  different  choices 
of  this  quantity,  the  result  leads  to  the  generalized  Lorentz  -  Lorenz  formula,  the  general¬ 
ized  Polder  -  van  Santen  formula,  and  the  generalized  coherent  potential  -  quasicrystalline 
approximation  formula. 
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We  have  used  the  strong  fluctuation  theory  to  derive  the  backscattering  cross  sec¬ 
tions.  The  study  of  the  strong  fluctuation  theory  for  a  bounded  layer  of  random  discrete 
scatterers  is  further  extended  to  include  higher  order  co-polarized  2Lnd  cross-polarized  sec¬ 
ond  moments.  The  backscattering  cross  sections  per  unit  area  are  calculated  by  including 
the  mutual  coherence  of  the  fields  due  to  the  coincidental  ray  paths  and  that  due  to  the 
opposite  ray  paths  which  are  corresponding  to  the  ladder  and  cross  terms  in  the  Feynman 
diagrammatic  representation.  It  is  proved  that  the  contributions  from  ladder  and  cross 
terms  for  co-polarized  backscattering  cross  sections  are  the  same,  while  the  contributions 
for  the  cross-polarized  ones  are  of  the  same  order.  The  bistatic  scattering  coefficients  in 
the  second-order  approximation  for  both  the  ladder  and  cross  terms  are  also  obtzuned.  The 
enhancement  in  the  backscattering  direction  can  be  attributed  to  the  contributions  from 
the  cross  terms. 

A  two-layer  anisotropic  random  medium  model  is  developed  for  the  active  and  pas¬ 
sive  microwave  remote  sensing  of  ice  fields.  The  dyadic  Green’s  function  for  this  two-layer 
anisotropic  medium  is  derived.  With  a  specified  correlation  function  for  the  randomness 
of  the  dielectric  constant,  the  backscattering  cross  sections  are  crdculated  with  the  Born 
approximation.  It  is  shown  that  the  depolarization  effects  exist  in  the  single-scattering 
process  and  that  the  backscattering  cross  section  of  horizontal  polarization  can  be  greater 
than  that  of  vertical  polarization  even  in  the  half-space  case.  The  principle  of  reciprocity 
and  the  principle  of  energy  conservation  are  invoked  to  calculate  the  brightness  temper¬ 
atures.  The  bistatic  scattering  coefficients  are  first  calculated  and  then  integrated  over 
the  upper  hemisphere  to  be  subtracted  from  unity,  in  order  to  obtain  the  emissivity  for 
the  random  medium  layer.  It  is  shown  that  both  the  absorptive  and  randomly  fluctuat¬ 
ing  properties  of  the  anisotropic  medium  affect  the  behavior  of  the  resulting  brightness 
temperatures  both  in  theory  and  in  actual  controlled  field  measurements.  The  active  and 
passive  results  match  favorably  well  with  the  experimental  data  obtained  from  the  corn 
stalks  with  detailed  ground-truth  information. 
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The  vegetation  canopy  and  snow-covered  ice  field  have  been  studied  with  a  three- 
layer  model,  an  isotropic  random  medium  layer  overlying  an  anisotropic  random  medium. 
We  have  calculated  the  dyadic  Green’s  functions  of  the  three-layer  medium  and  the  scat¬ 
tered  electromagnetic  intensities  with  Born  approximation.  The  backscattering  cross  sec¬ 
tions  are  evaluated  for  active  microwave  remote  sensing.  The  theoretical  approach  can  be 
extended  to  derive  the  bistatic  scattering  coefficients. 

Backscattering  and  propagation  of  electromagnetic  pulses  with  different  shapes  in 
earth  terreiin  media  such  as  snow-ice,  vegetation,  and  forest  are  studied  with  the  random 
medium  model.  The  backscattering  cross  sections  are  computed  with  both  the  Born  and 
the  distorted  Born  approximations.  We  study  the  pulse  shape  distortion  by  considering  a 
two-layer  random  medium  model  with  different  correlation  lengths  and  varying  variances. 
The  degree  of  distortion  is  inversely  proportional  to  the  pulse  width  and  determined  by 

strength.  Both  the  like-polarized 
examined  and  illustrated. 


the  magnitude  of  the  randomly  fluctuating  permittivity 
and  cross-polarized  radar  backscattering  coefficients  are 


32 


PUBLICATIONS  SUPPORT  BY  WATERWAYS  EXPERIMENTAL  STATION 

Applications  of  EMSARS  background  model  to  MMW  vegetation  scattering 
simulation  (C.  C.  Hsu,  Y.  E.  Yang,  R,  T.  Shin,  J.  A.  Kong,  C.  Kohler,  T. 
Nguyen,  H.  Nguyen,  and  J.  Ho),  3rd  Ground  Target  Modeling  and  Validation 
Conference,  Michigan  Technical  University,  Houghton,  MI,  August  18-20,  1992. 

External  calibration  of  polarimetric  radars  using  point  and  distributed  targets 
(S.  H.  Yueh,  J.  A.  Kong,  and  R.  T.  Shin),  Journal  of  Electromagnetic  Waves  and 
Applications^  Vol.  6,  No.  7,  1992. 

Simulation  of  forest  backscatter  as  a  function  of  forest  and  ground  parame¬ 
ters  (A.  Beaudoin,  T.  Le  Toan,  C.  C.  Hsu,  H.  C.  Han,  J.  A.  Kong,  and  R. 
T.  Shin),  International  Geoscience  and  Remote  Sensing  Symposium,  Houston, 
Texas,  USA,  May  26-29,  1992. 

Branching  model  for  vegetation  (S.  H.  Yueh,  J.  A.  Kong,  J.  K.  Jao,  R.  T.  Shin, 
and  T.  Le  Toan),  IEEE  Transactions  on  Geoscience  and  Remote  Sensing,  Vol.  30,  No. 
2,  390-402,  March  1992. 

Theoretical  modeling  for  passive  microwave  remote  sensing  of  earth  Terrain  (J. 
A.  Kong),  URSI  Specisdists  Meeting  on  Microwave  Radiometry  and  Remote 
Sensing,  Boulder,  Colorado,  January  24-26,  1992. 

Theoretical  models  for  polarimetric  microwave  remote  sensing  of  earth  terreun 
(M.  Borgeaud,  J.  A.  Kong,  R.  T.  Shin,  and  S.  V.  Nghiem),  in  Direct  and 
Inverse  Methods  in  Radar  Polarimetry,  Part  2,  1139-1190,  Kluwer  Academic 
Publishers,  The  Netherlands,  1992. 

K-distribution  and  multi-frequency  polarimetric  terrain  radar  clutter  (H.  A. 
Yueh,  J.  A.  Kong,  J.  K.  Jao,  R.  T.  Shin,  H.  A.  Zebker,  and  T.  Le  Toan),  Journal 
of  Electromagnetic  Waves  and  Applications,  Vol.  5,  No.  1,  pp.  1—15,  1991. 

Polarimetric  passive  remote  sensing  of  periodic  surfaces  (M.  E.  Veysoglu,  H.  A. 
Yueh,  R.  T.  Shin  and  J.  A.  Kong),  Journal  of  Electromagnetic  Waves  and  Applications, 
Vol.  5,  No.  3,  267-280,  1991. 

External  cjJibration  of  polarimetric  radars  using  point  and  distributed  targets 
(S.  H.  Yueh,  J.  A.  Kong,  and  R.  T.  Shin),  Proceedings  of  the  Third  Airborne  Synthetic 
Aperture  Radar  (AIRSAR)  Workshop,  JPL,  Pasadena,  CA,  May  23-24,  1991 

Application  of  neural  networks  to  sea  ice  classification  using  polarimetric  sar  im¬ 
ages  (R.  Kwok,  Y.  Hara,  R.  G.  Atkins,  S.  H.  Yueh,  R.  T.  Shin  and  J.  A.  Kong), 
International  Geoscience  and  Remote  Sensing  Symposium  (IGARSS'91), 
Helsinki  University  of  Technology,  Espoo,  Finland,  June  3-6,  1991. 

Electromagnetic  scattering  of  waves  by  random  rough  surface:  a  finite-difference 
time-domain  approach  (C.  H.  Chan,  S.  H.  Lou,  L.  Tsang,  and  J.  A.  Kong), 
Microwave  and  Optical  Technology  Letter*. Vol. 4,  No. 9,  355  —  359,  August  1991. 


33 


Application  of  random  medium  model  to  remote  sensing  of  vegetation  (T.  Le 
Toan,  S.  V.  Nghiem,  J.  A.  Kong,  H.  C.  Han,)  International  Geoscience  and  Remote 
Sensing  Symposium  (IGARSS’91),  Helsinki  University  of  Technology,  Espoo,  Fin¬ 
land,  June  3-6,  1991. 

Multiple  ellipsoidal  species  in  layer  random  medium  model  (S.  V.  Nghiem,  J. 
A.  Kong,  R.  T.  Shin,  T.  Le  Toan),  Progress  in  Electromagnetics  Research  Symposium, 
Cambridge,  Massachusetts,  July  1-5,  1991. 

Polarimetric  passive  remote  sensing  of  an  azimuthally  asymmetric  periodic  soil 
surface  (S.  V.  Nghiem,  M.  E.  Veysoglu,  J.  A.  Kong,  R.  T.  Shin,  K.  O’Neill, 
A.  W.  Lohanick),  Progress  in  Electromagnetics  Research  Symposium,  Cambridge,  Mas¬ 
sachusetts,  July  1-5,  1991. 

Radiative  transfer  theory  for  passive  microwave  remote  sensing  of  a  vegetation 
canopy  (H.  C.  Han,  W.  Au,  S.  H.  Yueh,  R.  T.  Shin,  J.  A.  Kong),  Progress  in 
Electromagnetics  Research  Symposium,  Cambridge,  Massachusetts,  July  1-5,  1991. 

Analytical  solution  of  the  vector  radiative  transfer  equation  with  rough  surface 
boundary  condition  (H.  C.  Han,  J.  A.  Kong,  S.  V.  Nghiem,  T.  Le  Toan),  Progress 
in  Electromagnetics  Research  Symposium,  Cambridge,  Massachusetts  ,  July  1-5,  1991. 

Branching  model  for  polarimetric  remote  sensing  of  vegetation  (S.  H.  Yueh,  J. 
A.  Kong,  J.  K.  Jao,  S.  Ayasli,  R.  T.  Shin),  Progress  in  Electromagnetics  Research 
Symposium,  Cambridge,  Massachusetts  ,  July  1-5,  1991. 

Radiative  transfer  modeling  of  vegetation  clusters  (C.  C.  Hsu,  S.  H.  Yueh,  H. 
C.  Han,  R.  T.  Shin,  J.  A.  Kong),  Progress  in  Electromagnetics  Research  Symposium, 
Cambridge,  Massachusetts,  July  1-5,  1991. 

Polarimetric  passive  remote  sensing  of  periodic  structures  (M.  E.  Veysoglu,  S. 
H.  Yueh,  R.  T.  Shin,  J.  A.  Kong),  Progress  in  Electromagnetics  Research  Symposium, 
Cambridge,  Massachusetts,  July  1-5,  1991. 

A  radar  scattering  model  for  structured  vegetation  (S.  H.  Yueh,  J.  A.  Kong, 
J.  K.  Jao,  R.  T.  Shin,  T.  Le  Toan),  Progress  in  Electromagnetics  Research  Symposium, 
Cambridge,  Massachusetts,  July  1-5,  1991. 

Application  of  neural  networks  to  polarimetric  sar  image  classification  (Y.  Hara, 
R.  G.  Atkins,  S.  H.  Yueh,  R.  T.  Shin,  J.  A.  Kong,  R.  Kwok),  Progress  in  Electro¬ 
magnetics  Research  Symposium,  Cambridge,  Massachusetts,  July  1-5,  1991. 

Covariance  of  phase  and  amplitude  fluctuations  of  electromagnetic  waves  prop¬ 
agating  through  a  random  medium  (J.  G.  Fleischman,  S.  Ayasli,  R.  T.  Shin,  N. 
C.  Chu,  S.  H.  Yueh,  S.  V.  Nghiem),  Progress  in  Electromagnetics  Research  Symposium, 
Cambridge,  Massachusetts,  July  1-5,  1991. 

Observation  and  modelling  of  polarimetric  sar  signatures  of  forests  (T.  Le  Toan, 
A.  Beaudoin,  S.  Goze,  S.  V.  Nghiem,  J.  A.  Kong),  Progress  in  Electromagnetics 
Research  Symposium,  Cambridge,  Massachusetts,  July  1-5,  1991. 


34 


Electromagnetic  wave  modeling  for  remote  sensing  (S.  V.  Nghiem,  J.  A.  Kong, 
and  T.  Le  Toan),  International  Conference  on  Directions  in  Electromagnetic  Wave  Model¬ 
ing,  New  York,  October  22-24, 1990. 

Simulation  of  electromagnetic  phenomema  using  a  finite  difference-time  domain 
technique  (K.  Li,  M.  A.  Tassoudji,  R.  T.  Shin,  and  J.  A.  Kong),  7th  Annual  Re¬ 
view  of  Progress  in  Applied  Computational  Electromagnetics,  Naval  Postgraduate  School, 
Monterey,  pp.  38-55,  March  18-22,  1991. 

Layer  model  with  random  spheroidal  scatterers  for  remote  sensing  of  vegetation 
canopy  (S.  V.  Nghiem,  J.  A.  Kong,  H.  C.  Han,  T.  Le  Toan,  and  M.  Borgeaud), 
to  be  submitted  to  Journal  of  Electromagnetic  Waves  and  Applications. 

Finite  difference  method  for  electromagnetic  scattering  problems  (C.  F.  Lee, 
R.  T.  Shin,  and  J.  A.  Kong),  Progress  In  Electromagnetics  Research,  edited  by  J.  A. 
Kong,  Vol.  4,  Chapter  11,  p.  373-442,  Elsevier,  New  York,  1990. 

Calibration  of  polarimetric  radars  using  in-scene  reflectors,  (S.  H.  Yueh,  J.  A. 
Kong,  and  R.  T.  Shin),  Progress  In  Electromagnetic  Research,  edited  by  J.  A.  Kong, 
Chapter  9,  Vol.  3,  p.  451-510,  Elsevier,  New  York,  1990. 

Classification  and  majcimum  contrast  of  earth  terrain  using  polarimetric  syn¬ 
thetic  aperture  radar  images,  (J.  A.  Kong,  S.  H.  Yueh,  H.  H.  Lim,  R.  T.  Shin, 
and  J.  J.  van  Zyl),  Progress  In  Electromagnetic  Research,  edited  by  J.  A.  Kong,  Chap¬ 
ter  6,  Vol.  3,  p.  327-370,  Elsevier,  New  York,  1990. 

K-distribution  and  polarimetric  terrain  radar  clutter,  (S.  H.  Yueh,  J.  A.  Kong, 
J.  K.  Jao,  R.  T.  Shin,  H.  A.  Zebker,  T.  Le  Toan,  and  H.  Ottl),  Progress  In 
Electromagnetic  Research,  edited  by  J.  A.  Kong,  Chapter  4,  Vol.  3,  p.  237-275, 
Elsevier,  New  York,  1990. 

Polarimetric  remote  sensing  of  geophysical  media  with  layer  random  medium 
model  (S.  V.  Nghiem,  M.  Borgeaud,  J.  A.  Kong,  and  R.  T.  Shin),  Progress 
In  Electromagnetic  Research,  edited  by  J.  A.  Kong,  Chapter  1,  Vol.  3,  p.  1—73, 
Elsevier,  New  York,  1990. 

Classification  of  earth  terrain  using  synthetic  aperture  radar  images  (H.  Lim, 
A.  A.  Swartz,  H.  A.  Yueh,  J.  A.  Kong,  R.  T.  Shin,  and  J.  J.  Van  Zyl),  Journal 
of  Geophysical  Research,  Vol.  94,  No.  B6,  7049—7057,  June  10,  1989. 

K-distribution  and  polarimetric  terrain  radar  clutter  (H.  A.  Yueh,  J.  A.  Kong, 
J.  K.  Jao,  R.  T.  Shin,  and  L.  M.  Novak),  Journal  of  Electromagnetic  Waves  and 
Applications,  Vol.  3,  No.  8,  747—768,  1989. 

Calibration  of  polarimetric  radar  using  in-scene  reflectors  (H.  A.  Yueh,  J.  A. 
Kong,  R.  M.  Barnes,  and  R.  T.  Shin),  Journal  of  Electromagnetic  Waves  and  Applica¬ 
tions,  Vol.  4,  No.  1,  27-49,  1990. 


35 


A  neural  net  method  for  high  range  resolution  target  classification  (R.  G. 
Atkins,  R.  T.  Shin,  and  J.  A.  Kong),  Progress  In  Electromagnetics  Research,  edited 
by  J.  A.  Kong,  Chapter  7,  Vol.  4,  p.  255-292,  Elsevier,  New  York,  1990. 

Scattering  from  randomly  oriented  scatterers  with  strong  permittivity  fluctua¬ 
tions,  (H.  A.  Yueh,  R.  T.  Shin,  and  J.  A.  Kong),  Journal  of  Electromagnetic  Waves 
and  Applications,  Vol.  4,  No.  10,  p.  983—1004,  1990. 

Polarimetric  remote  sensing  of  earth  terrain  with  two-layer  random  medium 
model,  (M.  Borgeaud,  J.A.  Kong,  R.T.  Shin,  and  S.  V.  Nghiem),  Progress  in 
Electromagnetics  Research  Symposium,  Boston,  July  25—26,  1989. 

Three-layer  random  medium  model  for  fully  polarimetric  remote  sensing  of 
geophysical  media  (S.  V.  Nghiem,  F.  C.  Lin,  J.  A.  Kong,  R.  T.  Shin,  and  H.  A. 
Yueh),  Progress  in  Electromagnetics  Research  Symposium,  Boston,  July  25—26,  1989. 

Faraday  polarization  fluctuations  in  transionspheric  polarimetric  VLF  waves  (S. 
V.  Nghiem  and  J.  A.  Kong),  Progress  in  Electromagnetics  Research  Symposium,  Boston, 
July  25-26,  1989. 

Calibration  of  polarimetric  radars  using  in-scene  reflectors  (H.  A.  Yueh,  J.  A. 
Kong,  R.  M.  Barnes  and  R.  T.  Shin),  Progress  in  Electromagnetics  Research  Symposium, 
Boston,  July  25-26,  1989. 

K-distribution  and  polarimetric  terrain  radar  clutter  (H.  A.  Yueh,  J.  A.  Kong, 
J.  K.  Jao,  R.  T.  Shin,  and  L.  M.  Novak),  Progress  in  Electromagnetics  Research 
Symposium,  Boston,  July  25—26,  1989. 

Theoretical  models  for  polarimetric  microwave  remote  sensing  of  earth  terrain 
(M.  Borgeaud,  S.  V.  Nghiem,  R.  T.  Shin,  and  J.  A.  Kong),  Journal  of  Electromag¬ 
netic  Waves  and  Application,  Vol.  3,  No.  1,  61—81,  1989. 

The  optimal  polarizations  for  achieving  maximum  contrast  in  radar  polarimetry 
(A.  A.  Swartz,  H.  A.  Yueh,  J.  A.  Kong,  L.  M.  Novak,  and  R.  T.  Shin),  Journal 
of  Geophysical  Research,  Vol.  93,  No.  B12,  15252-15260,  Dec.  1988. 

Theoretical  models  for  polarimetric  microwave  remote  sensing  of  earth  terrain 
(M.  Borgeaud,  J.  A.  Kong,  R.  T.  Shin,  and  S.  V.  Nghiem)  Proceedings  of  the  1988 
NATO  Advanced  Research  Workshop,  Nuremberg,  F.  R.  Germany,  September,  1988. 

Polarimetric  clutter  modeling:  theorv  and  application  (J.  A.  Kong,  F.  C.  Lin, 
M.  Borgeaud,  H.  A.  Yueh,  A.  A.  Swartz,  H.  H.  Lim,  L.  M.  Novak,  and  R.  T. 
Shin),  GACIAC PR-88-0S,  The  Polarimetric  Technology  Workshop,  Rocket  Auditorium, 
Redstone  Arsenal,  U.  S.  Army  Missile  Command,  Huntsville,  Alabama,  August 
16-19,  1988. 

Polarimetric  remote  sensing  of  earth  terrain  with  three-layer  random  medium 
model  (S.  V.  Nghiem,  F.  C.  Lin,  J.  A.  Kong,  R.  T.  Shin,  and  H.  A.  Yueh),  GA¬ 
CIAC  PR-88-0S,  The  Polarimetric  Technology  Workshop,  Rocket  Auditorium,  Redstone 
Arsenal,  U.  S.  Army  Missile  Command,  Huntsville,  AL,  Aug.  16-19,  1988. 


36 


Radax  range  profile  simulation  of  terrain  clutter  using  the  random  medium 
model  (A.  A.  Swartz,  L.  M.  Novak,  R.  T.  Shin,  D.  A.  McPherson,  A.  V.  Saylor, 
F.  C.  Lin,  H.  A.  Yueh,  J.  A.  Kong)  GACIAC  PR-88-0S,  The  Polarimetric  Technology 
Workshop,  Rocket  Auditorium,  Redstone  Arsenal,  U.  S.  Army  Missile  Command, 
Huntsville,  AL,  Aug.  16-19,  1988. 

Three-layer  random  medium  model  for  polarimetric  remote  sensing  of  earth 
terrain  (H.  A.  Yueh,  S.  V.  Nghiem,  F.  C.  Lin,  J.  A.  Kong,  and  R.  T.  Shin), 
IEEE  APS  International  Symposium  and  URSI  Radio  Science  Meeting,  Syracuse,  June  6 

-  10,  1988. 

The  optimal  polarizations  for  achieving  maximum  contrast  in  rad2Lr  polarimetry 
(A.  A.  Swartz,  L.  M.  Novak,  R.  T.  Shin,  H.  A.  Yueh,  and  J.  A.  Kong),  IEEE 
APS  International  Symposium  and  URSI  Radio  Science  Meeting,  Syracuse,  June  6  —  10, 
1988. 

K-distribution  and  polarimetric  terrain  radar  clutter  (H.  A.  Yueh,  J.  A.  Kong, 
J.  K.  Jao,  and  R.  T.  Shin),  IEEE  APS  International  Symposium  and  URSI  Radio  Science 
Meeting,  Syracuse,  June  6  -  10,  1988. 

Classification  of  earth  terrain  using  polarimetric  synthetic  aperture  radar  im¬ 
agery  (H.  H.  Lim,  A.  A.  Swartz,  H.  A.  Yueh,  J.  A.  Kong  and  J.  J.  van  Zyl), 
IEEE  APS  International  Symposium  and  URSI  Radio  Science  Meeting,  Syracuse,  June  6 

-  10,  1988. 

Theoretical  model  for  snow-covered  sea  ice  (F.  C.  Lin,  J.  A.  Kong,  R.  T.  Shin, 
A.  J.  Gow,  and  D.  Perovich),  IEEE  APS  International  Symposium  and  URSI  Radio 
Science  Meeting,  Syracuse,  June  6  -  10,  1988. 

Effect  of  multiple  scattering  on  the  radar  cross  section  of  polygonal  plate  struc¬ 
tures  (R.  G.  Atkins  and  R.  T.  Shin),  IEEE  APS  International  Symposium  and  URSI 
Radio  Science  Meeting,  Syracuse,  June  6  —  10,  1988. 

Bayes  classification  of  terrain  cover  using  norm^llized  polarimetric  data  (H.  A. 
Yueh,  A.  A.  Swartz,  J.  A.  Kong,  R.  T.  Shin,  and  L.  M.  Novak),  Journal  of 
Geophysical  Research,  Vol.  93,  No.  B12,  15261-15267,  December  1988. 

Scattering  of  electromagnetic  waves  from  a  periodic  surface  with  random  rough¬ 
ness  (H.  A.  Yueh,  R.  T.  Shin,  and  J.  A.  Kong),  Journal  of  Applied  Physics,  vol. 64, 
no.4,  1657-1670,  August  1988. 

Scattering  of  electromagnetic  waves  from  a  periodic  surface  with  random  rough¬ 
ness  (H.  A.  Yueh,  R.  T.  Shin,  and  J.  A.  Kong),  SPIE  Proceedings, vol.  927,  Florida, 
April  6  -  8,  1988. 

Polarimetric  microwave  remote  sensing  of  anisotropic  earth  terrain  with  strong 
fluctuation  theory  (M.  Borgeaud,  J.  A.  Kong,  and  R.  T.  Shin),  IGARSS  87  and 
URSI  Meeting,  University  of  Michigan,  Ann  Arbor,  MI,  May  18-21,  1987. 


37 


Theoretical  models  for  active  and  passive  microwave  remote  sensing  of  snow- 
covered  sea  ice  (F.  C.  Lin,  J.  A.  Kong,  and  R.  T.  Shin),  IGARSS  87  and  URSI 
Meeting,  University  of  Michigan,  Ann  Arbor,  MI,  May  18-21,  1987. 

Theoretical  models  for  polarimetric  radar  clutter  (M.  Borgeaud,  R.  T.  Shin, 
and  J.  A.  Kong),  Journal  of  Electromagnetic  Waves  and  Applications,  vol.  1,  no.  1,  pp. 
73-89,  1987. 


38 


Applications  of  EMSARS  BackgroTind  Model  to 
MMW  Vegetation  Scattering  Simulation 

C.  C.  Esti,  Y.  E.  Yang,  R.  T.  Shin,  J.  A.  Kcng 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

C.  Kohler,  T.  Nguyen,  H.  Nguyen,  J.  Ho 
US  Army  CECOM  Night  Vision  and  Electro-Optics  Directorate  (NVEOD) 
Ft.  Bel  voir,  Virginia  22060 


Abstract 

EMSARS  (Electromagnetic  Model  of  Scattering  Applied  to  Remote  Sensing)  is  an  electromagnetic  background 
scattering  model  developed  at  the  Massachusetts  Institute  of  Technology.  The  physical  characterization  of  the 
terrain  background  is  in  terms  of  layered  continuous  random  medium  and  discrete  scatterer  models.  Various 
theories  and  mathematical  techniques  such  as  the  analytical  Wave  Theory  with  the  Born  Approximation  and  the 
Radiative  Transfer  Theory  are  used  to  solve  the  integral  equations  characterizing  the  electromagnetic  scattering 
effects. 

In  this  paper,  applications  of  the  EMSARS  model  to  millimeter  wave  scattering  by  vegetation  will  be  presented. 
IVee  is  modeled  as  discrete  scatterers  with  multiple  sizes  for  the  trunk,  the  branches  and  leaves,  respectively. 
Currently,  the  leaves  of  coniferous  trees,  the  branches,  and  the  trunks  are  modeled  as  circular  cylinders.  For 
deciduous  trees,  the  leaves  are  modeled  as  thin  circular  disks.  The  physical  optics  approximation  is  used  to 
calculate  the  millimeter  wave  returns  &om  circular  disks  while  firrite  cylinder  approximation  is  used  for  circular 
cylinders.  Then,  the  branching  structure  and  their  orientation  distributions  are  applied  to  obtain  the  radar  returns 
from  vegetation.  The  trunk  is  model  as  a  circular  cylinder  with  random  surface  roughness  where  geometrical  optics 
approximation  is  used  to  calculate  the  radar  return  from  the  surface  of  trunk.  Except  for  the  case  in  which  the 
ground  is  illuminated,  we  only  consider  single  scattering  returns.  Transmission  loss  factors,  due  to  absorption  and 
scattering,  are  included  when  computing  the  contribution  from  interior  branches  and  leaves.  Theoretical  predictions 
will  be  compared  with  the  measurement  data  recently  collected  in  Bend,  Oregon  by  NVEOD  using  a  94GHz  radar. 


L  Introduction 

Electromagnetic  scattering  from  earth  terrains  is  a  very  sophisticated  physical  phenomenon.  During  the  past 
two  decades,  researchers  at  the  Massachusetts  Institute  of  Technology  have  been  engaged  in  the  development 
and  study  of  various  theoretical  models  for  the  microwave  remote  sensing  of  earth  terrain.  The  collection  of 
research  efforts  on  these  models  and  associated  computer  programs  is  called  EMSARS  (Electromagnetic  Modeling 
of  Scattering  Applied  to  Remote  Sensing).  As  summarized  in  the  monograph  “Theory  of  Microwave  Remote 
Sensing”  [1],  electromagnetics  models  have  been  developed  to  account  for  the  effects  of  absorption,  layering,  volume 
scattering  and  rough  interfaces,  and  have  been  applied  in  the  interpretation  of  active  and  passive  remote  sensing 
data  coUected  from  various  types  of  earth  terrain,  such  as  forests,  vegetation  canopies,  road  surfaces,  lava,  snow-ice 
fields,  under  environmental  and  seasonal  changes. 

In  dealing  with  radar  retiuns  from  background  media,  we  have  to  consider  the  effects  of  both  surface  and 
%oIume  scattering.  The  surface  scattering  effect  can  be  dealt  with  by  the  random  rough  surface  model  in  EMSARS. 
This  model  provides  a  simple  approach  to  account  for  returns  from  highly  reflective  surfaces  such  as  road,  sand- 
covered  terrain,  and  lake.  We  have  implemented  composte  rough  surface  model,  which  includes  the  Kirchhoff 
Approximation  (KA)  and  the  Small  Perturbational  Method  (SPM)  algorithms  for  computing  scattering  from  large- 
and  small-  scale  surface  roughnesses.  In  addition,  we  have  developed  the  extended  boundary  condition  (EBC)  and 
integral  equation  techniques  for  dealing  with  the  randomly  perturbed  periodic  rough  surface  model,  which  accounts 
for  row  structures  in  plowed  fields. 


Th«  volume  scattering  effect  is  important  for  returns  from  bare  soils,  vegetation,  forest,  snow-ice  fields,  in  which 
the  underlying  scatterers  play  important  roles.  The  volume  scattering  effect  is  also  useful  in  explaining  diurnal 
changes.  A  layered  random  medium  model  has  been  developed.  Each  layer  of  the  medium  is  described  using  both 
dctcrn’.ir.i'tic  and  statistical  paianictcrs.  For  example,  layer  thickness  can  be  used  to  model  vegetation  heights, 
snow  depth,  etc.,  and  correlation  length  and  variance  are  closely  related  to  the  shape  and  density  of  leaves  and 
water  contents.  It  has  been  shown  that  a  two-layer  random  medium  model  works  very  well  in  the  interpretation  of 
data  collected  from  snow  covered  terrain,  and  vegetation  fields  such  as  corn,  alfalfa,  milo,  and  soybeans.  Three-layer 
random  medium  model  has  also  been  developed  and  can  be  used  to  model  more  sophisticated  background. 

In  the  three-dimensional  mm  wave  radar  simulation,  it  is  sometimes  necessary  to  consider  isolated  trees  and 
shrubs  as  part  of  the  background.  In  this  case,  we  may  have  to  add  some  other  parameters,  such  as  the  size  of  the 
canopy,  when  using  the  layered  random  medium  model.  Alternatively,  a  discrete  random  scatterer  model  can  be 
applied.  This  model  provides  means  of  summing  the  contributions  from  individual  scatterers,  including  multiple 
scattering  effects,  thus  allows  us  to  work  on  a  tree-by-trec  basis.  Basically,  the  tree  trunk  can  be  treated  as  a 
deterministic  scatterer  much  like  targets,  whereas  for  branches  and  leaves  we  shall  apply  the  random  scatterer 
model  to  calculate  returns. 

One  popular  way  of  solving  the  scattering  problem  is  through  wave  theory  technique  by  expanding  the  solution 
to  the  wave  integral  equation  in  the  form  of  Born  scries.  The  first  order  term  in  the  Born  series  represents  the 
single  scattering  effects  and  the  higher  order  terms  contribute  to  multiple  scattering  effects  in  the  radar  returns. 
The  Born  series  technique  produces  solution  in  analytical  form,  but  the  expression  becomes  very  complicated  if 
multiple  scattering  effects  are  to  be  considered.  In  contrast,  the  radiative  transfer  (RT)  theory  has  been  commonly 
used  owing  to  its  simplicity  and  ease  of  handling  multiple  scattering  to  some  extent.  The  radiative  transfer  theory 
consists  of  the  radiative  transfer  equations,  which  govern  the  electromagnetic  energy  propagation  through  scattering 
media.  Various  models  have  been  developed  based  on  this  theory  [2—4].  The  conventional  RT  theory  ignores  the 
relative  phase  information  associated  with  structured  scatterers,  which  may  play  an  important  role  in  the  overall 
scattering  behavior  [5]  because  vegetation  generally  consists  of  structures  of  many  different  scale  lengths.  For 
example,  the  trunk,  primary  branches,  secondary  branches,  and  leaves,  give  rise  to  different  scattering  effects  at 
any  given  wavelength.  Vegetation  elements  of  each  scale  are  connected  to  elements  of  other  scales  in  a  fashion 
statistically  described  by  the  unique  architecture  pertaining  to  each  tree  species  group.  For  the  microwave  remote 
sensing  of  forest,  the  vegetation  structures  not  only  cause  the  separation  of  scattering  centers,  but  also  provide 
partially  coherent  scattering  by  different  scatterers  with  statistically  prescribed  relative  positions. 

In  this  paper,  a  RT  model  is  constructed  for  the  modeling  of  the  vegetation.  The  model  makes  use  of  the 
branching  model  for  vegetation  [5]  to  account  for  the  scattering  properties  of  trees.  Section  II  summarizes  the 
basics  of  radiative  transfer  theory.  In  Section  III,  the  phase  matrix  based  on  branching  model  is  formulated.  In 
Section  IV,  formula  for  the  calculation  of  range  profile  of  trees  is  given.  Section  V  presents  the  comparison  of  model 
with  experimental  data  collected  by  NVEOD  using  a  94  GHz  radar. 


U.  Radiative  Transfer  Theory 

The  canopy  of  a  tree  can  be  modeled  as  a  random  medium  of  certain  shape,  e.g.  a  cone.  The  random  medium 
is  composed  of  discrete  scatterers  embedded  in  homogeneous  background,  as  shown  in  Figure  1. 


Fig.  1  Tree  model. 


The  vector  radi&tive  transfer  equation  for  the  specific  intensity  in  the  scattering  region  is  of  the  form 
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where  the  Stokes  vector  I  contains  information  regarding  field  intensity  and  phase  relation  of  the  two  orthogonal 
polarizations  is  defined  as 
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In  (2),  the  subscripts  h  and  v  represent  the  horizontal  and  vertical  polarizations,  respectively.  The  angular  bracket 
{  )  denotes  ensemble  average  over  the  size  and  orientation  distributions  of  scatterers;  and  ij  =  y/fJo/^o  *s  the 
fiee-space  characteristic  impedance. 


The  extinction  matrix  k,  represents  the  attenuation  due  to  both  the  scattering  and  absorption,  and  can  be 
obtained  through  the  optical  theorem  with  forward  scattering  functions. 
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The  phase  matrix  P(6,  <p;  6\  <f>')  characterizes  the  scattering  of  the  Stokes  vector  from  (6',  <j>')  direction  into 
(6,  direction.  It  can  be  formulated  in  terms  of  scattering  functions  of  the  randomly  distributed  discrete  scatterers. 
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where  no  is  number  density  of  scatterers,  k  is  the  free  space  wave  number,  fab  is  the  element  of  scattering  matrix 
with  incident  field  having  polarization  b,  and  scattered  field  has  polarization  a.  a,  b  can  be  horizontal  or  vertical 
polarization. 

The  scattering  matrix  which  relates  the  incident  electric  field  to  the  scattered  electric  field  is  defined  as 
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where  e''‘^/r  is  the  spherical  wave  transformation. 


Along  with  the  boundary  conditions  that  require  the  Stokes  vector  to  be  continuous,  we  can  solve  the  radiative 
transfer  equations  iteratively  for  the  backscattering  radar  cross-section  [1]. 


rig.  2  Cylinder  cluster. 


in.  Scattering  lYinction  for  Clustered  Structure 


To  take  into  account  the  partially  coherent  effect  due  to  clustered  vegetation  structures,  we  formulate  the 
phase  matrix  based  on  the  branching  model  for  vegetation.  Take  a  cylinder  cluster  for  example.  For  a  cluster  that 
has  one  center  cylinder  and  N  branching  cylinders  as  shown  in  Figure  2,  the  total  backscattering  function  fap  is 

N 

fap  =  foafi  +  ^  fnap  (6) 

fl  =  l 

where  foap  is  the  a/3-th  clement  of  the  scattering  matrix  for  the  center  cylinder,  and  fnap  is  a^-th.  element  of  the 
scattering  matrix  for  the  n-th  branching  cylinder.  a,P,-f,6  represent  horizontal  or  vertical  polarization. 

Assuming  all  the  branching  cylinders  are  identical  and  independent  of  each  other,  and  their  relative  position 
is  independent  of  their  scattering  properties  such  as  size  and  orientation,  the  correlation  of  /  is 
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The  relative  phase  of  the  n-th  branch  with  respect  to  the  center  cylinder  is  defined  as 
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where  fci  and  k,  arc  the  incident  and  scattered  wave  vectors,  respectively.  ¥„  is  the  location  of  the  n-th  branching 
cylinder  relative  to  the  center  cylinder. 

In  (6),  the  third  and  fourth  terms  are  the  coherent  terms.  It  can  be  seen  that  the  incoherent  approximation 
is  valid  when  the  average  of  the  random  phase  factor  {e*^^)  is  so  small  that  the  coherent  terms  are  negligible 
as  compared  wth  the  incoherent  terms.  However,  this  is  not  true  for  vegetation  structures  with  scale  lengths 
comparable  to  the  wavelength,  as  was  demonstrated  with  a  two-scale  clustcr[5]. 

For  the  calculation  of  scattered  fields  of  different  cluster  elements,  dielectric  cylinders  can  be  used  to  model 
trunks,  branches,  and  coniferous  leaves  [6,7].  Leaves  of  deciduous  tree  can  be  modelled  as  disks  [6,3].  The 
physical  optics  approximation  is  used  for  the  calculation  of  scattering  from  disks.  The  truncated  infinite  cylinder 
approximation  [9]  is  employed  in  this  study  to  calculate  scattered  fields  from  cylinders. 


rV.  Calculation  of  Range  Profile 


In  cases  when  the  elevation  angle  is  small,  ignoring  the  ground  effect,  the  first-order  iterative  solution  of  the 
RT  equation,  which  specifies  the  scattering  Stokes  vector  I,  in  terms  of  incident  Stokes  vector  7©  is  as  follows: 
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where  D(P{6,  d)>')  is  the  4x4  diagonal  matrix  with  )  as  its  tt-th  elements.  Here  Pi  is  the  t-th  eigenvalue 

of  Kf  s*  is  in  the  direction  of  the  propagation  and  li ,  Ij  are  the  boundaries  of  the  random  medium. 

Equation  (8)  represents  the  single  scattering  by  the  discrete  scatterers  inside  the  random  medium.  D  is  the 
attenuation  factor  and  P  can  be  considered  as  the  scattering  intenaty.  M  is  Mueller  matrix.  Mu  and  M32 
radar  cross-sections  per  unit  area  for  E—  and  V—  polarizations,  respectively. 

If  the  random  medium  is  azimuthally  symmetrical  around  the  propagation  direction,  the  backscattering  radar 
cross  section  per  unit  area  for  the  co-polarized  returns  can  be  further  simplified  as 


(9) 


where  Ka  =  •  2/m  <  /»«  >  is  the  attenuation  constant  for  polarisation  a.  a  could  be  h  or  v. 

Therefore  the  radar  cross-section  per  unit  length  in  the  propagation  direction  for  polarization  a  can  be  expressed 

as 
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In  the  computation  of  trunk  scattering  contributions,  surface  roughness  is  an  important  factor.  At  94  GHz 
the  wavelength  is  much  smaller  than  the  radius  of  trunk  and  thus  the  surface  roughness,  no  matter  how  small,  can 
still  be  comparable  to  the  wavelength.  The  radar  cross-section  is  calculated  by  breaking  up  the  trunk  surface  into 
small  pieces  and  summing  up  the  response  from  individual  pieces.  The  radar  cross-section  of  each  small  surface  is 
calculated  by  multiplying  the  area  with  the  scattering  coefficient.  Thus  the  backscattering  radar  cross-section  per 
unit  length  is 
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where7  is  the  scattering  coefficient  for  the  reflected  wave  from  the  trunk  surface.  It  is  a  function  of  incident  angle 
and  the  local  curvature  of  trunk. 

To  calculate  the  range  profile,  we  need  to  convolve  the  radar  cross-  section  per  unit  length  (from  Equations  10 
and  11)  with  the  radar  response  of  the  point  target.  This  impulse  response  is  determined  by  the  bandwidth  and 
the  filter  function  of  the  radar  system. 


V.  Comparison  with  Experimental  Data  and  Sensitivity  Study 


Theoretical  prediction  and  computer  simulation  of  mm  wave  radar  range  profile  are  based  on  the  experiment 
conducted  by  NVEOD  in  Bend,  Oregon  with  a  truck-mounted  94  GHz  FM  CW  radar.  A  high  range  resolution  of 
about  l(one)  foot  is  achieved  by  stepping  the  radar  frequency  through  the  470MHz  bandwidth.  With  frequency 
step  size  of  lOMHz  and  total  of  48  steps,  the  system  provides  range  ambiguity  to  about  15  meters. 

The  NVEOD  radar  was  aimed  at  two  isolated  coniferous  trees,  as  shown  in  Figure  3.  The  elevation  angles 
are  around  5  degrees  (looking  up);  therefore,  we  can  safely  neglect  the  ground  effect  in  our  calculation.  Within 
the  cluster,  the  contribution  to  scattered  field  at  94  GHz  mainly  comes  from  the  leaves.  Thus  the  branches  were 
omitted  from  our  calculation.  The  scattered  field  from  trunk  has  also  been  computed  separately.  The  small  needles 
of  the  tree  leaves  inside  the  canopy  are  modeled  as  randomly  oriented  dielectric  cylinders. 


The  canopy  of  both  trees  is  considered  as  a  random  medium  of  cone  shape.  The  radius  and  length  of  the 
dielectric  cylinders,  of  which  the  random  medium  is  composed,  are  assumed  to  be  0.1cm  and  2cm,  respectively. 
The  number  density  of  leaves  is  12000/m*  and  dielectric  constant  is  7  +  tl.5.  Using  the  photographs  taken  at 
the  scene,  in  which  a  yaidstick  is  inserted  in  front  of  each  tree,  we  can  estimate  the  tree  height.  The  first  one  is 
estimated  at  9  meters  high  and  crown  length  of  7  meters  (calculated  from  the  the  first  live  branch  to  the  top  of 
tree)  The  radii  of  top  and  bottom  circles  are  0.7Sm  and  3.0m,  respectively.  As  for  the  second  tree,  the  height  of 
canopy  is  estimated  at  8  meters  and  crown  length  of  6.5  meters.  The  upper  and  lower  circle  radii  are  0.2m  and 
2.7m,  respectively. 

The  comparisons  between  the  predicted  and  measured  range  profiles  from  the  crown  of  both  trees  are  given  in 
Figure  4.  The  measurement  and  prediction  results  are  given  for  the  HH  polarization.  Because  we  did  not  have  the 
ground  truth  information  for  nearby  terrains,  the  comparison  is  limited  to  the  extent  of  the  trees. 


Tree  f  1  without  trunt; 


Tree  f2  without  trunk 


Fig.  4  Comparison  of  range  profile  from  experimental  data  and  theoretical  simulation  with  crown  only. 
Solid  curve  is  simiJatcd  result,  dotted  emve  is  data  measured  with  1  foot  resolution. 


The  scattering  fcom  trunk  is  calculated  vrith  the  trunk  surface  having  tms  height  0.2cm  and  correlation  length 
1cm.  The  radius  of  first  tree  is  estimated  to  be  0.2m  and  the  radius  of  second  tree  to  be  0.15m.  Using  the 
geometrical  optics  approximation  and  assuming  the  slope  of  surface  roughness  is  Gaussian  distributed,  the  range 
profile  calculated  for  the  whole  tree  is  shown  in  Figure  5.  It  is  shown  that  the  scattering  of  trunk  will  add  a 
small  peak  on  the  range  profile.  The  agreement  between  the  theoretical  prediction  and  the  experimental  data  are 
generally  good  for  the  second  tree.  As  for  the  first  tree,  the  theoretical  simulation  is  about  2  to  3  dB  lower  than  the 
data.  The  discrepancy  may  due  to  that  the  multiple  scattering  effect  is  ignored  in  our  calculation,  or  that  perhaps 
the  finite  cylinder  approximation  is  not  accurate  enough. 


Tree  with  trunk 


Tree  §2  with  trunk 


Fig.  5  Comparison  of  range  profile  from  experimental  data  and  theoretical  simulation  with  both  crown 
and  trunk.  Solid  curve  is  rimulated  result,  dotted  curve  is  data  measured  with  1  foot  resolution. 


As  part  of  this  study,  we  are  also  interested  in  the  quantitative  dependence  of  the  range  profile  on  selective 
parameters  of  the  tree  model.  These  parameters  tend  to  vary  from  trees  to  trees,  and  hence  should  provide  more 
insights  to  the  modeling  effort.  The  following  is  a  sensitivity  study  ba.sed  on  the  tree  size,  number  density  of  leaves 
and  the  scattenng/attcnuation  ratio  is  presented.  The  variations  ate  limited  to  the  leaves,  so  the  contribution  from 
trunk  ate  not  shown  in  the  comparison. 

Figure  6(a)  shows  the  range  profile  of  two  trees  with  different  sizes.  One  has  radii  of  top  and  bottom  circles 
at  0.7Sm  and  3.0m,  respectively,  and  is  9m  high.  Another  has  radii  of  top  and  bottom  circles  at  1.125m  and 
4.5m,  respectively,  and  is  13.5m  high.  It  was  shown  that  the  peak  radar  cross-section  of  the  second  tree  is  about 
2  to  3  dB  higher  and  the  range  of  significant  radar  return  is  about  60%  wider.  This  matches  the  increase  in 
the  tree  size.  The  simulation  on  different  number  density  of  the  leaves  is  shown  in  Fig.  6(b).  As  the  number 
density  increases,  the  radar  cross-section  from  the  front  part  of  range  profile  increase,  but  it  also  decrease  faster 
with  increasing  range  at  the  backside  of  the  tree.  This  is  because  as  the  number  density  of  leaves  increase,  the 
scattering  and  attenuation  both  increase.  The  radar  return  from  the  front  part  of  tree  would  increase  but  the 
return  from  the  back  side  of  tree  would  decrease  thanks  to  higher  attenuation.  A  few  simplifications  in  our  model 
may  have  affected  the  accuracy.  Firstly,  the  scattering  and  attenuation  from  the  leaves  are  calculated  from  the 
finite  cylinder  approximation.  Also,  the  multiple  scattering  effect  is  ignored.  These  two  likely  impact  upon  the 
scattering/attenuation  ratio.  We  performed  a  simulation  in  which  the  scattering  coefficient  is  kept  intact  but  the 
attenuation  coefficient  is  decreased.  As  the  results  shown  in  fig.  6(c),  the  front  part  of  the  range  profile  does 
not  change,  as  was  expected  since  the  attenuation  docs  not  affect  this  part  significantly.  But  for  the  backside  of 
the  tree,  the  lower  attenuation,  the  bigger  the  radar  cross-section  and  the  wider  the  range  with  appreciated  radar 
return.  Therefore  in  the  model,  the  correct  estimate  of  attenuation  is  very  important. 


Range  (meter) 


Fig.  fl(a)  Simulated  rarrge  profile  of  trees  with  different  tree  size,  r  b  the  bottom  circle  radius  of  canopy, 
h  b  tree  height. 


Tree  with  different  density  of  needles 


Fig.  6(b)  Simulated  range  profile  of  trees  with  different  number  density  of  leaves,  n  is  the  number  density 


of  leaves  per  m*. 


Tree  with  different  ottenuation 


Ronge  (meter) 


Fig.  6(c)  simulated  range  profile  of  trees  with  different  attenuation  and  fixed  scattering  amplitude  of 
unit  vohuuc. 


VI.  Conclusion 


In  this  pp.prT.  an  clcclrorr.acnctic  scattering  model  of  vegetation  based  on  tb.e  radiative  transfer  ti.c^r;  n 

presented.  The  phase  matrix  is  derived  with  the  coherent  effects  from  the  clustered  structures  of  leaves.  Using  this 
simple  model,  theoretical  prediction  has  been  derived  for  the  range  profile  simulation  and  the  RCS  prediction  of 
single  trees.  The  slight  discrepancy  with  experimental  data  may  be  due  to  the  approximation  in  calculating  the 
attenuation  coefficient.  A  simulation  is  also  presented  to  examine  the  sensitivity  of  different  parameters. 
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where  Q,  uiid  Qi  are  2  by  2  matrices  corresponding  to  vectors  Q,  and  Qx,  From  (8),  it  can  be  sliown  that 

respectively,  through  the  mapping  defined  by  (2).  ‘  *  ’  denotes  complex  conjugate. 

Hence,  T  =  c''Q[n-'Q‘,(QtQ‘,) 
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When  aziinutlially  symmetric  distributed  targets  arc  not  available  in  the  im-  (a) 

aged  scene,  calibration  by  particular  solutions  may  prove  to  be  effective  if  the  i 

system  has  good  cross-talk  isolation  like  the  radar  described  by  (51)  and  (52). 

Particular  solutions  of  the  case  considered  are  given  by  Figure  2.  Cu-pularization  signatures  of  a  distributed  target  with  az- 

imutiial  symmetry:  (a)  original. 


Ca/i6ra(>on  of  Potarimetric  Rada 


os  of  a  distributed  target  with  az- 
calibrated  by  the  porticidar  solution, 
calibrated  by  the  particular  solution, 
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SIMULATION  OF  FOREST  BACKSCATTER  AS  A  FUNCTION  OF  FOREST  AND  GROUND  PARAMETERS 
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MIT,  Cambridge  MA02139  .USA  /  Tel:  (613)253-5625  -  Fas:  (613)253-0987 


ABSTRACT 

To  develop  an  algorithm  for  retrieving  forest  biomass  using 
SAR  data,  a  study  is  performed,  which  includes:  1) 
observations  on  the  relation  between  backscatter  and  forest 
biomass,  2)  interpretation  of  the  observations  by  modelling 
and  3)  simulation  of  radar  response  for  a  wide  range  of 
forest  and  ground  conditions.  This  paper  presents 
preliminary  results  of  P-band  backscatter  simulations  as 
(unctions  of  radar  incidence  angle,  ground  slope, 
undergrowth  vegetation  state,  vegetation  moisture  content, 
and  crown  structure. 

Keywords:  P-band  SAR  backscatter,  forest  biomass, 
modellmg,  simulations,  forest  and  ground  parameters. 

1.  INTRODUCTION 

The  analysis  of  JPL  SAR  data  over  the  Landes  forest 
(MAESTRO-1  campaign)  revealed  strong  correlation 
between  L-  and  especially  P-band  backscatter  o*  and  the 
pine  forest  parameters  including  trunk  biomass  {1]. 
However,  it  was  not  obvious  to  conclude  about  the 
particular  effect  of  a  given  forest  parameter  on  since 
those  parameters  were  found  strongly  correlated.  To 
explain  the  physical  link  between  the  P-band  backscatter 
and  biomass,  a  new  polarimetric  backscatlering  model, 
based  on  the  vector  radiative  transfer  theory  including  a 
branching  model,  was  developed  and  successfully  validated 
with  the  measurements  for  P-band  data  [2]. 

2.  RELATING  P-BAND  ff*  TO  BIOMASS 

Fig.  1  presents  the  main  scattering  contributions  to  the 
P-band  as  functions  of  pine  forest  age.  The  major 
contribution  to  the  HH  return  is  the  crown  scattering  for 
young  forest  and  trunk-ground  double  scattering  for 
mature  forest.  The  VV  and  HV  returns  are  dominated  by 
crown  scattering,  which  comes  mainly  from  the  primary 
branches,  modelled  as  dielectric  cylinders.  Thus,  W  and 
HV  o*  are  directly  linked  to  diameter,  length,  fractional 
volume  and  moisture  content  of  primary  branches.  For  the 
available  age  range  (8  to  46  yrs),  the  primary  branch  mean 
length  and  radius  were  found  to  vary  from  1.2  to  2.5m  and 
1.2  to  3cm  respectively.  Therefore  o*  increase  with  branch 
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fresh  biomass  (total  volume  x  dry  density  *■  water  mass)  is 
obvious. 

According  to  this  result,  high  correlation  observed 
between  HV  and  VV  returns  and  trunk  biomass  {!]  is 
indirect  and  arises  from  a  correlation  between  biomass  in 
each  tree  parts,  i.e.  the  crown  and  trunk.  Fig.2  shows  the 
trunk,  branch  and  total  above-ground  biomass  at  the  stand 
level  as  a  function  of  forest  age,  derived  using  an  aliometric 
model  [3],  which  estimates  tree  biomass  by  parts  using  the 
trunk  section  (tree  basal  area).  Biomass  by  parts  increase 
quasi-linearly  with  age  and  are  highly  intercorrelated 
through  proportional  biomass  production  relative  to  total 
above-ground  biomass. 


Figure  1 :  P-band  main  scanering  contributlont  to  HH,  HV 
and  W  backscatter  as  a  function  of  forest  age 
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The  above  results  provide  a  significant  improvement 
of  our  knowledge  on  the  scattering  mechanisms  of  forest 
canopy.  Notice  that  the  VV  and  HV  returns  should  not  be 
directly  affected  by  environmental  conditions  in  bottom 
forest  layers  (bush  density,  soil  moisture,  ground  slope,...) 
compared  to  HH  return,  but  should  be  affected  by  crown 
structure  which  is  highly  species  dependent.  These  points 
are  addressed  in  the  following. 

3.  SIMULATION  STUDY 

3.1  Changing  Incidence  angle 

As  P-band  yields  best  sensitivity  to  biomass,  we  study 
the  effect  of  radar  incidence  angle  d  on  tf*  relative  to  forest 
age  (hence  forest  biomass).  Rg  3.  shows  or*  as  a  function  of 
6  for  pine  stands  at  different  ages.  The  peak  values  have 
been  found  in  the  range  30*-50°  depending  on  the 
polarization.  This  also  should  depend  clearly  upon  the 
crown  structure  at  HV  and  W.  Here,  an  incidence  of  40*  is 
particularly  adapted,  for  both  o*  sensitivity  to  forest  age 
(thus  biomass)  and  for  the  interesting  high  HV  <»*  values. 

3.2  Changing  forest  and  ground  parameters 
Effect  of  ground  slone 

Fig.4  presents  o*  as  a  function  of  local  ground  slope 
facing  the  radar.  As  the  local  slope  affects  mainly  the 
ground,  ground-trunk,  crown-ground  scattering  terms,  it  is 
expected  that  W  and  HV  will  not  be  changed  with  varying 
local  slope,  whereas  HH  return  should  be  significantly 
affected.  This  is  confirmed  in  Fig  4,  where  the  HH  return  is 
fastly  decreasing  with  local  slope  for  old  forests  where  the 
trunk-ground  term  dominates,  and  is  approaching  to  the 
crown  scattering  term.  No  significant  changes  occur  for  HV 
and  VV. 

Effect  of  undergrowth  vegetation 

To  account  for  varying  understory  conditions,  o*  arc 
computed  for  a  canopy  without  understory  vegetation  and 
with  a  relatively  dense  bush  layer.  Since  HV  and  VV  return 
results  from  the  crown  scattering,  no  change  is  observed. 
The  decrease  of  o*  for  HH,  caused  by  the  attenuation  of 
trunk-ground  term,  docs  not  exceed  1  dQ  with  the  addition 
of  the  bush  layer  (height  =  Im,  LAI  «  l.O) 

Effect  of  vegetation  moisture  condition 

To  account  for  variations  in  vegetation  moisture 
conditions  related  to  seasonal  nr  diurnal  effects,  the 
backscattcr  coefficients  arc  calculated  as  a  function  of 
vegetation  gravimetric  moisture  content  from  55%  to  65%, 
HH  is  not  significantly  sensitive  to  this  variation,  whereas 
VV  and  HV  returns  arc  more  affected,  the  variation  in 
backscattcr  coefficient  is  of  the  order  of  IdB. 


Effect  of  crown  itnictiire 

Relations  between  the  0*  and  forest  biomass  should 
be  species  dependent  through  the  dependence  on  crown 
structure.  As  P-band  HV  and  VV  come  mainly  from 
primary  branches,  a  major  structural  parameter  to  consider 
is  the  mean  elevation  angle  0  of  the  primary  branch 
orientation.  Fig.5  presents  the  simulated  0*  as  a  function  of 
0,  for  different  forest  ages  observed  at  45*  incidence.  At 
HV  and  W,  the  maximum  sensitivity  to  forest  age  is  found 
respectively  at  60*  (0  value  for  maritime  pine)  and  45* 
(plagiophile  species).  For  low  0  (crcctophile  species)  there 
is  no  significant  sensitivity  to  forest  age.  This  is  also  the 
case  for  W  for  high  0  (planophile  species).  However,  for  a 
majority  of  tree  species,  0  ranges  from  45*  to  90*.  In  this 
range,  the  sensitivity  observed  for  forest  from  6  to  46  years 
is  from  7  to  10  dB  in  this  simulation  case.  For  W, 
incidence  angles  complementary  to  0  should  yield  the  best 
sensitivity  of  o*  to  forest  biomass.  For  HH  returns,  no 
effect  of  0  on  the  0*  sensitivity  to  biomass  is  observed,  as 
the  return  results  mainly  from  trunk-ground  interaction. 
HH  is  therefore  less  sensitive  to  species  diversity. 

4.  DISCUSSION  AND  CONCLUDING  REMARKS 

At  present,  the  modelling  and  simulation  studies  were 
performed  at  P-band,  which  was  found  optimal  for  retrieval 
of  forest  biomass.  HH  return  was  found  physically  related 
to  both  trunk  and  crowTj  biomass,  whereas  W  and 
especially  HV  returns  were  found  tightly  linked  to  crown 
biomass  and  unaffected  directly  by  environmental 
parameters  such  as  soil  moisture,  ground  local  slope  and 
presence  of  an  understory  bush  layer.  Therefore,  it  is 
possible  to  derive  total  above-ground  biomass  from  P-band 
HV  or  VV  SAR  data,  if  the  relations  between  biomass 
found  in  different  tree  parts  arc  known.  On  the  other  hand; 
HH  return  should  be  less  species  dependent  as 
contribution  of  the  crown  is  lower  compared  to  trunk- 
ground  interaction,  but  is  more  affected  by  environment 
conditions  stated  above. 
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Figure  3  ;  Effect  of  Incidence  artgie  on  pine  backscatter  at  P>band  HH,W  and  HV 


Figure  4  ;  effect  of  ground  slope  on  pine  backscatter  at  HH,  HV  and  W  pol. 
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Abstract —  A  branching  model  is  proposed  for  the  remote 
sensing  of  vegetation.  The  frequency  and  angular  responses  of 
a  two>scale  cylinder  cluster  are  calculated  to  demonstrate  the 
signihcance  of  vegetation  architecture.  The  results  indicate  that 
it  is  necessary  for  theoretical  models  to  take  into  account  the 
architecture  of  vegetation  which  plays  an  important  role  in 
determining  the  observed  coherent  effects.  A  two-scale  branching 
model  is  implemented  for  soybean  with  its  internal  structure  and 
the  resulting  clustering  effects  considered.  Furthermore,  at  the 
scale  of  soybean  fields,  the  relative  location  of  soybean  plants  is 
described  by  a  pair  of  distribution  function.  The  polarimetric 
backscattering  coefficients  are  obtained  in  terms  of  the  scat¬ 
tering  properties  of  soybean  plants  and  the  pair  distribution 
function.  Theoretical  backscattering  coefficients  evaluated  using 
the  hole-correction  pair  distribution  are  in  good  agreement  with 
extensive  data  collected  from  soy  bean  fields.  Compared  with  the 
independent-scatterer  pair  distribution,  it  is  found  that  the  hole- 
correction  approximation,  which  prevents  two  soybean  plants 
from  overlapping  each  other,  is  more  realistic  and  improves  the 
agreement  between  the  model  calculation  and  experimental  data 
near  normal  incidence.  Extension  to  a  multiscale  branching  model 
can  be  achieved  by  recursion  of  the  two-scale  modeling  approach 
presented  in  this  paper. 


I.  Introduction 

There  has  been  growing  interest  in  the  investigation  of 
vegetation  using  polarimetric  remote  sensing  measure¬ 
ments.  In  the  past,  theoretical  models,  including  continuous 
random  medium  and  discrete  scatterer  approaches,  have  been 
applied  for  the  electromagnetic  modeling  of  vegetation  {!]. 

Continuous  random  medium  approach  utilizes  a  random 
process  described  by  a  correlation  function  to  represent  the 
permittivity  distribution  of  terrain  medium.  In  this  approach, 
the  effective  permittiv  ity  is  evaluated  to  account  for  the  attenu¬ 
ation  and  phase  delay  of  w  ave  propagation  inside  the  medium. 
Subsequently,  the  polarimetric  backscatter  coefficients  are 
calculated  and  used  to  interpret  the  scattering  from  various 
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geophysical  media  [1],  [2].  This  approach  has  been  applied  for 
the  interpretation  of  polarimetric  backscattering  measurement 
from  vegetation  with  seleaed  conelation  function  parameters 
[2] 

Discrete  scatterer  approach  models  the  medium  as  a  random 
collection  of  discrete  scatterers  with  different  size  and  shapes. 
The  relative  location  of  scatterers  is  described  by  a  pair 
distribution  function,  and  the  backscattering  coefficients  and 
effective  propagation  constants  have  been  evaluated  under 
various  approximations  [1].  With  the  application  of  discrete 
scatterer  theory  to  vegetation,  trunk  and  branches  have  been 
modeled  as  dielectric  cylinders  [3],  [4],  and  the  internal 
and  scattered  fields  were  estimated  by  invoking  the  finite 
cylinder  approximation  [5].  To  take  into  account  the  leaves 
of  deciduous  trees,  leaves  were  modeled  as  dielearic  disks 
(3],  (6],  and  the  internal  and  scattered  fields  were  calculated 
by  using  the  generalized  Rayleigh  Cans  approximation  [3]  and 
the  physical  optics  approximation  (7].  In  this  discrete  scatterer 
approach,  the  radiative  transfer  theory  was  applied  for  the 
backscattering  cross  section  of  randomly  embedded  scatterers 
18).  19). 

Durden  ei  al.  [4]  modeled  the  forest  as  the  following  three 
layers:  the  crown  layer  contained  randomly  oriented  cylinders 
representing  branches;  the  trunk  layer  below  the  crown  region 
contained  nearly  vertical  cylinders  representing  tree  trunks; 
and  the  bottom  layer  was  a  rough  surface  representing  the 
forest  floor.  A  radiative  transfer  approach  was  also  used  to 
evaluate  the  effective  propagation  constant  and  the  transfer 
of  the  Stokes  matrix  when  the  waves  traverse  through  the 
crown  or  trunk  layer.  Because  the  crown  and  trunk  layers 
are  assumed  to  be  uniform,  the  propagation  effects  were  then 
represented  by  a  transmission  matrix.  Hence  the  backscatiered 
Stokes  matrix  was  approximated  as  the  sum  of  five  terms 
corresponding  to  five  different  scattering  mechanisms.  The 
calculated  polarization  signatures  were  shown  to  be  consistent 
with  the  observed  polarization  signatures  from  forested  areas. 

However,  the  form  or  structure  of  vegetation  has  not  been 
addressed  in  these  models.  In  other  words,  none  of  the  current 
models  takes  ito  account  the  differing  growth  pattern  of 
vegetation  or  v..e  architectures.  In  the  book  Fractal  Geometry 
of  Nature  (10),  Mandelbrot  utilized  the  “pipe  model."  w  hich 
describes  trees  as  bundles  of  nonbranching  vessels  of  fixed 
diameter,  connecting  roots  to  the  leaves,  to  simulate  botanical 
trees.  The  resulting  vegetation  patterns  are  very  appealing 
due  to  the  consideration  of  vegetation  architecture.  Besides 
its  significant  impact  on  the  vivid  simulation  of  vegetation 
patterns,  architecture  or  structure  is  also  very  important  in  the 
computation  of  electromagnetic  scattering  and  radiation.  For 
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Fig.  1.  Configuration  of  a  two-scale  cluster  with  a  center  cylinder  and 
perpendicularly  attached  cylinders. 

example,  the  structure  of  antenna  array  has  significant  effects 
on  the  radiation  pattern  and  polarizations  of  fields.  Further¬ 
more,  by  including  the  concept  of  multiple  scale  structure  into 
the  antenna,  such  as  the  log-periodic  array,  the  antenna  array 
would  cover  a  broad  frequency  band  [11]. 

In  general,  vegetation  consists  of  structures  of  many  length 
scales,  and  each  type  of  vegetation  has  its  own  particular 
structure  and  form  [12].  Crops  have  simpler  structures  and 
fewer  scales,  whereas  trees  have  more  scales  and  are  more 
complicated  in  form.  By  making  an  analogy  between  vege¬ 
tation  and  antenna  anays.  it  can  be  anticipated  that  the 
information  pertaining  to  the  structure  and  biomass  of  veg¬ 
etation  is  embedded  in  the  frequency,  angle,  and  polarization 
responses  of  electromagnetic  backscatter.  Hence,  structure, 
playing  an  imponant  role  in  determining  the  frequency  and 
angular  responses  of  radiation  and  scattering,  should  not  be 
ignored  in  the  electromagnetic  modeling  of  vegetation. 

In  this  paper,  a  branching  model  is  proposed  for  the  remote 
sensing  of  vegetation.  Section  11  illustrates  the  importance 
of  coherence  effects  observed  in  the  frequency  and  angu¬ 
lar  electromagnetic  responses  of  a  vegetaiion-Iike  cluster. 
Section  III  presents  the  modeling  approach  of  a  two-scale 
branching  model  for  soybean.  In  Section  IV,  theoretical  results 
are  compared  with  experimental  soybean  data  to  validate 
the  proposed  model.  The  approach  of  a  multiscale  branching 
model  and  related  issues  for  the  modeling  of  vegetation  are 
discussed  in  Section  V. 


II.  Branching  Model  for  a  Two-Scale  Cluster 

This  section  is  designed  to  illustrate  the  approach  of  branch¬ 
ing  model  and  its  significance.  The  approach  will  be  described 
for  the  scattering  of  a  scaler  field  from  a  two-scale  cluster. 
Examples  will  be  given  to  demonstrate  the  effects  of  coherence 
on  the  frequency  and  angular  responses  of  scattered  power 
which  is  the  quantity  of  interest. 

Fig.  1  illustrates  a  typical  two-scale  vegetation  cluster  with 
some  smaller  cylinders  attached  to  a  large  cylinder.  To  simplify 
the  computation  all  the  branch  elements  will  be  assumed  to  be 
mutually  independent  and  statistically  identical,  independent 
of  their  growth  locations,  and  independent  of  the  center  trunk. 

Scattering  from  such  a  cluster  can  be  approximated  by  a 


sum  of  single  scattering  terms. 

5  =  So-i-^Sjexp(i^j)  (1) 

j=i 

where  So  is  the  return  from  the  center  trunk  and  Sj  represents 
the  response  of  branch  element  j.  N  represents  the  number 
of  branches.  A  phase  factor,  (j,  is  used  to  account  for  the 
additional  phase  delay  with  respect  to  the  reference  point  “O", 
and  is  given  as 

ij  =  (k,-k.)  Tj  (2) 

Here  ry  represents  the  branching  location  of  element  j.  Defi¬ 
nitions  of  wave  vectors  are  given  in  Appendix  A. 

Tlie  scattered  power  can  then  be  expressed  as 

\  S’ 

|S|"  =  |So|*  +  Yi  «xp(t^,)] 

J=1 

s 

3=1 

where  denotes  complex  conjugate.  The  third  and  fourth 
terms,  not  having  been  accounted  for  in  cunent  radiative 
transfer  models  of  vegetation,  take  into  consideration  the  phase 
interference  between  scattering  elements  and  affect  the  overall 
scattering  pattern  of  the  cluster. 

Employing  the  assumptions  that  all  the  branch  elements 
are  statistically  identical,  mutually  independent,  uniformly 
distributed  on  the  trunk,  and  independent  of  the  center  trunk, 
yields  the  ensemble  average  of  scattered  power, 

(iSl^)  =  {\Sof)^N(\S,f)  +  2N  Re[(S5>(S,)(exp(t^)>] 
N[N  -  l)(S*)(5,>(exp(-iC,)>(exp(iC,)>.  (4) 

In  current  radiative  transfer  models  of  vegetation,  the  third 
and  fourth  terms  in  (4)  have  not  been  considered  in  the 
theoretical  models  due  to  the  fact  that  the  total  phase  matrix 
is  calculated  as  an  incoherent  sum  of  the  phase  matrices  of 
individual  elements.  Incoherent  approximation  is  valid  when 
the  average  of  a  random  phase  factor,  (exp(i(j)),  is  so  small 
that  the  magnitudes  of  the  last  two  terms  are  negligible  as 
compared  with  the  first  two  terms.  However,  that  may  not 
be  valid  at  lower  frequencies  and  various  incident  angles  for 
such  a  cluster  like  that  shown  in  Fig.  1.  Moreover,  it  can  be 
expected  that  when  AT  is  a  large  number,  the  third  and  fourth 
terms  in  (4)  will  become  significant. 

Figs.  2  and  3  illustrate  the  radar  cross  section  in  the 
backscattering  direction  versus  angles  of  incidence  and 
frequencies  for  the  two-scale  cluster  shown  in  Fig.  1.  The 
center  cylinder  in  Fig.  1  is  aligned  with  the  z-axis,  i.e., 
fj  =  zzj.  All  the  branches  are  located  on  y-z  plane,  equally 
probably  distributed  on  both  sides  of  th^  center  cylinder, 
and  perpendicular  to  the  center  cylinder,  k,  and  k,  are  on 
T-z  plane.  The  finite  cylinder  approximation  (3),  [9]  is  used 
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Fig  2.  HH  and  VV  backscaliered  responses  from  a  ruo-scale  cluster  versus 
incident  angle  ai  frequency  5.3  GHe.  The  dielectric  constant  of  cylinders  is 
(29.0  *  (9.0).  The  length  and  radius  of  center  cylinder  are  10  cm  and  0.1  cm. 
The  length  and  radius  of  branches  are  2  cm  and  0.1  cm.  .V  =  10. 


OUSTER  or  CYLINDERS 


Fig  3.  HH  and  VV  backscaliered  responses  from  a  two-scale  cluster  versus 
frequency  at  r*  =  00  degree  incidence.  The  dielectric  constant  of  cylinders  is 
(29  0  +  i9  0)  The  length  and  radius  of  center  cylinder  are  10  cm  and  0.1  cm. 
The  length  and  radius  of  branches  are  2  cm  and  0.1  cm.  .V  =  10. 


to  calculate  the  scattering  matrix  of  center  and  branching 
cylinders. 

By  assuming  the  branching  elements  are  randomly  and 
independent  positioned,  the  average  of  random  phase  factor 
can  be  carried  out  and  is  given  as 


d 

(exp(i^J)  =  j  J  exp 
0 


exp{i2kt,d)  —  1 

i2k,id 


(5) 


where  d  is  the  length  of  the  center  cylinder. 

Fig.  2  plots  the  HH  and  W  backscatter  responses  versus 
the  angle  of  incidence.  Note  that  the  incoherent  sum  for  HH 
(the  first  two  terms  in  (4)]  is  a  flat  curve  showing  no  angular 
variation.  That  is  due  to  the  fact  that  HH  response  is  dominated 
by  the  branching  cylinders  and  the  response  from  each  one 
has  no  angular  dependence.  However,  the  phase  interference 
which  exhibits  the  structure  of  this  cluster  docs  present  itself 


in  the  coherent  sum.  and  results  in  an  observable  10-dB 
difference  between  the  full  coherent  sum  and  the  incoherent 
sum  at  incident  angles  near  90°  where  the  phase  angle  is 
small  and  the  fields  are  essentially  summed  coherently.  The 
10-dB  increase  observed  is  due  to  the  number  of  branches, 
N,  being  ten  for  this  case.  If  N  is  further  increased,  the 
difference  will  also  increase.  In  the  case  studied  where  all 
the  branching  elements  are  perpendicular  to  V  polarization 
and  have  small  W  response,  the  total  W  backscatter  is 
dominated  by  the  center  cylinder.  Consequently,  there  is  no 
significant  difference  between  coherent  and  incoherent  sums. 
However,  if  the  branches  are  not  perpendicular  to  the  trunk,  it 
can  be  expected  that  the  difference  tetween  the  coherent  and 
incoherent  sums  of  W  will  be  significant. 

Fig.  3  plots  the  HH  and  W  backscatter  versus  the  frequenc. 
at  the  incident  angle  6  =  60°.  Note  that  the  full  and  incoherent 
sums  of  HH  show  an  approximately  10-dB  difference  at 
1  GHz.  This  is  due  to  the  fact  that  the  phase  decreases  as 
the  frequency  decreases  and  all  the  branching  elements  star 
to  act  coherently. 

Results  of  Figs.  2  and  3  indicate  that  it  is  important  to 
incorporate  coherent  phase  interaction  in  the  calculation  of 
scattering  from  a  collection  of  scatterers  exhibiting  clustering 
structure,  which  is  a  typical  characteristic  of  vegetation.  As 
found  in  this  section,  the  architecture  of  a  cluster  clearly 
displays  its  own  signature  in  the  frequency,  angle,  and  po¬ 
larization  responses.  In  the  following,  we  will  use  soybean 
as  an  example  to  illustrate  the  procedure  of  implementing  a 
branching  model  for  vegetation  and  how  the  phase  interaction 
terms  can  be  properly  incorporated. 


111.  Branching  Model  for  Soybean 

Soybean.  Glycine,  nowadays  constitutes  a  very  most  im¬ 
portant  food  legume  in  China.  Korea.  Japan,  aiKl  Malaysia. 
Large  quantities  of  soybean  are  also  grown  in  the  U.S.  [13].  A 
soybean  is  an  erect  branching  plant  having  a  stem  and  terminal 
branches  attached  with  leaves.  The  leaf  of  soybean  is  so  called 
a  compound  leaf  which  is  comprised  of  three  leaflets  (most 
species)  and  a  petiole,  which  is  a  slender  stem  supporting 
the  leaflets.  A  few  species  have  five  or  seven  leaflets  on  one 
compound  leaf.  The  leaves  usually  fall  before  the  fruits  (pods) 
are  mature  and  the  pods  in  general  contain  two  or  three  seeds. 
An  excellent  color  illustration  showing  the  architecture  and 
length  scale  of  soybean  can  be  found  in  [13]. 

In  this  section,  soybean  will  be  modeled  by  a  two-scale 
branching  model  shown  in  Fig.  4.  A  leaf  will  be  modeled  by 
a  cylinder  (petiole)  associated  with  a  dielectric  disk  (leaflets). 
For  convenience,  we  will  use  "branch”  to  represent  the  cylin¬ 
der  so  that  the  effects  of  terminal  branches  (ignored  in  this 
simplified  soybean  picture.  Fig.  4)  and  petioles  are  both  ac¬ 
counted  for  by  this  modeling  element  “branch."  Leaf  elements 
and  fruits  are  then  connected  to  a  cylinder  representing  the 
stem  to  form  the  next  order  cluster — a  plant.  Eventually, 
polarimetric  backscattering  coefficients  for  soybean  fields  (a 
colleaion  of  soybean  plants)  will  be  derived  for  correlated 
soybean  locations  described  by  a  pair  distribution  function.  In 
the  following,  scattering  from  these  two  scales  regarding  the 
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Fig.  4  Configuralion  of  the  two-scale  branching  model  for  soybean. 


The  phase  factor  takes  into  account  the  fact  that  the  circular 
disk  (leaflets)  is  attached  to  the  end  of  the  branch.  A  translation 
in  the  position  results  in  a  phase  term  for  the  field  received 
by  the  antenna. 

Note  that  Somp  (Sjmp)  is  not  the  scattering  matrix  elements 
of  the  branch  (disk).  In  fact,  5omp  and  Sjmj,  consist  of  a  term 
corresponding  to  the  scattering  matrix  element  of  the  object, 
and  other  terms  accounting  for  the  effects  of  propagation 
through  the  other  vegetation  material  and  reflection  of  incident 
and  scattered  fields  by  ground  surfaces.  Explicit  forms  for 
5omp  or  Sjmp  are  given  by  (B20)-(B23),  Appendix  B. 

The  phase  angle  6m  can  be  calculated  by  using  (2)  for  four 
scattering  mechanisms. 


Fig.  5.  Four  main  mechanisms  for  the  scattering  from  an  object  above  a 
surface. 

inner  structure  of  a  soybean  plant  and  the  scale  related  to  the 
distribution  of  soybean  plants  will  be  described  separately. 

A.  Leaf 

Scattering  from  a  leaf  element  in  the  presence  of  a  flat 
ground  surface  will  be  considered  in  this  section.  Due  to  the 
presence  of  ground  surface,  the  incident  wave  may  be  reflected 
from  the  scatterer  via  four  main  mechanisms  shown  in  Fig.  S. 
Hence,  after  neglecting  the  reflection  from  the  upper  boundary 
and  the  multiple  scattering  between  scatterers  and  ground,  the 
backscattering  matrix  from  a  leaf  element  consists  of  four  main 
terms,  direct  (1),  scatterer-ground  (2),  ground-scatterer  (3),  and 
ground-scatterer-ground  (4),  fig.  5. 

Neglecting  the  multiple  scattering  between  the  branch  and 
the  circular  disk,  the  total  backscattering  matrix  5  is 

Sp  =  Sop  +  Sjp  (6) 

where  So  is  the  return  from  the  branch  and  S,  is  the  return 
from  the  disk.  Subscript  “p”  denotes  one  of  the  four  antenna 
polarization  states,  HH,  HV,  VH,  and  W,  where  the  first  and 
the  second  letters  denote  antenna  transmitting  and  receiving 
polarizations,  respectively.  “//”  represents  horizontal  polariza¬ 
tion  and  “V”  vertical  polarization  (Appendix  A).  Furthermore, 
each  term  is  constituted  by  four  terms  corresponding  to  four 
scattering  mechanisms  (Appendix  B), 

4 

•^Op  —  ^  "  Somp 
m=l 

4 

Sjp  =  ^  "  Sjmp  f'xp(idm).  (7) 


=  2{— kzi  cos 'r sin  0  +  kyi  sin 7 sin/? -I-  k,iCOs/?)L  (8a) 
62  =  2(— kx,  cos sin  0  + kyi  sin-f  sin 0)L  (8b) 

6s  =  62  (8c) 

64  =  2{-kxi  cos  "y  sin  0  +  kyi  sin  ■ysin0  -  k,i  cos  0)L  (8d) 


where  a,  0,  and  7  are  the  Eulerian  angles  describing  the 
orienution  of  the  branch  ([1],  p.  158),  and  L  is  the  length 
of  branch. 

Note  that  6m  is  a  function  of  the  length  and  orientation 
angles  of  the  branch.  Assuming  the  size  and  orientation  of 
branches  are  independent  of  the  scattering  properties  of  leaves, 
ensemble  average  of  Sp  yields 

4  4 

Sl,.f,p  =  (Sp)  =  Y,  (Somp)  +  Y,  (5;mp>(exp(iem))  (9) 

msl  msl 


and  the  correlation  between  “p"  and  “q"  polarizations 


Clear, pg  —  (■^p^g) 


=  (^op55,)-F(Sj,s;,)  +  (SopS;,) 

+  (5;p5o%) 


(10) 


where 

(SopSo,)  =  Y,  ^  {^OmpSoj,g) 

T71  =  l  nsl 

{Sjp^jij)  —  ^  Y>  (SjmpSjny) 
m=l  nsl 

■  (exp(ttfm)exp(-t^„)> 

4  4 

(5op5j*g)  =  51  (^Ompexp(-t«„))(Sj*„,) 

m=l  n  =  l 
4  4 

(■SjpSog)  =  5Z  51  (5jmp)(5^„,exp(l0m)).  (11) 

msl nsl 

Each  double  sum  can  be  decomposed  into  a  sum  of  four 
self  terms  (nt  =  n)  and  12  interaction  terms  between  four 
scattering  mechanisms  (m  ^  n).  Note  that  the  reciprocity 
between  scattering  mechanisms  2  and  3  demands  that 

5o2a6  =  5o36a,  (1-3) 

Sj2ab  =  Sj3hn',  n.b=-h.v.  (l?bl 
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Thus,  the  interaction  terms  between  mechanisms  2  and  3  are 
essentially  equal  to  the  self-terms  of  mechanisms  2  or  3.  In 
other  words,  the  contribution  from  mechanisms  2  and  3  should 
be  added  coherently.  However,  these  interaction  terms  have 
been  ignored  in  the  radiative  transfer  theory. 

In  the  following,  the  mean  and  correlation  of  scattered  fields 
given  in  (9)  and  (10)  of  a  leaf  element  will  be  combined  with 
the  scattering  characteristics  of  the  center  stem  to  calculate  the 
polarimetric  scattering  from  a  soybean  plant. 


where  N,  signifies  the  total  number  of  soybean  plants  within 
the  illuminated  region.  represents  the  iocatioa  of  soybean 
plant  “m”. 


Pm  =  ymV 

(17) 

k.0  —  "F  kytV- 

(18) 

B.  Soybean  Plant 

Again  by  invoking  the  single  scattering  approximation,  the 
scattering  matrix  with  polarization  “p”  from  a  soybean  plant 
can  be  written  as  a  sum  of  the  returns  from  the  stem  (Sq),  leaf 
elements  (5^),  and  fruits  (S[ ).  For  simplicity  Sop  is  assumed 
to  be  independent  of  and  1  arc  statistically 

identical  and  mutually  independen; 

.V,  -'■/ 

=  +  +  (13) 

j«i  j=i 

where  .V(  represents  the  total  number  of  leaf  elements  Nf 
the  number  of  fruits.  Note  that  propagation  factors  have  been 
combined  with  the  scattering  matrix  of  the  object  to  define 
Sop  and  Sjp.  Appendix  B  (B20)-(B2j). 

The  mean  scattered  field  from  a  soybean  plant  is  given  by 

5„vb.an.p  =  (Sp)  =  (Sop)  +  .V,5|„f,p  +  yf(S[p}.  (14) 

Taking  the  ensemble  average  of  the  second  moments  of 
scattered  field  Sp  gives  the  correlation  between  “p”  and  “g” 
polarizations,  as 

C’»oybeaii.pq_=  ('^Op^g^)  +  ^lO\ftf,pq  +  f(S^jpS^^  ) 

+  (5op)(.v,s,;.f.,  +  .\v(sf;)) 

+  +  SU,{S[p}) 

+  AV(,\v-i)(sf,)(sf;) 

+  .v,(,v,-i)s,„,.pS,;.,,,  (15) 


Assuming  that  all  the  soybean  plants  are  mutually  indepen¬ 
dent  and  statistically  identical,  we  obtain  the  following  field 
correlation. 


(■^p^^)  —  y$Ctoyhrhn,pq  "I"  ^soybean. (19) 


where 

.'"s 

(exp[2tfcp  •  (p^  -  P„)]  )•  (20) 

m  =  l  n#m 

The  assumption  of  statistically  identical  soybean  plants  leads 
to 

F  =  N.(N.  -l)j dp„  I dp„P2(p„.p„) 

A  A 

•  exp[2tfcp  •  (p„  -  p„)]  (21) 

where  .4  is  the  illuminated  area  and  P2  is  the  probability  dis¬ 
tribution  of  the  locations  of  two  soybean  plants.  According  to 
Bayes  rule,  P2(p„  .  pj  can  be  expressed  as  P(p„lp„)P(Pr,). 
Assume  a  uniform  distribution  for  p„,  thus.  P(p„)  =  1/-4. 
Subsequently,  a  normalized  two-point  distribution  p;  can  be 
defined  by  g2{p„.pn)  =  AP(p„\p„).  For  symmetric  and 
space-invariant  relative  position  distribution,  p2(p„.p„)  = 
g(Pjj,  -  p„)  is  called  the  pair  distribution  function  [1].  Thus, 

F  =  ~  j  dpgip)  exp(2tfcp  •  p) .  (22) 

A 


Further  assuming  TV,  :?>  1  and  substituting  F  into  (19)  yield 
(^SpSg'*^  —  rtjQSsoybean.p^goybean.?  )  (23) 

where 


where  all  the  leaf  elements  and  fruits  have  been  assumed  to 
be  mutually  independent  and  statistically  identical.  Note  that 
in  the  above  expression  (5^p)  and  (Sj^Sj")  corresponding 
to  the  contribution  of  leaf  elements  have  been  replaced  by 
^i*af.p  znd  respectively'. 

C.  Scattering  from  Soybean  Fields 

In  this  section,  polarimetric  backscattering  coefficients  from 
a  collection  of  soybean  plants  within  the  radar  beam  footprint 
will  be  derived  with  the  correlated  locations  of  soybean  plants 
described  by  a  pair  distribution  function  [1].  Applying  single 
scattering  assumption  gives  the  total  backscattered  field  Sp  as 

•Vy 

Sp  =  ^  S„p  exp(2tfcp  •  pm)  (16) 

m  =  l 


A 


(24) 


and  n,  =  N,/A  is  the  density  of  stems  (number  per  unit  area) 
Several  pair  distribution  functions  have  been  used  to  eval¬ 
uate  the  integral  for  Q  [1).  For  independent  scatterer  distri¬ 
bution,  g(p)  equals  one  everywhere.  For  a  hole  correction 
approximation,  the  pair-distribution  function  has  the  form 


9{p)  = 


p  <  R 
p>  R. 


(25) 


That  means  each  soybean  plant  occupies  a  certain  area  with  a 
radius  R  where  the  stems  of  other  soybean  plants  are  not  al¬ 
lowed  to  grow.  In  other  words,  the  minimum  spacing  between 
two  stems  is  R.  Outside  of  that  region,  stems  are  uniformly 
distributed  like  the  independent  scatterer  assumption. 
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Let  the  illuminated  area.  A,  approach  infinity.  The  integral 
for  Q  is  given  by 


Q  =  4n^6{kp) 


d p  exp(2tfcp  • 


P)> 


for  independent  scatterer 
for  hole  conection 


(26) 


where  the  integral  car  be  carried  out  and  given  as 
j  dp  exp(2iJp  •  p)  =  kR^ 

t<R  " 


To  obtain  the  expression  of  R  for  the  hole  correction,  we 
assume  that  N,nR^  =  pA  which  relates  the  area  of  soybean 
canopy  to  the  total  illuminated  area,  where  /r  is  a  constant 
determining  the  minimum  spacing  between  stems.  Thus, 


R  = 


(28) 


The  cor£elation  of  measured  polarizations  for  oblique  inci* 
dence  (fcp  #  0)  reduces  to 


=  A’, 


L  soybean. p9  ^  ^soybean.p'^sovbean.y 


where 


j  =  2koRsin6,. 


(29) 

(30) 


To  ensure  the  power  of  each  polarization  (5p5*)  is  always 
greater  than  zero,  parameter  p  must  be  smaller  than  one.  Also 
note  that  when  p  =  0,  the  hole  correction  reduces  to  the 
independent  scatterer  assumption. 

The  polarimetric  covariance  matrix  [2]  having  the  unit  as 
the  backscattering  cross  section  per  unit  area  is  defined  by 

cTpj  =  —  (SpS, ).  (31) 


Substituting  (29)  into  (31)  and  summing  up  the  contribution 
from  a  rough  vegetation  floor  give  the  total  polarimetric 
backscattering  coefficients  as 


<Tp,  =  47rn, 


^^•oybean.pfl  2p  ^soyb«an,p^«oyb«an,9 


Ofpq(7ground.p7 


(32) 


where  cTground.pf  is  the  contribution  from  a  rough  terrain 
surface.  Wrchhoff  approximation  will  be  used  in  our  soybean 
model  to  compute  (Tground.p^  (1]-  Furthermore,  due  to  the 
attenuation  introduced  by  the  vegetation  cover,  <Tgjound.p« 
is  multiplied  by  a  factor  Op,  which  is  calculated  in  the 
following  manner.  The  propagation  factor  due  to  the  vegetation 
cover  appearing  in  (B20)  for  each  polarization  is  evaluated  at 
the  bottom  of  the  vegetation  layer  (z,  =  -dz).  Multiplying 
the  propagation  factor  of  "p"  polarization  with  the  complex 
conjugate  of  that  of  “(/”  polarization  yields  the  factor  r»p,. 


D.  Models  of  Vegetation  Elements 

In  the  implementation  of  the  above  branching  model  for  soy¬ 
bean,  dielectric  cylinders  are  used  to  represent  stems,  branches, 
and  fruits,  whereas  leaflets  are  represented  by  dielectric  disks. 
Scattering  matrices  are  obtained  by  using  the  finite  cylinder 
approximation  for  cylinder  [3],  [9]  and  the  physical  optics 
approximation  for  disk  [6],  [7].  Note  that  in  [3],  (6j,  (7], 
and  [9]  the  exp{iuit)  convention  was  used  and  the  vector 
definitions  of  horizontal  and  vertical  polarizations  are  also 
different  from  those  defined  by  (A3a)  and  (A3b).  Therefore, 
scattering  matrices  are  rederived  using  the  exp(-twf)  notation 
and  definitions  of  polarization  vectors  (Appendix  A)  for  these 
two  elements  [14].  These  formulas  were  shown  to  be  applica¬ 
ble  over  a  wide  range  of  frequencies.  However,  it  is  known 
that  both  approximations  give  nonreciprocal  bistatic  scattering 
matrices.  To  ensure  the  reciprocity,  the  scanering  matrices  in 
reciprocal  directions  are  averaged.  For  oriented  cylinders  or 
disks,  three  Euler  angles  are  used  to  specify  the  orientation 
of  these  scatterers.  The  scattering  matrix  of  each  element  is 
calculated  in  its  local  coordinates  and  then  transformed  back 
to  the  antenna  coordinates. 

E.  Treatment  of  Coherent  Wave  Propagation 

The  propagation  is  accounted  for  by  using  the  Foldy’s 
approximation,  Appendix  C.  Note  that  the  forward  scattering 
matrix  of  a  two-scale  or  multiscale  system  is  simply  the  sum 
of  forward  scattering  matrices  from  every  components  if  the 
multiple  scattering  between  scattering  elements  is  ignored. 
Note  that  the  averaging  processes  for  the  propagation  and 
scattering  of  the  waves  are  carried  out  independently  in  this 
paper,  which  is  an  approximation. 

IV.  Comparison  of  Theoretical  Results  with 
Experimental  Data 

A.  Ground  Truth  of  Soybea>!  Data 

The  soybean  data  analyzed  were  backscattering  coeflicients 
of  HH,  HV,  and  W  measured  at  5.3  GHz  by  the  Centre 
d’Etude  Spatiale  des  Rayonnements  in  France  over  a  2  1/2- 
month  period  in  1986.  The  height  of  radar  was  fixed  at  14.5  m. 
and  the  3-dB  beam  width  of  radar  is  9.2°.  Each  of  the  data  set 
was  measured  at  three  polarizations  for  six  incident  angles. 
Along  with  backscattering  coefficients,  extensive  ground  truth, 
including  vegetation  fractional  volume  (/),  water  content, 
canopy  height  (h),  ground  surface  roughness  (correlation 
length  I  and  standard  deviation  o),  and  soil  moisture  [16], 
is  summarized  in  Table  I  with  the  calculated  permittivity  for 
vegetation  and  soil.  The  soil  at  the  experimental  site  contains 
11.1%  sand,  61.7%  loam,  and  27.2%  clay  [16].  The  first 
column  (date)  corresponds  to  the  time  when  the  data  were 
collected.  For  example,  0717  represents  July  17,  1986.  The 
permittivity  of  vegetation  material  was  deduced  [16]  using  the 
formula  given  in  [17]  based  on  the  measured  volumetric  water 
content.  The  permittivity  of  vegetation  material  showing  a 
decreasing  trend  is  due  to  the  faa  that  the  measured  volumetric 
water  content  of  soyhc.in  decreased  during  the  experimental 
period  [16].  The  periiiiuivity  of  ground  soil  was  e.ilciilatcd 
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TABUl 

Measured  and  Calculaied  Ground  Truih  Parameters 


h 

/ 

Vegelitioii 

a 

1 

Soil 

Plam  Density 

Date 

(m) 

(*) 

f 

(cm) 

(cm) 

n.  (No./m^) 

0717  0.20  0.11  07.84  +  / 20.75  1.7  14  13.09 +  / 2.03  18 

0725  0.29  0.18  60.07  +  1  18.44  1.7  14  9.33 +  »  1.04  18 

0807  0.44  0.29  48.93  + 1  15.13  1.6  14  11.93  +  /  2.34  18 

0814  0i2  0.60  45.44  +  /  14.10  1.6  14  13.57  +  /  3.34  18 

0820  0.59  0.87  40.22  + 1  12.55  1.5  13  7.02  +  1  1.04  18 

0828  0.68  0.89  32.43  +  /  10.23  1.4  13  10.30  + 1  3.00  18 

0903  0.745  0.88  28.77  +  / 9.14  1.3  12  8.70  +  /  1.49  18 

0912  0.825  0.68  24.62  +  /  7.90  1.2  12  5.59  +  /  0.68  18 


TABLE  II 

Esttmatid  Fractionai.  Volumes  of  Vegetation  Components 


Leaves 

Branches  and 

Fruits 

Date 

(%) 

Stems  (%) 

(%) 

0717 

0.055 

0.055 

0 

0725 

0.09 

0.09 

0 

0807 

0.135 

0.149 

0.006 

0814 

0.247 

0.310 

0.043 

0820 

0.32 

0.45 

0.1 

0828 

0.293 

0.36 

0.237 

0903 

0.264 

0.29 

0.326 

0912 

0.175 

0.196 

0.309 

in  accordance  with  the  soil 

texture  and  moisture 

using  the 

formula  given  in 

[18]  and  has  been  reported  in 

[16].  Note 

that  the  vegetation  fractional  volume  increases  at  the  early 
stage  and  stays  constant  at  the  fully  grown  stage.  Afterwards, 
a  decrease  of  the  fractional  volume  is  due  to  the  fall  of  leaves 
and  steady  dry  ing  of  vegetation  material. 

In  addition,  the  weight  of  vegetation  components  was  mea¬ 
sured  at  5  different  days  (0804,  0820,  0902,  0915,  and  1003) 
and  the  fractional  weight  was  plotted  versus  the  day  in 
Fig.  6  where  day  1  represents  July  17  and  day  80  represents 
October  4.  Note  that  the  weights  of  branches  and  stems  were 
not  measured  separately.  The  fruits  appeared  at  August  7  and 
their  weight  became  dominant  at  the  end  of  the  season  reaching 
approximately  50%  of  the  total  vegetation  weight.  Due  to 
the  fact  that  weight  measurements  in  general  were  carried 
out  at  different  days  from  the  time  when  the  backscatter 
was  measured,  the  fractional  weight  of  each  component  at 
the  days  (Table  I,  Column  1)  when  the  radar  measurements 
were  carried  out  was  estimated  by  linearly  interpolating  the 
measurements.  Further  assuming  the  densities  (weight  per 
unit  volume)  of  all  the  vegetation  elements  are  the  same,  the 
fractional  volumes  of  vegetation  components  were  estimated 
as  the  product  of  their  fractional  weight  and  the  total  vegetation 
fractional  volume  (Table  I,  Column  3)  and  listed  in  Table  II. 

Some  other  parameters  necessary  for  theoretical  calculation 
are  selected  in  the  following  manner.  The  length  of  stem  is  set 
to  be  the  measured  canopy  height  and  one  vegetation  layer  is 
assumed  so  that  d\  =  d2-  Other  parameters  shown  in  Table  III 
regarding  the  size  and  density  of  vegetation  components  are 
adjusted  so  that  their  fractional  volumes  match  those  shown 
in  Table  II.  Hence,  there  are  only  six  independent  parameters 
(one  for  disks,  three  for  branches  and  stem,  and  two  for  fruits) 


Fig.  6.  Fractional  weight  of  soybean  components  at  5  different  davs  (0804. 
0820.  0902.  0915.  and  1003). 


for  each  date  in  Table  III.  which  are  selected  in  such  a  way  that 
a  good  agreement  between  theoretical  results  and  experimental 
data  is  reached.  Note  that  Ni{Nf)  is  calculated  by  multiplying 
the  density  of  branches  (fruits)  with  the  height  of  canopy. 

In  selecting  the  parameters  listed  in  Table  III,  the  sizes  of 
vegetation  elements  are  ensured  to  agree  with  botany  studies. 
For  example,  the  thickness  of  disk  (~0.2-0.3  mm)  agrees  well 
with  the  usual  thickness  of  soybean  leaflets.  The  radius  of  disks 
is  chosen  so  that  the  area  of  disk  corresponds  to  the  total  area 
of  three  leaflets.  It  can  be  seen  that  the  radii  of  disk,  stem, 
branch,  and  fruit  are  in  good  agreement  with  the  sizes  shown 
in  the  illustration  of  soybean  elements  [13].  Moreover,  the 
radius  of  the  disk  agrees  very  well  with  the  measured  average 
area  of  a  leaf  at  4  days  as  shown  in  Fig.  7  where  the  radius 
for  a  circular  disk  with  the  same  area  is  also  indicated. 

In  addition,  the  probability  distributions  of  orientation  an¬ 
gles  of  disk,  branch,  and  fruit  are  chosen  to  be 


p(g\  =  /  sf  if  0  5  ^  £  2  • 

\  0,  otherwise. 

p/ o  \  /  ®i^^di»k.  if  “o'  — 

P{(3d..k)-|o,  otherwise. 

f  |siD(2/dfruit)|.  if  y  $  /Jfruil  S 

[  0.  otherwise. 

Ph)  =  Phf,..,)  =  ^ 


sin/3di.k. 


lsin(2,9frui,)|, 


(33) 

(34) 

(35) 

(36) 
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TABLE  m 

Selected  Parameters  kr  Theoretical  CaijCULAik>n 


Dale 

Disk 

Radius 

(cm) 

Disk 

Thickness 

(mm) 

Branch 

Radius 

(cm) 

Branch 

Length 

(cm) 

Density  of 
Branches 
(No./m) 

Stem 

Radius 

(cm) 

Fruit 

Radius 

(cm) 

Fruit 

Length 

(cm) 

Density 
of  Fruits 
(No./m) 

0717 

40 

0.22 

0.12 

10 

27.5 

0.24 

0 

07^^ 

4.3 

0.24 

0.16 

9 

36 

0.275 

— 

— 

0 

4.8 

0.25 

0.2 

11 

41 

0.3 

0.25 

2.0 

8 

0814 

5.0 

0.27 

0.23 

11 

63 

0.44 

03 

3.0 

28 

5.1 

0.32 

0.25 

11 

69 

0.55 

0.4 

3.0 

37 

0828 

5.4 

0.31 

0.25 

11 

57 

0.5 

0.5 

3.6 

47 

5.2 

0.31 

0.24 

11 

54 

0.44 

0.5 

4.0 

57 

0912 

4.7 

0.31 

0.2 

11 

45 

0.4 

0.49 

4.0 

57 

Vorioiion  of  the  Average  Area  of  A  Leaf 


Fig.  7.  Measured  aserage  area  of  a  leaf  at  4  days  same  as  the  first  4  days 
for  the  data  shoun  in  Fig.  6.  The  radius  of  a  circular  disk  of  the  same  area 
is  indicated 

SO  that  branches  are  more  horizontally  distributed  and  disks 
(leaflets)  tend  to  bend  toward  the  ground  to  account  for  the 
effects  of  gravitational  force.  As  to  the  fruits,  most  of  them 
are  inclined  45°  downward  in  accordance  with  the  illustration 
shown  in  [16].  Uniform  distributions  for  7  and  7fnii,  are 
selected  in  order  to  have  an  azimuthally  symmetric  soybean 
canopy.  There  was  no  detailed  ground  truth  measurement  on 
the  structure  of  the  soybean  at  the  time  when  the  experiment 
was  carried  out.  Hence,  the  above  probability  density  functions 
were  selected  so  that  a  good  agreement  between  theory  and 
data  is  achieved.  To  avoid  introducing  additional  parameters, 
the  same  orientation  distributions  are  used  for  all  data  sets. 

B.  Comparison  of  Theoretical  Results  with  Experimental  Data 

Figs.  8(a)-(h)  compare  experimental  data  with  theoretical 
results  which  are  calculated  using  (33)  for  /l  =  1  (hole 
correction).  There  is  a  good  agreement  between  theoretical 
results  and  experimental  data  except  for  the  data  sets  where 
the  fractional  volumes  of  vegetation  are  small  (early  stages 
of  soybean).  For  the  first  data  set,  theoretical  contribution 
from  the  ground  surface  alone  is  3  dB  near  normal  incidence, 
which  is  larger  than  the  measured  backscattering  coefficients. 
That  may  be  due  to  either  the  surface  scattering  model 
(Kirchhoff  approximation)  is  not  appropriate  for  this  case  or 
the  nicnsiiiciJ  roughness  parameters  of  ground  are  not  accurate 


enough.  According  to  a  simulation  using  [1,  eq.  (59b)]  based 
on  Kirchhoff  approximation,  it  is  found  that  a  10%  error  in  <7 
and  I  can  result  in  a  change  of  1  dB  for  surface  scattering. 

Using  the  vegetation  parameters  for  date  0814,  Fig.  9  il¬ 
lustrates  the  deficiency  of  independent  scatterer  assumption 
(ft  =  1)  which  is  not  a  good  description  for  soybean  plant 
distribution.  As  observed  from  vegetation  like  trees  and  crops, 
there  is  always  a  certain  spacing  between  plants.  In  other 
words,  individual  plant  tends  to  exclude  other  plants  from  its 
territory.  Hence,  hole  conection  is  more  appropriate  for  this 
situation.  Note  that  the  peaks  observed  in  the  HH  and  W 
backscatter  for  independent  scatterer  assumption  near  normal 
incidence  are  greatly  reduced  by  using  the  hole  correction 
factor  p  -  1.  and  the  theoretical  results  at  higher  angles 
of  incidence  are  essentially  unaffected  where  both  2J\{t)Ii 
and  SMvb»»n.p  are  observed  to  be  small.  The  reduaion  of  the 
sharp  peaks  near  normal  incidence  is  due  to  the  destruaive 
interference  introduced  by  the  hold  correaion  approximation 
which  prevents  two  soybean  plants  from  overlaying  each 
other.  It  should  also  be  noticed  that  the  cross  polarized  return 
HV  is  not  affected  by  the  hole  correction.  This  is  due  to 
the  fact  that  the  coherent  backscatter  5K>yb»An.hv’  is  zero  for 
azimuthally  symmetric  canopy.  In  order  to  fully  demonstrate 
the  effects  of  coherence,  experimental  data  collected  at  many 
frequencies  in  conjunction  with  sufficient  ground  truth  data 
must  be  used.  The  soybean  data  presented  in  this  paper  has 
extensive  ground  truth,  but  the  measurement  was  limited  to  C- 
band.  For  soybean,  because  its  internal  scales  are  in  general 
of  the  order  or  greater  than  the  wavelength  at  5.3  GHz.  we  do 
not  expect  the  coherence  effect  to  be  significant  at  C-band.  We 
expect  that  coherence  will  be  significant  at  lower  frequencies 
such  as  L-  and  P-bands. 

V.  Discussion 

In  this  paper,  a  two-scale  branching  model  has  been  imple¬ 
mented  for  soybean  and  validated  with  experimental  data.  To 
generalize  the  above  model  to  vegetation  with  more  branching 
scales,  such  as  trees,  the  following  procedure  can  be  applied. 
The  model  for  trees  will  consist  of  several  scales  including 
trunk,  bough,  branch,  twig,  and  leaves.  In  the  branching 
process,  the  leaves  are  first  attached  to  the  smallest  branch 
(twig).  The  resulting  clustered  twig-leaves  pattern  is  then 
attached  to  a  larger  branch  lo  form  the  next  order  clustering. 
Repeating  this  process  with  branches  attached  to  bough  .ind 


»terii>9  Coafflclonte  (dB)  Bacl<«catt«rWig  Co«fflclenU  (dB)  Bocli«c<rtf«1ng  CiMtflclw)ts  (dB) 


398 


IEEE  TKANSACnONS  ON  OEOSODICE  AND  KEMOTE  SENSINC.  VOL.  3&  MX  Z  MARCH  1992 


(a)  SOYBEAN  8S0717 


0  20  40  80 

Inddant  Angla  (ddgnn) 


(b)  SOYBEAN  aeom 


0  20  40 


Incidant  Angla  (dagfaa) 


(e)  SOYBEAN  860807 


0  20  40  60 

Ineidant  Angla  (dagraa) 


(d)  SOYBEAN  860814 


0  20  40 


Ineidant  Angla  (dagraa) 


(a)  SOYBEAN  860820 


(f)  SOYBEAN  880828 


Ineidant  Angla  (dagraa) 


Ineidant  Angla  (dagraa) 


(g)  SOYBEAN  860003 


(h)  SOYBEAN  860012 


0  20  40  60 

Ineidant  Angla  (dagraa) 


Fig.  8.  Comparison  of  theoretical  results  and  experimental  data.  Hole  correction  with  /t  =  1  is  used  for  the  pair-distribution  function  of  soybean  plants. 
H,  1',  and  .V  represent  measured  backscallering  coefficients  of  HH.  VV,  and  HV  polarirations,  respectively. 
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Effects  of  It 


0  20  40  60 

Incidsnt  Angle  (degree) 


Fig.  9.  Comparison  of  theoretical  results  generated  by  independent  scatterer 
approximation  =  1)  and  hole  correction  p  =  1.  Parameters  for  date  0814 
are  used  for  calculation. 

boughs  attached  to  trunk  finally  forms  a  tree.  By  changing  the 
sizes  and  growth  direction  of  each  scale,  various  tree  patterns 
can  be  generated. 

A  recursive  procedure  can  be  used  to  calculate  the  backseat- 
tering  parameters  for  trees  by  taking  advantage  of  their  branch¬ 
ing  architecture.  For  a  tree  with  the  above  multiscaie  pattern, 
the  problem  of  scattering  from  a  multiscale  structure  can 
be  decoupled  into  a  cascade  of  scattering  from  two-scale 
structures  (subscale  and  primary  scale).  For  example,  in  the 
case  of  the  smallest  clustering,  subscale  corresponds  to  leaves 
and  primary  scale  represents  a  twig.  The  averages  polarimet- 
ric  scattering  matrix  and  covariance  matrix  of  a  two-scale 
structure  can  be  expressed  in  terms  of  those  of  subscale  and 
primary  scale.  The  resulting  averaged  polarimetric  scattering 
matrix  and  covariance  matrix  will  conespond  to  those  of  the 
subscale  of  the  next  order  two-scale  structures.  Repeating  the 
above  recursive  procedure  eventually  leads  to  the  averaged 
scattering  matrix  and  covariance  matrix  of  a  tree.  For  example, 
a  branching  model  for  soybean  is  implemented  by  following 
the  above  procedure. 

The  advantages  of  this  modeling  approach  over  the  conven¬ 
tional  methods  are  that  the  multiscale  structure  of  vegetation 
is  considered  and  the  coherence  or  phase  interference  between 
each  component  is  taken  into  account.  In  the  past,  theoretical 
models  did  not  consider  the  structure  of  vegetation.  The  phase 
interaction  between  vegetation  components  was  ignored  and  an 
incoherent  sum  was  used  to  compute  the  total  scattering  power 
from  every  components.  In  the  multiscale  branching  model, 
the  phase  interaction  or  coherence,  which  can  be  observed 
at  the  inner  scales  of  a  single  plant  and  the  spatial  scale 
of  plant  distribution,  is  retained.  Hence,  the  model  should 
be  applicable  over  a  larger  frequency  spectrum  and  range  of 
incident  angles  than  the  conventional  approaches  for  remote 
sensing  of  vegetation. 

appendix  a 

Notation  and  DEFiNmoN  of  Vectors 

The  notation  and  definition  of  wave  vectors  and  polarization 
vectors  are  given  in  this  appendix  (Fig.  10).  The  wave  vector 


Fig.  10.  Illustration  of  vector  definition  of  wave  vector  F  and  polarization 
vecton. 


Rh.Rv 

Fig  1 1 .  Configuration  of  a  scanering  object  located  inside  a  layered  medium. 


of  a  wave  propagating  in  the  direction  {9, 4>),  is  defined  as 


where 


k  =  xkg  +  yky  +  ikt 

(Al) 

kz  =  Itocos^sin^ 

(A2a) 

ky  =  A^sin^sintf 

(A2b) 

kz  =  k{)Cos9 

(A2c) 

where  ifco  is  the  free  space  wavenumber.  The  horizontal  and 
vertical  polarization  vectors  of  electric  fields  are  defined  as 


t)(*a)  = 


hxk 
h  X  k 


iky  -  yfci 

kp 

-xkxkf  -  ykyk,  -I-  ikj 

kckp 


(A3a) 

(A3b) 


The  above  definition  is  used  for  both  incident  and  scat¬ 
tered  fields.  For  an  incident  wave  with  direction  of  incidence 
{9l,|^>,),  its  wave  vector  is  denoted  as 


where 


i  =  ikzi  +  ykyi  -1-  ik,i 

(A4) 

kzi  =  kocoGipi  sin  6, 

(A5a) 

kyi  =  A:o  sin  ^  sin  9i 

(A5b) 

k,i  =  A-Qcostfi. 

(A5c) 
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propagation  effects,  and  ground  reflectkm,  are  the  resulting 
scattering  matrices  as  would  be  observed  by  the  receiver  for 
four  main  scattering  mechanisms. 

APPENDIX  C 

Propagation  of  Coherent  Fields 
Under  Foldy  s  approximahon 

For  sparse  media  like  vegetation,  Foldy's  approximation 
[1],  [15],  [19],  [20]  can  be  employed  to  account  for  the 
absorption  and  scattering  effects  introduced  by  the  medium. 
The  propagation  of  coherent  wave  vuh  £/,  and  £,.  as  the 
horizontal  and  vertical  components  of  the  electric  held  is 
governed  by  the  following  (1,  p.  139): 

=  (iko  +  Mhh)Eh  +  MhvEv  (Cl) 

^  =  hU  Eh +  (iko  +  A/,.,  )£,  (C2) 

where  s  is  the  distance  along  the  direction  of  propagation,  and 

A/,,=  jJ  =  v.h.  (C3) 

Here  {S.  o)  is  the  direction  of  propagation,  and  no  is  the  density 
of  scatterers  (number  per  unit  volume).  In  the  above  equation, 
the  average  is  carried  out  over  the  size  and  orientation  of  em¬ 
bedded  vegetation  elements.  For  the  soybean  model  presented 
in  Section  III. 

+  -A'/(/j/,br«nch(^-  <t>:  0.  <t>)} 

+  (C4) 

u  here  h  is  the  height  of  vegetation  canopy,  n,  is  the  number 
ol  stems  per  unit  area,  and  A'l  and  Nf  are.  respectively,  the 
.rage  number  of  branches  and  fruits  per  stem, 
if  the  canopy  structure  exhibits  statistically  azimuthal  sym¬ 
metry.  there  will  be  no  coupling  between  horizontal  and 
vertical  components  of  the  coherent  field.  It  can  be  shown 
that 


-V/bt  =  ^fvh  =  0.  (C5) 

Hence,  the  effective  propagation  constant  of  horizontally  and 
vertically  polarized  coherent  waves  are  given  by 


Rh  —  fco  +  — : — 

1 

(C6) 

J.  J.  .. 

f 

(C7) 

References 

|1|  L  Tsang.  J  A  Kong,  and  R  T.  Shin.  Theory  of  Microwave  Remote 
Sensin/i  New  York:  John  Wiley  &  Sons.  1985. 

(2|  M  Borgi  jud.  S  V  Nghiem.  R  t  Shin,  and  J.  A.  Kong.  “Theoreiicai 
models  ;  polanmeiric  microwave  remote  sensing  of  earth  terrain." 
J.  Electromagnetic  Waves  and  Appl.,  vol.  3.  no.  1,  pp.  61-81,  1989. 

[.Ij  .M  A  K.iram.  A  K  Fiinc.  and  Y  M  M  Anlar,  "Fleciromagnelic  wave 
scaiicnng  from  some  vccciaiion  samples."  IEEE  Trans.  Giosi  i.  Remote 
.Sefo/ny.  sol  y*.  pp  799  Sri.s,  Nov.  19S8. 


(4]  S.  L  Durden,  J.  J.  via  Zyl.  ind  H.  A  Zcbkei,  “Modeling  and  obiervaiior 
of  the  ndir  polarization  aignatuit  of  forested  areas,*  /£££  Tni, 
CeosCL  Remote  Sensutg,  vol.  27,  pp.  290-301,  May  1989. 

(5]  H.  C.  van  de  Hulst,  Light  Scaaeruig  by  Smalt  Partieles.  New  York. 
Dover,  1981. 

(6]  R.  H.  Lang  and  H.A.  Saleh,  “Microwave  inversion  of  leaf  area  and 
inclination  angle  distributions  from  backscattered  data."  IEEE  Trans. 
CeoscL  Remote  Sensing,  vol.  GE-23,  pp.  685  -694,  Sep.  1985. 

(7]  David  M.  Le  Vine,  A  ^hneider,  R.  H.  Lang,  and  H.  C.  Carter,  “Scat¬ 
tering  from  thin  dielectric  disks,"  IEEE  Trans  Antennas  Propagai., 
vol  AP-33.  pp.  1410-1413.  Dec.  1985. 

(8]  L.  Tsang,  M.  C.  Kubasci,  and  i.  A.  Kong.  “Radiative  transfer  theory  for 
active  remote  sensing  of  a  layer  of  small  ellipsoidal  scatterers,"  Radio 
Sci..  vol.  16,  no.  3.  pp.  321-329,  1981. 

(9]  M.  A.  Karam  and  A.  K.  Fung.  “Electromagnetic  scattering  from  a  layer 
of  finite  length,  randomly  oriented,  dielectric,  circular  cylinders  over  a 
rough  interface  with  application  to  vegetation,'  Ini.  J.  Remote  Sensing. 
vol.  9.  no.  6.  pp.  1109-1134,  1988. 

[10]  B.  Mandelbrot.  Fractal  Geometry  of  Nature.  New  York:  W.  H.  Free¬ 
man  Company,  1983. 

Ill]  E.  C.  Jordan  and  K.  G.  Balmain.  Electromagnetic  Waves  and  Radiating 
Systems  Englewood  Ciffs,  NJ:  Prentice-Hail.  1972. 

[12]  G.  L.  Hightshoc,  Native  Trees  Shrubs  and  Vines  for  Urban  and  Rural 
America.  New  York:  Van  Nostrand  Reinhold.  1988. 

[13]  M.  Neushul,  Botany.  Santa  Barbara.  CA:  Hamilton  Publishing  Com¬ 
pany,  1974.  p.  471. 

[14]  S.  H.  Yueh.  “Electromagnetic  and  suiistical  models  for  polarimetnc  re¬ 
mote  sensing  of  vegetation,"  Ph.O.  dissertation.  Massachusetts  Institute 
of  Technology,  Feb.  1991. 

|1S]  R.  H.  Lang.  “Bectromagnetic  backscanering  from  a  sparse  distribution 
of  lossy  dielectric  scatterers,'  Radio  Sci..  vol.  16.  pp.  15-33.  1981 

(16]  T  Le  toan.  A.  Beaudoin.  D.  L.  S.  Chong,  J.  A.  Kong,  S.  V.  Nghiem.  and 
H.  C.  Han.  “Sludv  of  microwave  interaction  with  the  Earth's  surface." 
ESA  Rep  ESTEC/No.  8447/89/NLPB(SC),  1990. 

(17]  F  T.  Ulaby  and  M.A.  El-Rays.  "Microwave  dielectric  spectrum  of 
vegetation — Pan  II:  Dual-dispersion  model."  IEEE  Trans.  Ceosci.  Re- 
moie  Sensing,  vol.  CE-25.  no.  5,  pp.  550-557.  1987. 

(18]  M  T  Hallikainen.  FT.  Ulaby,  M.A.  El-Rays,  and  L.K.  Wu.  “Mi¬ 
crowave  dielectric  behaviour  of  wet  soil — Pan  1:  Empirical  models 
and  experimenul  observations,"  IEEE  Trans.  CeoscL  Remote  Sensing. 
vol.  CE-:3,  pp.  25  -  34.  Jan.  1985. 

(19]  T  Oguchi.  “Attenuation  and  phase  rotation  of  radio  waves  due  to  rain: 
Calculation  at  19  3  and  34  8  GHz."  Radio  Sci.  vol.  8.  pp.  32-38.  19'’3. 

(20]  A.  Ishimaru  and  C.  W.  Yeh.  “Matrix  representation  of  the  vector  radia¬ 
tive  transfer  equations  for  randomly  distributed  nonsphcrical  particles." 
J.  Opt,  Soc.  Am,  vol.  1.  pp.  359-  364,  1984. 


Simoa  H.  Yueh  received  the  B.S.  and  M.S.  degrees 
in  electrical  engineering  from  the  National  Taiwan 
University,  Taiwan,  in  1982  and  1984,  respectively, 
and  the  Ph  D.  degree  in  electrical  engineering  from 
the  Massachusens  Institute  of  Technology,  Cam¬ 
bridge.  MA,  in  1991. 

He  is  now  on  the  technical  staff  in  the  SAR 
Systems  Development  and  Processing  Group  at  the 
Jet  Propulsion  Laboratory.  His  fields  of  interest  are 
electromagnetics  and  remote  sensing. 


J.A  Kong  (S'65-M'69-SM'74-F'851  is  Profes¬ 
sor  of  Electrical  Engineering  and  Chairman  of 
Aea  IV  on  Energy  and  Electromagnetic  Systems 
in  the  Department  of  Electrical  Engineering  and 
Computer  Science  at  the  Massachusetts  Instinite 
of  Technology  in  Cambridge.  MA.  His  research 
interest  is  in  the  field  of  electromagnetic  wave 
theory  and  applications.  He  has  published  7  books 
and  over  300  refereed  journal  and  conference 
papers,  and  is  ihe  Ediior  of  the  Wiley  Senes  in 
Remote  Sensine.  and  (^hief  Editor  of  the  F.lsesier 


N'nk  scries  of  Progress  In  riceironi.ipnclics  Resc.treh  (PIER). 


402 


IEEE -nUNSACnONS  ON 


AND  UMOIE  SENSBW,  VOL  30,  m  3.  MAKCH  tm 


jMKtaf  JA*(S’70-M’71-S’72-M’75)w»thor« 
on  October  31.  1943  in  Chim.  He  received  tbe 
Beccalmireaie  degree  in  1968  from  tbe  DepnitmeM 
of  Electrophysics  of  Ntlional  Chino-Tung  Univcr- 
«iy.  Ihiwin  and  tbe  Ph.D.  degree  in  1973  from  tbe 
Depanmem  of  Electrical  Engineering  and  Computer 
‘  ''=^7  Science  of  Massachusetts  Institute  of  Technology, 

Since  then,  he  has  worked  in  Riverside  Reseat 
institute  and  Raytheon  Company  before  joiiting  tbe 
MIT  Lincoln  L^ntory  as  a  technical  staff  mem¬ 
ber.  He  has  done  substantial  work  in  the  studies  and 
measurements  of  optical  and  radar  systems.  His  current  activities  concentrate 
on  airborne  radar  systems  and  clutter. 


Thaf  Ln  1hn«  aaeeivcd  ibr  CngiiMr  dcgnn  mi 
tbe  Ph.D.  degree  in  physics  ftona  tbe  Paul  Sabatier 
Uaivcniiy.  Toulonne,  France. 

In  1973,  sbe  joined  the  OeMie  d’Emde  Spaiiaie 
des  Rayonnemenis,  where  sbe  is  prcscnily  leader 
of  the  ReiiMtt  Sensirig  GroopL  Since  1980,  her 
research  activity  has  been  in  the  area  of  microwave 
remote  sensing  applied  to  vegetation.  Her  research 
interests  include  experimenution  and  modeling  of 
microwave  interaction  with  igriculturai  and  (chested 


media,  and  analysis  of  SAR  images.  She  has  served 
as  the  Project  Coordinator  and  Principal  Investigator  of  several  microwave 
airborne  experiments  in  Europe  and  is  the  Principal  Investigator  of  ERS-1 
and  SIR<r  projects. 


Robert  T.  Shin  (S  82-M’83-SM-90)  received  the 
B.S.,  M.S.,  and  Ph.D.  degrees  in  electncal  engineer¬ 
ing  from  the  Massachusetts  Institute  of  Technology 
in  1977,  1980,  and  1984,  respectively. 

He  has  been  with  the  Air  Defense  Techniques 
Group  at  the  MIT  Lincoln  Laboratory  as  a  research 
staff  member,  since  1984,  and  as  a  senior  staff 
member,  since  1989.  His  research  interests  are  in 
the  areas  of  electromagnetic  wave  scattering  and 
propagation  and  theoretical  model  development  and 
data  interpretation  for  microwave  remote  sensing. 
He  is  the  coauthor  of  Theory  of  Microwave  Remote  Sensing  (Wiley.  1985). 

Dr.  Shin  is  a  Senior  Member  of  IEEE  and  a  member  of  The  Electromag¬ 
netics  Academy.  aGL,  Tau  Beta  Pi.  Eta  Kappa  Nu,  and  Commission  F  of  the 
fniernational  Union  of  Radio  Science. 


THEORETICAL  MODELING  FOR  PASSIVE  MICROWAVE  REMOTE  SENSING  OF 

EARTH  TERRAIN 
by 

J.  A.  Kong 

Center  for  Electromagnetic  Theory  and  Applications 
Research  Laboratory  of  Electronics  and 
.  Department  of  Electrical  Engineering  and  Computer  Science 
Massachusetts  Institute  of  Technology 
Cambridge.  Massachusetts  02139,  USA 
TEL:  617-253-5625 
FAX:  617-253-0987 


Abstract 

In  this  paper  we  address  the  issue  of  theoretical  modeling  for  passive 
microwave  remote  sensing  of  terrain  media  such  as  snow,  ice,  vegetation, 
and  periodic  surfaces.  Historically,  the  volume  scattering  effects 
stimulated  the  development  of  the  continuous  random  medium  model  and 
the  random  discrete  scattering  model  for  the  description  of  the  media. 
Theoretical  treatments  were  developed  along  two  different  paths. 
Invoking  the  principle  of  reciprocity,  the  wave  theory  based  on  Maxwell's 
equations  has  been  used  to  calculate  the  emissivity.  The  other  approach 
was  to  start  with  the  radiative  transfer  equations  and  solved  for  the 
brightness  temperatures  directly.  Attempts  have  been  made  to  derive  the 
radiative  transfer  theory  from  the  wave  theory.  At  the  same  time,  both 
theoretical  approaches  have  been  used  to  calculate  the  radiometric 
emissions  and  to  interpret  experimentally  measured  data. 

The  successful  interpretation  of  the  Cosmos  243  data  was  perhaps  the 
first  most  important  step  towards  a  serious  development  of  the 
continuous  random  medium  model  to  account  for  the  volume  scattering 
effects  of  snow  ice  fields.  Subsequent  interpretation  of  measurement 
results  from  snow  field  with  both  passive  radiometers  and  active  radar 
systems  established  a  unique  position  for  its  description  of  earth  terrain 
media.  Recent  efforts  in  classifying  sea  ice  with  correlation  function 
characterization  are  demonstration  of  acceptance  of  this  model.  Future 
inverse  scattering  developments  will  perhaps  rely  heavily  on  this  model. 
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Abstract 

Earth  terrains  are  modeled  by  a  two-layer  configuration  with  a  random  permittivity  deKribed 
by  a  three-dimensional  correlation  function  with  horixontaJ  and  vertical  correlation  lengths  and 
variances.  Using  the  wave  theoo'  *ith  Bora  approximations  carried  to  the  second  order,  this  model  is 
applied  to  derive  the  polarimetric  backKattering  coefficients  of  the  Mueller  and  covariance  matrices. 
From  a  physical  point  of  view,  the  Born  first-  and  second-order  approximations  account  for  the  single 
and  double  scattering  processes,  respectively. 

For  the  isotropic  random  permittivity  configuration,  five  out  of  the  nine  elements  of  the 
covariance  matrix  are  tero  under  the  Born  first-order  approximation.  For  the  uniaxial  tilled  random 
permittivity  case,  the  covariance  matrix  does  not  contain  any  xero  elements.  To  account  for  the 
azimuthal  randomness  in  the  growth  direction  of  leaves  in  tree  and  grass  fields,  an  averaging  scheme  is 
developed  in  which  the  backscattenng  coefficients  are  averaged  over  the  azimuthal  duection.  In  this 
case,  the  covariance  matrix  is  described  by  four  zero  elements  though  the  tilt  angle  is  different  from 
zero.  Applying  (he  Born  approximation  to  the  second-order,  the  covariance  matrix  is  computed  for 
(he  isotropic  and  the  uniaxial  untilted  random  permittivity  configurations.  The  covariance  matra 
has  four  zero  elements,  and  a  depolarization  factor  is  obtained,  even  for  the  isotropic  case. 

To  describe  the  effect  of  the  random  medium  on  electromagnetic  waves,  the  strong  permittiv¬ 
ity  fluctuation  theory,  which  accounts  for  the  losses  due  to  both  of  the  absorption  and  the  Kittering, 
is  used  to  compute  the  effective  permittivity  of  the  medium.  For  a  mixture  of  two  components,  only 
the  frequency,  the  correlation  lengths,  the  fractional  volume,  and  the  permittivities  of  the  two  con¬ 
stituents  are  needed  to  obtain  the  polarimetric  baefcseattering  coefficients.  Theoretical  predictions 
ate  illustrated  by  comparing  the  results  with  experimental  data  for  vegetation  fields  and  sea  ice. 
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Introduction 

For  many  types  of  earth  terrain,  the  scattering  effects  due  to  medium  inhomo¬ 
geneities  play  an  important  role  in  the  determination  of  radar  backscattering  coeffi¬ 
cients.  The  volume  scattering  properties  of  a  medium  have  been  successfully  modeled 
by  the  random  permittivity  fluctuations  [Ij.  The  random  medium  is  characterized  by 
a  background  permittivity  and  a  random  permittivity  fluctuation  modeled  by  a  three- 
dimensional  correlation  function  with  horizontal  and  vertical  correlation  lengths  and 
variances. 

The  objective  of  this  paper  is  to  develop  a  mathematically  rigorous,  fully  po- 
larimetric  model  to  compute  the  covariance  matrices  in  the  backscattering  direction 
for  different  kinds  of  earth  terrains  with  a  two-layer  configuration  medium  vnth  ei¬ 
ther  an  isotropic  or  an  anisotropic  random  permittivity.  Of  all  the  analytical  wave 
approaches  to  the  study  of  scattering  by  random  medium,  the  Born  approximation 
is  the  simplest.  An  integral  equation  is  first  formed  for  the  electric  field  using  the 
unperturbed  Green's  function  in  the  absence  of  permittivity  fluctuations.  The  inte¬ 
gral  equation  is  then  solved  by  iteration  with  the  iteration  series  knoam  as  the  Born 
series.  With  both  the  first  and  second  terms  of  this  series  which  respectively  describes 
a  single  and  a  double  scattering  processes,  the  polarimetric  backscattering  coefficients 
are  computed. 


Scattering,  Mueller,  and  Covariance  Matrices 

Consider  an  electromagnetic  plane  wave  propagating  in  the  direction  k. 

£  =  (£,A +  (1) 


where  Eh  is  the  horizontal  component.  £,  is  the  vertical  component,  h  and  f  are  tw’o 
orthogonal  unit  vectors  satisfying  the  relation  c  x  ft  =  k.  Consider  an  incident  plane 
wave  £,  impinging  on  a  particle.  This  particle  will  radiate  and  the  scattered  electric 
field  E,  can  be  related  to  the  incident  electric  field  £,  by  the  scattering  matrix  [2]; 


(Eh.\_(HH  HV\  (EkA 
\E.J~\VH  VVJ'KE^J 


(2) 


The  elements  of  the  scattering  matrix  are  complex  and  the  relation  VH  =  HV  holds 
[1]  in  the  backscattering  direction  and  for  reciprocal  media. 

Instead  of  relating  the  scattered  electric  field  to  the  incident  electric  field,  the 
modified  Stokes  vector,  defined  by  [1],  can  also  be  used  to  relate  the  scattered  Stokes 
vector  to  the  incident  Stokes  vector: 

/«  =  —  /v=—  C'  =  -£e(£,£;)  V  =  -/m(£*£;) 

I?  n  tj  V 


(3) 


II6I 


This  definition  is  chosen  since  the  first  two  modified  Stokes  parameters  represent  the 
iiitensit}'  of  the  wave  in  the  horizontal  and  vertical  polarizations: 


In  terms  of  the  elements  of  the  scattering  matrix  (2),  the  Mueller  matrix  M  relates 
the  scattered  to  the  incident  Stokes  vectors  according  to 

1.  (5) 

is  found  to  be 

IWV'I’  Re(HHHV] 

\VH\^  |VV|’  Re(VH  VV]  -Im(VH  VV') 

2Re{HHVH‘)  2Rt(HVVV)  Re(HJi  ■  VV  ^  HV  VH’ )  VV'  -  HV  V  H‘ ) 

2lm{HH  VH')  2J7n(HV  VV)  lm(HH  VV'  ■¥  MV  VH*)  Rt(HH -W  -  HV  V H‘ ) 

(6) 

where  all  the  Mueller  matrix  elements  are  real  numbers.  From  the  four  ele¬ 
ments  of  the  scattering  matrix  (2),  the  covariance  matrix  may  also  be  defined.  In 
general  the  covariance  is  a  complex  matrix  of  sixteen  elements.  In  the  special  case 
of  the  backscattering  direction  and  for  a  reciprocal  medium,  the  covariance  matrix  is 
transformed  to  a  matrix  of  nine  elements  as  follows 

^  (HH\  {  \HH\^  HH  HV  HH  VV\ 

C==\hv\{HH’  HV  VV)-\HV  HH'  \HV\^  HV  •  VV 
\VV;  \VV  HH'  VV  HV'  |VV?  j 

(7) 

By  comparing  the  scattering  (2).  the  Mueller  (6).  and  the  coN’ariance  matrices  (7). 
it  is  easily  seen  that  all  the  polarimetric  information  may  be  expressed  in  any  of 
these  matrices.  The  covariance  and  Mueller  matrices  can  be  transformed  to  each 
other  since  they  contain  the  same  information.  However,  from  these  two  matrices, 
the  scattering  matrix  cannot  be  fully  retrieved  because  the  absolute  phase  factor  is 
not  taken  into  account  in  either  the  covariance  matrix  or  in  the  Mueller  matrix 


In  our  study  of  microwave  remote  sensing,  the  random  medium  model  is  used 
to  simulate  the  permittivity  of  the  earth  terrain.  This  implies  that  the  simple  case 
of  a  single  panicle  can  no  longer  be  considered  but  the  ensemble  average  of  the 
electric  field  needs  to  be  calculated.  The  polarimetric  bistati-  scattering  coefficient 
is  defined  by  (3j 


—  hm 


4trr* 

A  cos  6o, 


Er,E:. 


lim 

>4 


4»rr^ 

A  cos  6o, 


<  fiT  ■  vk‘  >  (8) 
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where  A  is  the  area  illummated  by  the  radar,  r  is  the  distance  between  the_radar 
and  the  observed  terrain,  and  Soi  »  the  polar  incident  angle.  Purthetmore,  fcj  and 
k,  are  the  incident  and  scattered  wave  vectors,  respectively  The  subscripts  s  and 
i  stand  for  scattered  and  incident  waves,  respectively.  Finally.  iiT  and  vk  are  the 
four  components  of  the  scattering  matrix  where  n.  t.  v,  and  k  represent  either  the 
horizontal  h  or  vertical  v  polarizations. 

The  polarimetric  backscattering  coefficient  is  defined  13)  from  the  bistatic 
coefficient  when  k,  =  -k. 


(k,. 


k,)  cos  00, 


E  E’ 


<  /IT  -  l/K*  > 

(9) 


The  polarimetric  backscattering  coefficients  are  complex  except  the  three  "tra¬ 

ditional”  real-valued  backscattering  coefficients 


OhhKh  Okh 

Okxkx  =>  ^kx 

Oxxxx  ^  ^XX  (101 

Similarly,  the  covariance  matrix  (7)  constituted  by  the  different  polarimetric  backscat¬ 
tering  coefficients  can  be  derived 

,  2  ( (  ^kk  Okkkx  Okkxx\ 

e  =  lim  ^  \  HV  )  =  oj**.  o,.  (111 

in  ^  [Vvj  O'kxrx  ‘’rx  I 

Therefore,  the  complete  polarimetric  information  can  be  obtained  by  computing  only 
six  backscattering  coefficients  in  which  three  are  real  (o*/,.  o*,.  and  a„}  and  three 
are  complex  (OkkkT-  (^kkrr-  Okxxx)- 


Formulation 

Consider  an  electromagnetic  plane  wave  impinging  upon  a  two-layer  random 
medium  configuration  as  shown  in  Fig.  1.  The  top  layer  represents  free  space  with 
a  permittivity  <o  while  the  bottom  layer  has  a  deterministic  permittivity  In 
between,  a  layer  with  thickness  d  and  extending  infinitely  in  the  lateral  direction  is 
characterized  by  a  random  permittriity  <i(r) 

<l(?)  “  Clm  -t-  Cl/(f)  (12) 

This  random  variable  has  a  mean  %-alue  and  a  small  fluctuating  term  whose 
average  is  zero.  The  wave  equations  for  both  regions  may  be  written  as 

V  X  V  X  £o{f)  -  ko^o(r)  =  0 


(13) 


V  .  r  X  £)(r)  -  ^,„£,(r)  =  P(r)£,(r)  (14) 

'^aere  £o(^)  iad  £nf;  denote  the  electric  fields  fluctuation  in  region  0  and  region  1, 
respectively. 


2 


Layer  With  Random  Permitlivrty 


2  « -  d 


Figure  1;  Isotropic  two-layer  configuration. 


Furthermore,  the  follottnng  variables  have  been  defined 

^'o  =  <?(f)  =  )  (15) 

With  the  effective  source  term  (?(r)£i(r),  (13)  and  (14)  can  be  expressed  in  integral 
form  as  follows 


£o(f) 

=  ^0 

J 

1  d^ri  Go:(f,f,) 

!<?(?!  ) 

■  £'i(’'i) 

(16 

^l(c) 

0  J .  _  . 

=  f,  (r)-  ^ 

f  Gjilr.ri) 

Q(r,) 

(17 

where  £i'’\r)  and  £','”{?)  are  th**  specular  or  co^herent  electnc_fields  in  the  absence  of 
random  permittivity  in  region  1.  Furthermore  Goi{f,fj)  and  Gii(r,fi)  are  the  dyadic 
Green's  functions  for  observers  in  regions  0  and  1,  respectively,  with  the  source  in 
region  1.  Solutions  to  these  two  coupled  equations  can  be  obtained  by  iteration 
Substituting  (17)  into  (16),  the  total  electnc  field  in  region  0  is  found  to  be  in  form 
of  the  Neumann  series 
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where  the  n'*  order  field  is  given  by 

=  J d’n-  5o.(f.?,)<?(?i)  -5„(l',,T,)C?(M  -5n(f.-,,F.)<?(f,)  r<,®Vn)  (19) 

The  so-called  “Born  first-order  approximation"  is  obtained  by  letting  n  =  1  in 
the  previous  equation.  Hence 

E['\f)  =  J  d'f,  Goi(r.r.)  (20) 

and  the  “Born  second-order  approximation”  is  derived  by  letting  n  =  2 

fo'V)  =  /rf'n  5o,(f.f, )<?(?,)  S„{f,,f,)  (21) 
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Born  First-Order  Approximation 

a)  Isotropic  random  medium  model 


Consider  an  electromagnetic  plane  wave  impinging  upon  a  two-layer  random 
medium  configuration  as  shown  in  Fig.  1.  For  an  incident  electric  field  in  region  0  of 
the  following  form 

£o.(f}  =  e’*'-"  (22) 

where  h{-ko„)  and  if-ton)  respectively  the  horizontal  and  vertical  unit  vectors 
in  region  0  associated  with  down-going  waves,  the  coherent  electric  field  in  the  random 
medium  is  obtained  j4j 

^i\r)  =  {£o*.  -1- 

-f  £o„  (X..t;(il-u.)e**''‘'  +  (23) 

in  which  the  following  coefficients  are  defined  j3] 


a.  —  o,, 

^2i 


Z?2,  =  1  -t-  R(iu£\7>  e'^*‘*’^ 

^  kcii 

~  7  .  I  I. 

'*€11  ' 


m  *^121  ^ 

A:i  F2, 

n  ^Olt 

£21  =  1  +  £oif5i2i  e’****'** 

«*•  —  J  ^7  J 

+  <oKk»i 

1  «Si  “  1  "t* 


As  can  be  seen  in  (23),  the  A  coefficients  are  associated  with  the  up-going  waves  while 
the  B  coefficients  describe  the  down-going  waves  Under  the  far-field  approximation, 
the  Green's  function  in  (20)  takes  the  following  form  in  the  backscattering  direction 

!5l 

Co)(r,fi)  =  h(-Uo»i)vjoiA(r.rj )  —  i'(-Ioii)Coi»{r. f, )  (25) 

Co,*(r.f.)=^  U/..MU,=.)e'*“-'=  (26) 

■lirr  t  J 


Gou(r,f,)=  ^  [.4..r(A:i,.)c’*>-^-  +  e’*- ' 


For  each  polarization,  both  Goi(r,fi)  and  £^i\f)  are  constituted  by  a  dol^•n-  and  an 
up-going  waves.  Therefore,  £o**(r)  is  formed  by  the  sum  of  four  different  terms  as 
shown  in  Fig  2.  All  the  multiple  reflections  occurring  at  the  boundaries  are  also 
incorporated  in  the  model  and  the  backscattered  wave  is  due  to  the  single  scattering 
process  under  the  Bom  first-order  approximation. 
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Figure  2;  Physical  interpretation  for  Born  first-order  approximation 

To  compute  the  polarimetnc  backscattering  coefficients  (9).  only  the  components 
of  the  backscattered  electric  field  are  needed.  Furthermore,  the  polarization  of  the 
incident  wave  is  only  along  one  direction,  either  horizontal  or  vertical  With 

the  Fourier  transform  of  the  autocorrelation  of  the  permittivity  fluctuation  for  a 
statistical  homogeneous  random  medium 

<  Q{h)Q(f2y  >=  r  (28) 

V-oc 


«r- 
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the  ensemble  average  for  the  first-order  backscattered  intensity  takes  the  following 
form 

So,,(f.rO  .  £<“’(?,)]  [5oi.(r.r,) -  £!“’(?:)]’  (29) 

where  r  and  k  represent  either  the  horizontal  or  vertical  component  of  the  transmitted 
wave  and  //  and  u  represent  either  the  horizontal  or  vertical  component  of  the  received 
wave.  From  (23),  the  following  term  has  been  defined 

=  E  a(pl:,„)e-'*>-^e**'' (30) 

p=-i 

where  a  stands  either  for  h  or  v,  and  p  describes  either  up  or  down-going  waves.  In 
a  similar  manner, 

=  ^  t,  /*'•">  (31) 

where  the  following  definitions  have  been  used 

A"ai,  =  j4*,  =  B*,  =  Ax,  =  Bx,  (32) 

Substituting  (30)  and  (31)  into  (29).  a  sum  of  16  terms  is  obtained 


<  B;i’(r)£'i’(r)*  >  =  ! 


■£ot,£o„  ^  A(p^'ui)  • 


(AV,.A'„.  u(rki„)  «(sl,,„)]’  j  d^fid^fj  j  d^3  *(J) 

As  can  be  seen,  all  the  polarization  information  is  contained  outside  the  integral. 
Also,  the  intensity  is  expressed  as  a  sum  of  sixteen  terms  since  the  electric  field  is 
composed  of  four  elements.  With  the  two  foDowing  scalar  quantities 

*«’•  _  Y  V  Y*  V*  (341 

j’(9li;.)i  [i'(rki„)  k\ski,,}'‘  (2b) 

(33)  can  be  wntten  as 

<  B;'’(r)£ii’(f)*  >  =  E  ^  *1:^'-  ^7::^  J 
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The  integrals  in  (36)  are  solved  by  decomposing  the  variables  of  integration  into  their 
lateral  and  vertical  components.  By  using  the  property  of  the  Delta  function  and 
introducing  the  illuminated  area  A,  the  integration  over  the  lateral  variables  p,,  P2< 


and  Jp  yields 


<  E  4';r«  n«:*.  Jdj.dzidz, 

In  order  to  solve  the  last  set  of  equations,  let  us  define  first  the  two  complex  quantities 


a  —  pk'i-,  + 


b  =  [rku,  +  sku.]’ 


Hence,  the  integral  pan  of  (37)  may  be  aiitten 


I dz.dz^  e“---  e-'‘«  j d3.  =  2k^,3.)  (39) 

Due  to  Cauchy's  theorem,  the  integral  over  3:  >ields  the  following  expression  by 
distinguishing  two  cases  depending  upon  the  sign  of  zi  ~  Z2 

f  d3.  i[2k^,3.)  if  rj  <  ^ 

^  l-2jri  fles4>(2JLv, ,3; le--*- <*>-■>'  if  r,  > 

where  it  is  assumed  that  the  spectrum  density  )  has  only  one  pair  of  complex 

conjugate  poles  at  3..  =  3f.  and  i(3^.3,)  t-anishes  everj-where  on  the  infinite  circle. 
For  the  random  medium  of  infinite  lateral  dimension  and  vertical  extension  from 
r  =  -  d  to  r  =  0,  (39)  becomes 

2tri  Res^J^  =  2k^,3:)  f^dzij^dz^  e»>(«-'87)e--->(t-3r) 

-  2x1  Res^l^  =  2k^,3:)  /’/-‘/'j  e*--. (•-«.- )e-.*>(»-^- )  (41) 

For  ResO^  =  2k^,3~)  =  -Res(5,  =  2kf„3*],  these  integrals  may  be  easily  solved 
to  get 


2Trj  Resi{3p  =  2kp,.3*) 


j  _  1  _ 

(a- 3r){b- 3:)'  (a-b)(b-3:) 


^(a-b){b-3:)  ■  {a-0:){b-3:)  ^  {a~0:){b-3:) 
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To  simplify  the  expression,  let  us  define  the  following  term 

"  [(a  -  "  {a-b){b-l3t)  (a-  b)(b-  ^7) 

-  (a-  i3:)(b- (a- 0:)(b-  n:) 

From  (38),  a  is  different  from  b  for  most  cases  because  of  the  complex  conjugation. 
However,  for  the  special  case  when  a  =  b  =  Q,  (43)  may  be  written  as 


Substituting  (43)  and  (44)  into  (37)  and  using  (9),  the  polarimetric  backscattering 
coefficients  are  found  to  be 

c,^.a  =  21^^^6u;*n^\(UiRes^0,^2k^,0:)  ^  («) 

Only  three  polarimetric  backscattering  coefficients  are  different  from  zero.  This  comes 
from  the  factor  (35)-  They  are  zero  due  to  the  definition  of  the  unit  vectors 

[3].  Hence,  the  final  results  are 

chh  =  iHes<t(3,  =  ^  /«'• 

c„  =  2x'8-;V'|f iRes*(3,  =  2k^.3: )  L  ^^rlr  /”'* 

r.*„  =  2,r'WiCa^l^/?es*(3,  =  2i^,^;)  f  KhW  m 


^hv  ”  ^hThr  *  ^hhhv  *“  ~  0 

For  an  autocorrelation  function  for  the  random  permittivity  (12)  of  the  form 


_  iti-t 


<  <i/(fi)ei/(?j)*  >=  c"  '»  c'  ''  e' 


where  f  is  the  variance,  and  1,  are  the  lateral  and  vertical  correlation  lengths, 
respectively,  the  spectrum  density  (28)  is  given  by 


(49) 
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The  spectral  density  goes  to  zero  as  ^  tends  to  infinity  satisf)'ing  the  requirement  iur 
the  complex  integrauon  (40).  The  poles  0,  of  the  spectrum  density  are  cbtained 


3*  = 


t 


Furthermore,  the  associated  residues  are  found  to  be 


Res^(3,,3y,3f)  =  ± 


‘I 

2ir3i  (l+ 5.^/7)  (l+JJ/2) 


50' 


‘  51 


b)  Anisotropic  random  medium  model 

In  order  to  describe  earth  terrains  exhibiting  an  anisotropic  behavior  the  ran¬ 
dom  permittivity  is  chosen  to  be  uniaxial  in  a  tilted  coordinate  system  ,'6j  Earth 
terrains  such  as  sea  ice  and  com  fields  display  such  a  property.  From  (12),  the  random 
permittivity  is  tiiitteu  in  a  tensor  form 

^l(r)  =  ^Im  +Il/{f)  ;52) 

where  ?i/(r)  is  a  function  of  position  characterizing  the  randomly  fluctuating  com¬ 
ponent  of  the  permittivity.  For  the  random  medium,  it  is  assumed  that  both 
and  ?i/(r)  are  uniaxial  in  a  coordinate  system  ly'z'  lilted  respectively  at  an  angle  v. 
t'f  around  the  r  axis  as  represented  m  Fig.  3.  In  this  tilted  coordinate  system,  the 
permittivity  tensor  takes  the  following  form 

(u  0  0\ 

lim  =  0  f,  0  (53) 

Vo  0  ful 

In  the  coordinate  system  xy:,  the  permittivity  tensor  can  be  expressed  as 

/eu  0  0  \ 

tlm  ~  I  0  <22  <23  I 

V  0  <32  <33/ 

where 

<11  =  <1 

<22  =  <1  COS*  W  +  €i,  sin*  t’ 

<23  =  <32  =  (<w  -  <i)sin  VCOSV 
<33  =  <1  sin*  +  <1,  cos*  v 


(54) 


(55) 


r- 


IP, 


7 


Figure  3;  Geometr;ca]  configuration  of  the  permittiMty  tensor  in  an  amsotrjp.t 
randon:  medium. 


Similar  equations  may  be  written  for  ?i/(r)  with  the  tilt  angle  V/  From  (20).  the 
“Born  first-order  approximation"  yields  for  the  anisotropic  case 

^o''(n  =  / dh,  (56) 

where  Q(f)  =  u,’-/i(i/(f ).  Using  Euler  notation  and  decomposing  the  electric  field  in 
horizontal  and  vertical  components,  the  ensemble  average  for  the  first-order  backseat- 
tered  intensity  m  ly  be  written  as 

<  £o|,*(T)£iJ,'(f)’  >  =  Y.  f  [Coi^;(f,T,)£;“*(f,)] 

[Goiw(?.  >  (oT  i 

where  the  average  of  the  components  of  the  electric  field  has  been  computed  by  defin¬ 
ing  ij  and  1/  as  the  received  polarizations  and  r  and  k  as  the  transmitted  polarizations 
As  shown  in  Fig  4,  a  horizontally  or  a  vertically  polarized  incident  electric  field 
will  generate  both  an  ordinary  and  an  extraordinary  waves,  due  to  the  anisotropic 
permittivity  in  region  1,  defined  by  their  respective  wave  numbers  [6] 


k^,  =  -  kl 
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=  -  —ky,  +  lit  it  1,(33  -  f  1(33^1.  -  (Kukl. 

*33  fj3  * 

^W.  =  -  ~^v.  -  ~  ^/•*>^^ttl(  1,(33  -  (itszkli  -  fjfui'J, 


(58) 


a) 


z 


b) 


z 


Figure  4:  Ordinary  and  extraordinary  electromagnetic  waves  excited  by  the  ani¬ 
sotropic  random  medium,  (a)  Horizontally  polarized  incident  electric 
field,  (b)  Vertically  polarized  electric  field. 
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The  components  of  the  unperturbed  electric  field  E*,***(f)  and  the  Green’s  function 
Ggi(r,ri)  may  be  written  as  (7] 

rl«(r)=  E  '  (59) 


Goio(r.fi)=^  E  (60) 

where  o  stands  for  either  A  or  t,  and  p  describes  either  an  ordinary  up-goii^g  wave  (o). 
an  ordinary  down-going  wave  (-o).  an  extraordinary  up-going  (eu).  or  an  extraordi¬ 
nary  down-going  wave  (ed).  Both  the  incident  electric  field  and  the  Green’s  function 
are  formed  each  by  four  terms  meaning  that  backscattered  electric  field  computed 
with  Born  first-order  approximation  will  be  constituted  by  sixteen  terms.  Further¬ 
more,  the  coefficients  A,,,,,  the  unit  vectors  pfH*?,,),  and  the  wave  numbers  Jlf„  take 
on  the  following  expressions  depending  on  the  value  of  p 


p  —  0  =»  Af.,  —  kjj; 

P(kL)  = 

P=  ~o^  irf.,  =  -k°,i 

Xapi  —  Boot 

p  =  eu  •  Af..  =  AJ* 

p{kL)  =  «(Jtn.) 

^api  *  ^o«i 

p  =  ed=>.  k^„  =  k\< 

Pikli)  =  iikrii) 

^api  Bf^ti 

(61) 

Hence,  substituting  (59)  and  (60)  into  (57)  yields 


uiVVil’ 


Y. 


lA^p.A',^  P( *•{’,.);  •  9(*i..)i.]  [A^iA'„.  ■  s(*I.i)m]‘  j 


(62) 


As  can  be  seen,  all  the  polarization  information  is  contained  outside  the  integral 
The  parameters  p.  q.  r.  and  s  can  take  the  values  o,  -o,  eu,  or  ed  as  described  in 
(61).  The  intensity  is  hence  given  b>-  a  sum  of  256  terms  For  each  of  these  terms.  81 
products  (sum  over  jklm)  have  to  be  computed  due  to  the  anisotropic  beha^^or  of 
the  random  medium.  In  order  to  simplify  this  result,  let  us  define  the  two  following 
scalar  quantities 

.x:.,  (63) 


=  A,pA,^A'l 


=  lF(A-f,.),  ■  diklM  ir(i[..),  •  Hkuxr 


(64) 
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Hence,  (62)  may  be  written 

^  ^Om  \^)  16jr"r^  ^or*C(j^ 


Eor,EL  z  ii  ♦r-. 


(65) 


By  following  a  similar  procedure  as  in  the  isotropic  case  described  in  (37)  through 
(45),  the  parameters  a  and  b  are  first  defined  as 


^  =  kl,  +  i-?,. 


b  =  -  klX 


(66) 


By  assuming  a  spectrum  density  <P(J^,J^)  with  n  pairs  of  complex  conjugate  poles 
at  3i  =  3fn  defining  the  followng  factors 


/r;*  = 


1  _  g-'io-^rn)^ 


J  _ 


(a-J;j(b-3;j  (a-b)(b-J;j 


/nT'  = 


1  _ 


f '(»-«.■.)- 


-  i)  (6  -  /3.-n)  (a  -  5rj  (6  -  ;3- )  (o  -  3:„){b  -  JrJ 
the  six  polarimetric  backscattering  coefficients  (11)  are  finally  found  to  be  [8] 


(67) 


(68) 


_  0-2  ,  2.,  2 
^ itrun  —  ^  ^ 


Z  k::.  t 

p,q,r,a=o  ;>./.m  =  l 

tZ  =  2k^.3:j  /«"  -  /?es«*>,n„(3,  =  2^,.,  /T*}  (69) 


Compared  to  the  results  of  the  isotropic  case  shoi^m  in  (46)  and  (47).  there  is  no  zero 
elements  This  is  due  to  the  anisotropic  behavior  and  the  excitation  of  both  ordinary 
and  extraordinary  waves. 
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Born  Second-Order  Approximation 


For  this  case,  only  the  isotropic  random  permittivity  will  be  considered.  Prom 
(21),  the  second-order  scattered  electric  field  is  given  by 

d^fj  5oi(r,  ri)  (?(fi)  -  SuCn.rj)  Qifj)  ■  E[°\f7)  (70) 

where  and  Goi(f.f,)  are  given  by  (30)  and  (31),  respectively.  The  Green’s 

function  Gi](fi,f2)  relates  the  field  observed  in  region  1  at  a  position  fj  due  to  a 
source  in  region  1  located  at  fj.  It  takes  the  following  form  from  (9) 


tih~h> 


I  -'1  >  -3 

1  21  <  i. 


where 


t,  e”*‘"‘e''*‘--’  Yi  6(9k,,)6(-rlt,,)  (72) 


5n{i;„^i,22)  =  ~  E  E  a(qkM-rku)  (73) 


Here  a  denotes  the  polarization  direction  along  h  or  v,  and  the  coefficients  C  take 
the  foUo«ing  values  [3] 


=  C*<  ,(!:,)  = 

CLiiik,)  =  C*<,.,(k,)  = 

Cn(tp)  =  =  5.je’^*‘'<'/i^2 

Cn(A>)  =  C<.i(X-,)  =  S,oSi3f'^*-V-F2 


~  ^h-li(^p)  ~  l/'^2 

^A-i-i(^e)  “  G*_,_j(A:p)  =  Rio/D^ 

^rl-l(^v)  ~  =  1/^2 

r>.,.,(Xp)  =  c<...,(Xp)  =  s.o/-F2 

(74) 


The  physical  meaning  o^the  Born  second-order  approximation  (70)  are  interpreted 
based  on  the  fact  that  Goi(r.fi).  Giilrj.fj).  and  E^^\r)  are  composed  of  a  down- 
and  an  up-going  wave  for  each  polarization  Therefore.  Eg\f)  is  formed  by  the  sum 
of  sixteen  different  terms  as  shown  in  Fig.  5.  All  the  multiple  reflections  occurring 
at  the  boundaries  are  accounted  for  and  the  backscattered  wave  is  due  to  the  double 
scattering  process  under  the  Bom  second-order  approximation. 

With  (30),  (31),  and  (70),  the  ensemble  average  for  the  second-order  intensity 
is  given  by 


<  £ip^(F)  ^^’(r)’  >=  J <  Q{f,)Q{fi)Q{f,r Q(f,y  > 
|Soip(r,ri)  ■  Gai(fi.f2)  •  £*,°’(r2)j  |5oi^(r,r3)  •  Guf/j.f^)  •  £'j°^(f4)j  (75) 
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where  ft,  v,  r,  and  k  represent  either  the  horizontal  or  vertical  polarization.  Further¬ 
more,  ft  and  V  are  associated  with  the  receiver  whereas  r  and  k  are  associated  with 
the  transmitter.  By  assuming  gaussiao  random  variables  and  considering  only  the 
backscattering  direction,  the  random  term  of  (75)  may  be  written  as  a  sum  of  two 
terms 


<  Q(r})Q(r2)Q{r3)* QiuY  >  =<  (?(ri)Q(fj)*  ><  Q{r2)Q(f^y  > 

+  <  Q(riW(r4)’  ><  <?(rj)f?(r3)*  >  (76) 


Figure  5:  Ph3'sica]  interpretation  for  Born  second-order  approximation. 


1177 


Similarly  to  (28),  by  using  tb«  spectrum  density  of  the  autocorreUtion  function  and 
integrating  over  the  lateral  variables,  the  ensemble  average  takes  the  following  form 

<  4J’(r)  4l\rr  >=  Eor.ELld%  / dz,dzrdz,dz, 

1  1 

[ dj,id0,2  +  kp,  0,1 )  0.2 ) 

J  l^ul 

g-'(e'*r.. ■*•«'**.  1*5  e"***'****"^'  **'’'*  ^oT!£itj(kp) 

where  the  following  terms  have  been  defined 

Krt':'  =  (78) 

('9) 

=  lA(pki«)  ■  MM  d(-rku)  ■  T(«ku.)] 

It>(p'*i..)  •  7(«'kK-)  M-r'ku)  ■  k(s'ki.i)]*  (80) 

firrv«f^'''''("*e)  =  iMpki..,)  •  a(^ku)  M-rku)  ■  T(ak,.i)] 

[i>{p'ku.}-i(-<l'ki..)  i{r'ku).k(s%„)]-  (81) 

In  the  polar  coordinate  system,  (77)  becomes 

<  EZ\-r)  £Z\-rr  >=  ^or.E;^.J j dz,dz,dzydz, 

kp  ^  1=1  i=j  •"'(p'**,, +«’**. )*»  l=< 

1  Z- 

I A  li :  »  «•••*  -  -  J 

-'.s'  ” 

o.^  =  A.t 

,-.5.i(--.-.-ll  +  i:!rl%V'-''^\kp)  j  dJM2  ^(k>  +  k,.0,,] 

4>(k>  -  A>/?,2)e-’^“‘“-^’  (82) 

where  the  integration  over  <p  has  yielded  [3] 

1  r  do  =  /«;;':i^’'’''(A>) 

TT  JQ 


(83) 
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7i  «0 

Using  the  Cauchy  s  theorem,  the  integration  over  J,i  and  J;j  may  he  earned  ou: 
The  final  result  for  the  poiarimetric  backscattenug  coefficients  are  obtained  3| 


iRes^(k^- k,.3:)  t,  ^  j: 

a.'y*h.v 


w‘  •«' 


The  factors  M  and  .V  are  defined  in  [3]  It  is  noted  that  the  integr<ij  ovr:  v,  vannn 
be  solved  analytically  but  only  numerically  From  the  integration  over  :.te  azimuiLa. 
e'ariable  o  in  (83)  and  (84).  it  is  found  that  Okrrr  =  OhKhr  =  0.  However.  t;ir 
depolarization  terra  tr*,  is  different  from  zero,  unlike  the  result  obtained  with  the 
isotropic  Bom  first-order  approximation. 
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Results 

a)  Born  first-order  approximation 

With  the  Born  first -order  approximation,  the  polarimetric  backscattering  coeffi¬ 
cients  have  been  obtained  for  the  two-layer  isotropic  and  anisotropic  random  medium 
configurations.  For  the  anisotropic  random  medium,  the  depolaruation  term  Ok.  t-  0 
exists  even  with  Born  first-order  approxima;.ion.  In  general,  ail  the  backscattering 
coefficients  o*/,,  o*,.  are  different  from  zero  for  the  anisotropic 

case.  However,  when  f  =  V/  =  0*  but  <i  ^  ci,,  no  depolarization  term  is  obtained  for 
this  uniaxial  untilted  case  and  o*.  =  oaaa*  =  =  0-  Also,  for  =  0  but  V/  ^  0. 

the  non- zero  depolarization  terms  and  OArtr  are  obtained  Finally,  when  C  ^  0. 
all  the  backscattering  coefEciems  are  different  from  zero.  In  a  similar  manner,  when 
<p  =  90*.  implying  that  the  plane  of  incidence  contains  the  scatterers,  there  is  no  de¬ 
polarization  factor,  and  only  three  backscattering  coefficients  are  different  from  zero 

^ kvrv  “01. 

Fig.  6  illustrates  dependence  of  the  backscattering  coefficients  on  the  angle 
of  incidence  for  Oa*,  and  o^r  for  a  two-layer  configuration,  and  these  results 
are  matched  with  expenmentai  data  for  sea  ice  taken  in  1984  in  Point  Barrcm-  in 
Alaska  by  [lOj  T*''  input  parameters  <i,  cu.  6,  and  6,  are  obtained  using  the  strong 
fluctuation  theory  developed  in  (11).  The  tilt  angle  is  chosen  to  be  v>  =  v/  =  20*, 
and  the  azimuthal  angle  is  assumed  to  be  ©  =  20*.  The  good  correspondence  is 
noted  between  the  experimental  results  and  the  theoretical  predictions.  Fixing  the 
incident  angle  at  6o,  =  40”  in  Fig  6.  the  covariance  matrix  for  the  backscattering 
coefficient  (14)  is  shown  in  Fig  7.  As  can  be  seen,  this  matru  does  not  contain  any 
zero  elements. 

For  an  earth  terrain  medium  such  as  trees  or  grass,  leaves  grow  randomly  in  the 
azimuthal  direction.  This  can  be  modeled  by  assuming  that  one  patch  illuminated 
by  the  radar  is  oriented  in  a  certain  direction  whereas  the  next  patch  is  oriented  in  a 
different  direction  and  so  on  To  obtain  the  backscattering  coefficients,  the  average 
over  many  patches  should  be  taken.  This  model  is  used  in  the  theory  by  assuming 
a  fixed  direction  of  the  scatterers  but  modihnng  the  position  of  the  radar  o\Tr  the 
azimuthal  direcrion  This  is  shown  in  Fig  8  for  grass  and  in  Fig  9  for  trees  where 
both  these  vegetation  fields  are  simulated  by  averaging  the  backscattering  coefficients 
over  o  The  step  of  averaging  is  Ao  =  1°  Though  the  tilt  angle  i-  is  different  from 
zero,  there  is  a  depolanzation  factor  (Oa,  ^  0).  but  Oaaat  =  =  0.  This  is  due  to 

symmetrical  reasons  because  the  contribution  of  these  two  backscattering  coefficients 
at  0  is  cancelled  by  the  one  at  ISO”  -  o  Also  in  Fig.  S  and  in  Fig  9.  experimental 
polanmetric  radar  data  obtained  by  MIT  Lincoln  Laboratory  are  displayed  Il4].  The 
frequency  of  operation  was  35  GHz,  and  the  angle  of  incidence  was  82”  (depression 
angle  of  6”).  The  radar  illuminated  on  a  vegetation  field  consisting  of  either  grass  or 
trees  at  a  range  of  2  km  A  good  correspondence  between  the  theoretical  predictions 
and  the  experimeutal  results  is  obtained. 
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£,^  =  (3.71  +i  0.134  )eo 
V  =  20*  v,  =  20’  SEA  ICE 
6  =  0.269  5^  =  0.528 

Ip  =  0.7  mm  =  1.1  mm 
- -  z  =  - 1 .65  m 

e2  =  (45  +  i40)eo 


Covariance  matrix: 


10'3- 


8.43 

0.23-11.18 
9.40  +  i  0.76 


0.23 +  i  1.18 
0.39 

0.02 +  i  0.134 


9.40  -  i  0.76 
0.02 -i  1.34 
10.6 


Figure  7:  Co>-ariance  matrix  for  the  backscattering  coefficients  at  $o,  =  40°  for 
sea  ice. 
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e,  =  (1.0233 +  i  0.0013  )eo 
=  (1.0233 +  i  0.02  )£o 
V  =  25“  V,  =  25’  TREES 

5  =  0.0451  5^  =  0.9644 

Ip  =  0.5  mm  1^  =  3.5  mm 

-  z  =  -10m 
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0.12 

1.2 
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0.01 

Born 

-10.6 

0.12 
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Figure  9:  Covariance  matnx  elements  for  tree  region  using  the  average  over  the 
azimuthal  angle  o  with  ^0  =  1'  Experimental  measurements  and 
theoretical  calculations  are  shown. 
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b)  Born  second-order  approximation 

Unlike  the  isotropic  Bom  first-order  approximation,  a  depolarization  factor 
oa,  ^  0  is  obtained  with  the  isotropic  Born  second-order  approximation  [9].  Fur¬ 
thermore.  it  was  shown  that  o***,  =  =  0.  Hence,  the  covariance  matrix  can  be 

aTitten  from  (11) 


Although  a  mathematical  proof  has  been  used  to  show  that  five  elements  of  the  covari¬ 
ance  matrix  are  zero,  this  fact  may  also  be  illustrated  from  a  physical  point  of  \iew. 
Let  us  concentrate  on  the  term  a„Ar  which  is  proportional  to  <  VV  ■  HV*  >  from 
(11)  By  assuming  that  a  vertically  polarized  electric  field  is  transmitted,  symbolized 
by  V  in  Fig.  10.  the  product  of  the  a-aves  polarized  in  the  horizontal  and  vertical  di¬ 
rections  at  the  receiver  is  taken.  The  aave  hits  a  first  scatterer  at  position  ri  and  then 
hits  a  second  scatterer  at  position  At  that  location,  the  electric  field  ^  induces  a 
dipole  which  can  be  decomposed  into  horizontal  and  vertical  components.  From  that 
point,  the  wave  is  backscattered  and  the  receiver  collects  either  HV'  or  V'l'.  Due 
to  symmetry  reasons,  the  second  scatterer  may  be  located  at  position  fl(.  The  corre¬ 
sponding  electric  field  'E  is  decomposed  into  a  vertical  component  and  a  horizontal 
component  which  is  in  the  opposite  direction  from  the  component  created  by  At 
the  receiver,  either  HV*  or  V'V  is  collected.  In  order  to  obtain  o,,*,.  the  product 
of  these  two  quantities  must  be  taken  Ho»-ever.  the  product  created  by  the  location 
at  r*,  will  have  an  opposite  sign  from  the  one  created  at  the  location  f'j.  By  taking 
the  average,  these  two  contributions  cancel  and  =  0  is  obtained.  Similarly  for 
every  point,  there  is  a  symnietnc  scatterer  which  will  cancel  its  contribution. 

The  depolarization  factor  is  proportional  to  <  HV  ■  HV'  >,  i.e..  a  vertically 
pwlarized  electric  field  is  transmitted  and  the  horizontal  component  of  the  scattered 
field  is  received.  Unlike  the  term  <  V’V'  •  HV'  >.  both  quantities  HV  and  HV'  due 
to  the  dipole  at  location  f'J  are  in  the  opposite  direction  from  the  ones  created  by  the 
scatterer  at  f'j.  This  implies  that  the  product  of  the  two  quantities,  which  have  the 
same  sign  independent  of  the  position  of  the  second  scatterer  f'j  or  Fj,  renders  the 
average  backscattenng  coefficient  <7^,  different  from  zero  as  shown  in  Fig.  11. 
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Figure  10:  Interpreiation  of  the  terra  <  HV’  >  with  the  Born  second  order 
approxiraation. 
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Figure  11:  Interpretation  of  the  terra  <  HV  ■  HV'  >  with  the  Bom  second  order 
approximation. 
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c)  Strong  permittivity  fluctuation  theory 

The  \'ariance  of  permittivity  fluctuations  is  usually  large  in  geophysical  media 
consisting  of  mixtures  of  constituents  »-ith  very  different  properties.  With  the  di¬ 
agrammatic  approach  [1],  the  propagation  and  scattering  in  random  media  can  be 
studied.  The  diagrammatic  approach  leads  to  the  Dyson's  equation  for  the  mean 
field.  It  should  be  noted  that  in  soKing  Dyson’s  equation  with  the  bilocal  approxi¬ 
mation.  both  the  observation  point  and  the  source  point  within  the  random  medium 
can  coincide  with  each  other  in  the  domain  of  integration.  To  take  care  of  the  singu¬ 
lar  nature  of  the  dyadic  Green’s  functions,  a  strong  permittivity  fluctuation  theory 
that  applies  to  both  small  and  large  variances  of  the  permittivity  functions  has  been 
developed  [11]. 

The  mean  field  and  the  mean  Green  s  function  are  ctilculated  by  replacing  the 
medium  with  equivalent  medium  having  an  effective  permittivity  which  accounts 
for  losses  due  to  both  absorption  and  scattering. 

It  should  be  emphasized  that  the  effect  of  the  random  medium  on  electromag¬ 
netic  waves  IS  better  described  with  the  strong  permittivity  fluctuation  theory  than  in 
the  case  of  the  traditional  “wealc  permittivity  fluctuations”.  For  the  case  of  a  random 
medium  described  by  a  mixture  of  two  components,  only  the  fractional  volume,  the 
correlation  lengths,  the  permittivity  of  the  two  components,  and  the  frequency  are 
necessary  to  describe  the  effect  of  the  earth  terrain.  The  other  parameters,  that  is 
the  variances  and  the  mean  permittivities,  needed  for  the  Bom  approximation  are 
computed  with  the  strong  permittivity  fluctuation  theory. 

In  order  to  illustrate  the  results  of  the  strong  permittivity  fluctuation  theory, 
we  compare  in  Fig  12  obtained  with  the  weak  permittivity  fluctuation  theory 
compared  with  the  strong  permittivity  fluctuation  theory.  We  use  a  correlation  func¬ 
tion  exponentially  decaying  in  all  directions  with  Ip  =  L  =  0.3  mm.  The  physical 
parameters  chosen  in  Fig  12  characterized  a  layer  of  dry  snow  [13l.  In  Fig.  12. 
we  plot  at,t,  in  the  function  of  the  frequency  for  ffo.  -  60*.  As  shown,  the  results 
obtained  with  the  strong  permittivity  fluctuation  theory  are  smaller  than  those  ob¬ 
tained  with  the  weak  permittivity  fluctuation  theory.  This  difference  gets  larger  wnth 
increasing  frequencies  This  is  due  to  the  fact  that  the  losses  due  to  scattenng.  which 
increase  with  increasing  frequencies,  are  taken  into  account  in  the  strong  pemiittivity 
fluctuation  theory. 


Figure  12;  Backscattenng  coefficient  cr*,),  in  function  of  the  frequency  computed 
with  the  weak  and  strong  peniiitti\ity  fluctuation  theories. 


